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A8STRAC 


,  BRIK^is—*;  unique  nuclear  exchange  model 


<NEM> . 


H-/if  the 


first  NEM  to  be  both  totally  interactive  and  to  rely  entirely 
upon  preemptive  goal  programming  to  drive  it*  allocations. 


-The  program 


iV^written  in 


Fortran  7Zya*dT-*ix i st*  on  the  'MX 


11-783  using  the  UNIX  operating  system. \ /t  includes  an  in- 
tergal  allocation  subrout  i  ne?  ^end^,  is  -therefore)  compl  etel  y 
transpor tabl e. 

ERIK  Mill  allocate  weapons  to  targets  in  order  to  meet 
damage  expectancy  <DE>  goals  of  the  user.  The  decision 
variable  is  a  non-integer  geometric  mean  representing  the 

number  of  weapons  allocated  to  an  entire  target  class,  using 

f>.  U'cj 

any  of\thr<*e<  objective  functions  and  up  to) tbirteen> different 

allocation  rules.  The  program  is  useful  for  a  wide  range  of 
scenarios  including  meeting  DE  goals  with  the  available 
arsenal  or  creating  an  arsenal  to  meet  established  goals. 

The  programmed  model  uses  a  preemptive  linear  goal 
programming  algorithm,  allowing  explicit  preferential  treat¬ 
ment  of  target  classes,  which  can  be  placed  into  any  of  seven  ->  7 
priorities.  Target  values  can  be  input  to  differentiate 

target  classes  within  the  same  priority.  BRIK  computes  single 

,  X>  1  A’sy 

shot  probability  of  survival  using}  Che5  nuclear  targeting 
vulnerability  XVN>^  system^  of  the  Defense  Intelligence  Agency 
or  from  vulnerability  expressed  as  a  PSI  overpressure  number. 

As  currently  dimensioned,  the  program  can  handle  up  to  20 
target  classes,  23  weapon  classes,  and  28  hedging  constraints. 
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^  BRIK  X  unique  nuclear  exchange  model  (NEM),  s  the  first 
NEM  to  be  both  totally  interactive  and  to  rely  entirely  upon  preemptive 
goal  programming  to  drive  its  allocations.  -The- program -is>  written  in 
FORTRAN  7^4e<F>exists  on  the  VAX  11  -780 fusing  the  UNIX  operating  system, 
/ft  Includes  an  internal  allocation  subroutine^and^ls  tberefore^completely 
transportable.'  *  r~f 

BR]!K  will  allocate  weapons  to  targets  In  order  to  meet  damage 
expectancy  (0E)  goals  of  the  user.  The  decision  variable  is  a  non-integer 
geometric  mean  representing  the  number  of  weapons  allocated  to  an  entire 
target  class,  using  any  of  -thre^objective  functions  and  up  to  -thirteen^ 
different  allocation  rules.  The  program  is  useful  for  a  wide  range  of 
scenarios,  including  meeting  0E  goals  with  the  available  arsenal  or 
creating  an  arsenal  to  meet  established  goals. 

The  programmed  model  uses  a  preemptive  linear  goal  programming 
algorithm,  allowing  explicit  preferential  treatment  of- target  classes, 
which  car  be  placed  Into  any  of  •eoven/fpriorltles.  Target  values  can  be 
<nput  to  differentiate  target  classes  within  the  same  priority.  BRIK  j 
computes  single-shot  probability  of  survival  using  fchernucl ear  targeting  ~ 
vulnerability  X'VM-ysystem  the  Defers#- Intelligence  Agenc^or  from 
vulnerability  expressed  as  a  PSI  overpressure  number.  As  currently 
dimensioned,  the  program  can  handle  up  to  20  target  classes,  20  weapon 
classes,  and  20  hedging  constraints. 
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The  report  describes  the  physics]  end  mathematical  model , 
scenario  limitations,  input  requirements,  and  the  damage 
function.  It  includes  a  user  guido,  variable  and  subroutine 
definitions,  and  a  computer  listing. 
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The  existence  of  nuclear  arsenals  has  progressed  -from 
the  production  of  the  fir  it  atomic  bomb  to  the  development  of 
enhanced  "traditional"  triad  elements  such  as  the  B-1B  bomb¬ 
er,  the  MX  missile,  the  SS-18,  and  the  Backfire  bomber.  This 
activity  has  resulted  in  the  deployment  of  a  substantial 
number  of  delivery  systems  and  nuclear  weapons.  Also,  large 
expenditures  of  funds  on  research  and  development  have  pro¬ 
duced  systems  which  for  one  reason  or  another  never  reached 
deployment. 

When  nuclear  arsenals  began  to  grow,  analysts  attempted 
to  determine  the  "Nuclear  Balance."  In  the  early  l?50's,  the 
US  used  straight  inventory  comparisons  (Ref  13) .  With  time, 
these  static  measurements  proved  inadequate  and  failed  to 
give  a  true  indication  of  a  force's  capability. 

When  Robert  McNamara  became  Secretary  of  Defense  in 
1961,  he  introduced  "Defense  Economics"  (Ref  3:1.5).  The 
ingredients  of  this  dynamic  approach  to  nuclear  force  plan¬ 
ning  required  (1)  identifying  the  force's  mission  objectives 
and  (2)  evaluating  how  well  alternative  forces  could  meet 
these  objectives.  Thus,  force  assessment  techniques  changed 
from  static  to  dyn«mic  measurement. 

Today,  analysts  and  policy  makers  are  still  using  a 
dynamic  approach  to  assess  the  adequacy  of  the  US  strategic 
arsenal  to  deter  the  Soviets  and  to  assess  the  effect  of  new 
US  and  Soviet  weapon  systems  on  the  strategic  balance.  To 


assist  Analysis,  a  number  of  nuclear  exchange  models  (NEMs) 
have  been  built. 


These  models  generally  fall  into  two  broad  categories  — 
computer  simulation  models  and  analytical  models.  Computer 
simulation  models  are  used  to  establish  targeting  strategies 
for  various  contingencies  (Ref  27il).  These  models  consider 
many  variables.  Consequently,  they  require  many  hours  of 
computer  time  and  are  generally  infleMble  regarding  varying 
applications.  Unlike  simulation  models,  analytical  models 
use  expected  value  concepts.  They  are  used  within  the  000  to 
study  the  issue  of  optimizing  weapons  allocation. 

Analytical  allocation  models  found  in  the  literature 
differ  considerably  in  degrees  of  realism  and  flexibility. 
Realism  refers  to  the  model's  output  reflecting  an  actual 
representation  of  what  is  possible.  An  allocation  model  is 
realistic  if  the  assignment  of  a  weapon  to  a  target  is  in 
fact  feasible.  An  example  of  an  infeasible  assignment  is  an 
allocation  of  a  weapon  to  a  target  that  is  beyond  the 
weapon's  range.  Flexibility  refers  to  the  model's  ability  to 
handle  many  different  types  of  strategic  allocation  problems. 
This  is  achieved  by  allowing  a  variety  of  allocation  rules 
and  objectives. 

Problem  Statement 

A  review  of  the  analytical  allocation  models  indicated 
two  things.  First,  if  they  have  the  flexibility  to  handle 
many  different  scenarios,  they  are  large  and  their  design 
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■favors  batching,  thus  prohibiting  interactive  analysis. 
Second,  i-f  they  are  of  reasonable  size  to  permit  interactive 
analytics,  they  are  limited  in  scope  and  capabilities. 


The  objective  of  this  research  effort  was  to  build  an 
allocation  model  that  was  both  realistic  and  flexible,  yet  of 
reasonable  size  to  permit  interactive  analytics.  Milan 
Zelany  says  • . . .researchers  seem  to  be  evolving  a  consensus: 
let  the  human  being,  human  decision  maker,  become  a  core 
around  which  to  build  our  techniques,  adapting  them  to  his 
needs  and  amplifying  his  own  decision  making  powers"  (Ref  28:2). 
Nuclear  exchange  models,  with  their  batch  mode  designs,  prohi¬ 
bit  any  convenient  human  interaction  with  the  computer.  There¬ 
fore,  it  was  the  intent  of  this  effort  to  develop  a  model 
that  provided  the  qualities  of  the  larger  models  yet  allowed 
interaction  between  the  analyst  and  the  machine. 

A  model  of  this  nature  should  not  be  restricted  to  a 
single  capability.  Therefore,  building  this  model,  BRIK,  to 
perforfn  just  one  type  of  analysis  would  have  defeated  the 
purpose  of  this  effort.  An  analyst  who  would  potentially  use 
an  NEM  should  have  the  ability  to  perform  a  variety  of  tasks. 
Therefore,  the  capability  of  this  model  exceeds  the  require¬ 
ments  of  any  one  single  analysis.  Given  a  weapon  base  and  a 
target  base,  the  model  performs  an  allocation  based  upon  a 
variety  of  objectives.  For  example,  BRIK  includes  the  fol- 
1  owing: 
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1)  It  optimally  allocates  weapons  to  meet  targeting 

goals. 

2)  If  goals  are  under-determi ned ,  i.e.,  if  the  goals  are 
met  and  there  are  weapons  remaining,  the  model  could  minimize 
the  number  of  warheads,  megatonnage,  equivalent  megatonnage 
<EMT> ,  or  cou.n  term!  1  i  tary  potential  <CMP>. 

3)  It  determines  what  -force  is  necessary  to  meet  damage 
expectancy  <DE>  goals. 

To  install  realism  in  the  allocation  process,  numerous 
allocation  rules  are  allowed.  For  example,  the  model  per¬ 
mits  the  -following: 

1)  Designation  o-f  inappropriate  weapon/target  assign¬ 
ments. 

2)  Attaining  minimum  levels  o-f  DE  on  designated  target 
cl  asses. 

3)  Attaining  minimum  levels  of  DE  on  designated  target 
classes  using  a  specific  set  of  weapons. 

4)  Enforcing  an  upper  level  of  damage  on  designated 
target  classes  resulting  from  a  specific  set  of  weapons. 

5)  Restricting  the  number  of  a  class  of  weapons  which  can 
be  allocated  to  a  designated  target  class. 

6)  Restricting  the  number  of  weapons  which  can  be  allo¬ 
cated  to  each  target  in  a  particular  target  class. 

7)  Building  any  user-defined  allocation  rule. 

To  increase  the  interactive  capabilities  of  the  model, 
BRIK  allows  the  user  to  ask  the  following  "what  if"  type 
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questions  while  seated  at  a  terminal! 

1)  What  if  the  number  of  available  weapons  changed? 

2)  What  if  the  DE  goals  changed? 

3)  What  if  weapon  characteristics  changed? 

4)  What  if  target  characteristics  changed? 

This  concludes  the  description  of  the  objective  and  scope 
of  this  thesis.  The  next  section  describes  the  approach 
followed  by  the  authors  in  producing  this  model. 

Approach 

This  effort  was  a  continuation  of  a  thesis  completed  in 
December  1982,  by  Michael  C.  Wambsganss  (Ref  27) .  He  demon¬ 
strated  that  auxiliary  variables  in  linear  programming  con¬ 
straints  could  be  used  to  allocate  weapons  to  targets  ?n 
order  to  meet  user-defined  damage  expectancy  goals.  His  meth¬ 
odology  provides  the  foundation  for  this  effort  and  will  be 
reviewed  in  the  literature  review  section. 

After  studying  Wambsganss'  work,  it  was  necessary  to 
design  BRIK.  To  increase  transpor tabi 1 i ty ,  the  authors  de¬ 
cided  that  BRIK  should  not  use  a  library  package  for  its 
allocation  routine.  Therefore,  a  search  was  conducted  for  a 
suitable  linear  programming  package.  The  authors  finally 
decided  to  use  PAGP,  a  goat  programming  algorithm  acquired  from 
The  Association  for  Computing  Machinery.  Next,  time  was 
spent  determining  what  allocation  rules  should  be  included 
and  how  the  objective  functions  should  be  designed.  With  the 
basic  model  adequately  outlined,  the  coding  stage  began. 
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Aft  QRIK  was  being  written,  it*  design  -followed  an  evolu¬ 
tion  process.  While  its  basic  outline  is  unchanged,  some 
notable  -features  were  included  well  into  the  project.  The 
two  most  noteworthy  additions  are  the  VNTK  vulnerability 
system  and  the  -file  systems  -for  the  weapon  and  target  bases. 
The  addition  o-f  the  VNTK  system  proved  to  be  invaluable  as  i  t 
permitted  direct  comparison  with  the  Arsenal  Exchange  Model 
(AEM)  during  the  -final  portion  o-f  this  e-f-fort. 

Format 

Chapter  II  provides  a  description  o-f  the  elements  -found 
in  most  nuclear  exchange  model  si  target  base,  weapon  base, 
and  damage  -function.  This  is  -followed  by  a  review  o-f  many  o-f 
the  nuclear  exchange  models  -found  in  the  literature. 

Chapter  III  presents  some  characteristics  o-f  BRIK.  These 
include  model  restrictions,  scenario  limitations,  and  input 
requirements.  Also,  in  this  chapter  is  an  explanation  of  the 
single  shot  probability  of  survival  computations  and  a  discus¬ 
sion  of  the  damage  function. 

Chapter  IV  derives  all  of  the  constraints  used  by  BRIK  in 
the  linear  programming  matrix. 

Chapter  V  reviews  Goal  Programming.  This  is  followed  by 
a  description  of  the  mathematical  formulation  of  BRIK,  includ¬ 
ing  the  objective  functions  and  the  matrices.  Finally,  the 
assignment  algorithm,  PAGP,  is  reviewed. 

Chapter  VI  concerns  both  verification  and  demonstration. 
After  a  discussion  concerning  BRiK's  "correctness* ,  a 
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comparison  is  mads  between  the  solutions  from  a  sample  prob¬ 
lem  as  computed  by  6RIK  and  AEM.  This  comparison  is  followed 
by  numerical  examples  demonstrating  BRIK's  versatility. 

Chapter  VII  provides  a  summary  of  the  modeling  process, 
descriptions  of  potential  uses  for  BRIK,  and  a  discussion  of 
potential  improvements  to  the  model . 


7 


II  Literature  Review 


The  term  "model*  is  used  broadly  in  the  literature  with 
various  definitions.  In  this  study,  a  model  is  defined  as  a 
computer  program  which  converts  information  about  the  targets, 
weapons,  and  employment  plans  of  two  sides  into  a  prediction 
of  the  outcome  of  a  war.  This  is  commonly  called  an  arsenal 
exchange.  For  those  readers  who  are  not  familiar  with  nuclear 
exchange  models,  a  general  overview  of  the  model  components 
and  terminology  follows  along  with  a  discussion  of  several 
models. 

MWEl--ELEM.afl:g 

Nuclear  exchange  models  can  be  grouped  into  two  general 
categories!  computer  simulation  models  and  analytical  models 
<Ref  27t  1)  . 

Simulation  Models.  Computer  simulation  models  are  used  to 
develop  targeting  strategies.  These  models  require  extensive 
data  preparation,  contain  a  large  number  of  variables,  and 
often  use  hours  of  computer  time.  The  complex  structure  of 
the  model  frequently  precludes  sensitivity  analysis  because  of 
the  large  amount  of  runs  necessary  to  obtain  an  expected  value 
which  has  a  high  degree  of  certainty.  Simulation  models  are 
also  regarded  as  generally  inflexible  and  oriented  to  a 
specific  application  (Ref  27:1-2). 
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Analytical  modal  *  cantor  on  expected 


values  such  as  an  axpactad  damage  leva!  and  axpactad  surviving 
arsanal .  Data  praparation  timas,  computer  run  times,  and 
levels  of  complexity  vary  widely,  increasing  in  proportion  to 
the  level  of  realism  required  by  the  specific  model. 

Analytical  models  also  permit  sensitivity  analysis  of 
variabl as. 

The  analytical  model  generally  attempts  to  optimize  either 
a  weapon  allocation  or  cost.  To  optimize  weapon  allocation, 
the  model  searches  for  the  highest  level  of  expected  damage  to 
the  opponent  from  a  specified  weapon  arsenal .  To  optimize 
cost,  the  model  seeks  the  least  expensive  combination  of 
weapons  which  causes  a  specified  level  of  damage  to  the 
opponent.  Analytical  models  are  the  most  widely  used  nuclear 
exchange  models  (Ref  27i2) . 

Although  models  differ  widely  in  construction  and 
application,  both  simulation  and  analytical  models  share 
common  submodel  si  the  weapon  complex,  the  target  complex,  the 
engagement  and  allocation  rules,  the  damage  function,  and  the 
algorithm  or  solution  technique.  The  weapon  and  target 
submodels  generally  determine  the  complexity  of  the  other 
submodels  (Ref  27t9>.  Since  model  complexity  may  be  an  issue 
to  prospective  model  users,  the  content  of  the  weapon  and 
target  complex  submodels  is  given  in  more  detail,  followed  by 
a  discussion  of  terminology  that  is  used  in  nuclear  exchange 
literature. 


Weapon  Comp!  ex  .  The  weapon  compl  ex  addresses  three 
distinct  categories  (Ref  27:10):  the  scope,  the  number  of 
weapon  system  types  and  their  associated  penetration  aids} 
weapons  reach,  which  refers  to  the  ability  of  a  weapon  to 
reach  a  specific  target}  and  commitment  policy  of  each  weapon, 
which  refers  to  tho  amount  of  strikes  launched,  quality  of  the 
intelligence  and  bomb  damage  assessment,  and  weapon 
availability  uncertainties. 

a)  Scope:  The  scope  may  specify  single  or  multiple  weapon 
types  based  on  their  weapon  characteristics  such  as  circular 
er.' or  probable  <CEP)  ,  the  weapon's  explosive  yield,  weapon 
system  reliability,  and  defense  penetration  aids.  CEP  is  the 
radius  of  an  imaginary  circle  centered  on  the  target  in  which 
at  least  30/C  of  the  weapons  aimed  are  expected  to  land.  The 
weapon's  explosive  yield  is  the  size  of  the  weapon  usually 
measured  in  kiloton  or  megaton  equivalent  TNT  yields.  Weapon 
system  reliability  <HSR)  is  the  probabi 1 i ty  that  once  the 
weapon  system  is  launched  it  will  function  properly  during  its 
mission.  Penetration  aids  are  those  items  that  help  a  weapon 
system  defeat  a  defen  ive  system.  Chaff  and  decoys  are 
examples  of  penetration  aids.  In  many  models,  a  probability 
of  penetration  <PTP)  simulates  enemy  defenses.  Sometimes,  in 
order  to  simplify  the  input  requirements,  models  will  combine 
FTP  and  WSR  into  a  single  probability.  Models  either  assume 
that  all  weapons  are  available  {deterministic)  or  that  there 
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is  a  certain  amount  of  system  unreliability  in  launch 
{probabilistic). 

b)  Weapon  Reach i  Weapon  reach  may  limit  feasible  assign¬ 
ments  of  weapons  to  targets  or  degrade  weapon  accuracy  as  the 
range  traveled  increases.  A  0-1  matrix  can  normally  be  used 
to  represent  weapon-target  assignments  with  a  1  meaning  a 
feasible  assignment.  The  targets  that  a  MIRV  (multiple 
independently  targeted  re-entry  vehicles)  can  attack  are  not 
only  limited  by  the  range  of  the  carrier  but  also  by  the  MIRV 
footprint.  Most  models  do  not  consider  the  MIRV  footprint. 

Instead,  all  weapons  are  considered  independent.  This  assump¬ 
tion  may  provide  allocations  that  are  beyond  the  range  of  the 
actual  weapon  capability. 

Range  is  not  the  only  factor  considered  in  weapon  reach. 

The  BSWAOE  model  uses  time  as  a  factor  to  limit  feasible 
assignments.  If  a  target  is  time  urgent,  meaning  that  the 
target  must  be  hit  by  weapons  that  can  strike  the  target 
within  a  specified  time  limit,  an  allocation  will  be  allowed 
only  from  those  weapons  that  can  arrive  at  the  target  in  the 
specified  time  period. 

c)  Weapon  Commitment  Strategies!  Models  usually  assume 
that  the  weapons  will  be  allocated  against  existing  targets. 

For  example,  only  those  silos  that  contain  missiles  will  be 
attacked.  This  aspect  indirectly  assumes  that  the  attacker 
has  prior  intelligence  about  his  adversary's  battle  plan. 

Even  though  most  assume  some  prior  knowledge  about  an 
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adversary,  most  modal s  do  not.  consider  any  battle  damage 
assessment  <8DA) .  The  models  will  continue  to  allocate 
weapons  against  a  target  until  its  damage  expectancy  has  been 
attained.  There  is  no  checking  procedure  to  indicate  whether 
a  target  has  been  destroyed  by  an  earlier  arriving  warhead,  so 
that  the  next  warhead  can  be  allocated  against  a  target  that 
has  not  been  attacked  or  destroyed. 

A  model  can  also  be  categorized  on  the  basis  of  how  many 
strikes  it  can  compute.  Models  limited  to  a  single  strike  are 
usually  limited  to  evaluating  a  single  goal  by  a  single 
attacker.  A  two-strike  model  can  corttpute  the  results  o-f  an 

i 

initial  attack  -followed  by  a  retaliatory  strike  without  making 
two  runs.  Models  that  can  handle  three  or  more  strikes  can 
evaluate  the  optimum  reserve  -force  that  could  be  held  back 

I 

-from  a  strike  to  discourage  an  opponent  from  a  retaliatory 

strike.  j 

. 

■ 

Taroet  Complex.  The  target  complex  specifies  the  type  of 

!  ■■ 

target,  the  value  of  the  target,  and!  the  target  defenses. 

I 

a)  Target  Type*  Targets  are  classified  as  either  point, 
area,  or  collateral  targets.  The  first  two  classes  are  the 
most  commonly  used.  A  target  is  considered  a  point  target 
when  a  single  well -placed  weapon  can  kill  it,  an  example  being 
a  missile  silo.  If  more  than  one  weapon  is  required  to  cover 
the  target's  area,  such  as  a  military  base,  it  is  considered 
an  area  target. 

Collateral  targets,  however,  are  a  collection  of  point 
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targets  or  an  area  target  that  has  been  broken  up  into 
individual  aim-points.  In  contrast  to  a  collection  of  point 
targets  considered  to  be  independent  in  that  no  single  weapon 
can  kill  more  than  one  target,  collateral  targets  are  targets 
located  close  enough  together  that  all  of  them  can  be 
destroyed  by  a  single  weapon.  Including  collateral  targets 
significantly  increases  the  complexity  of  a  model  because  the 
point  targets  that  are  collateral  targets  are  not  statis¬ 
tically  independent. 

Targets  can  also  be  classified  according  to  target 
characteristics.  Two  commonly  usee  classes  of  targets  are 
Value  and  Force  type  targets.  Force  targets  are  those  tar¬ 
gets  which  if  not  destroyed  have  the  ability  to  do  immediate 

» 

damage  to  the  opponent.  These  include  such  targets  as  ICBMs 
and  bomber  bases.  Value  targets  do  not  threaten  the  opponent 
but  do  have  value  to  the  attacker.  These  include  targets  such 
as  factories  and  cities. 

b)  Target  Valuex  The  usual  measure  of  effectiveness  with 
which  different  weapon  allocations  are  compared  is  the 
expected  target  value  killed.  The  more  complex  the  calcu¬ 
lation  of  target  value,  the  more  complex  the  model.  The 
simpler  models  either  do  not  address  target  value  or  assume 
that  all  targets  are  given  the  same  value  of  one  unit.  More 
complex  models  prioritize  the  targets  as  to  their  importance 
of  destruction.  The  most  complex  models  assign  numerical 
values  to  the  targets;  however,  there  is  a  serious  lack  of 
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methodology  in  generating  target  values  (Ref  27:13). 


c)  Target  Defenses:  The  simplest  models  either  assume  no 
de-fenses  or  include  the  defenses  in  the  form  of  penetration 
probability.  More  advanced  models  treat  defenses  as  a 
specified  attrition  of  attacking  weapons  before  a  target  can 
be  damaged.  In  many  models  the  defenses  are  simulated  by  a 
specific  probability  of  penetration  for  each  target  complex 
versus  a  particular  type  of  weapon  (Ref  27:14). 

Model  Terminology.  The  following  terminology  is  used 
extensively  in  the  literature. 

a)  Damage  Expectancy  <DE) i  DE  is  the  statistical 
expected  value  that  the  weapon  will  arrive,  detonate,  and 
cause  a  certain  level  of  damage  (Ref  24:11). 

b)  Hedge:  A  hedge  is  an  input  rule  or  constraint  that 
allows  the  analyst  to  specify  auxiliary  side  goals  or  extra 
requirements  while  maintaining  the  ability  to  allocate  in 
accordance  with  the  main  objective. 

c)  Force  Optimization:  Force  optimization  is  concerned 
with  the  determination  of  the  appropriate  force  structure  to 
accomplish  a  set  of  specified  goals  subject  to  force 
availability,  flexibility,  and  budget  constraints. 

d)  Target  Optimization:  Target  optimization  consists  of 
allocating  weapons  against  targets  in  such  a  way  as  to  obtain 
an  optimal  employment  of  the  weapon  mix.  Normally  this 
involves  an  allocation  that  attempts  to  maximize  the  damage 
achieved  to  a  set  of  targets. 
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REVIEN  OF  EXISTING  MODELS 


The  following  list  of  analytical  models  and  various 
mathematical  formulations  were  examined  to  determine  the 
essential  elements  of  various  models.  This  examination  helped 
in  the  conceptual iztion  and  formulation  of  functions  and 
formulas  used  in  BRIK. 

Arsenal  Exchange  Model .  The  Arsenal  Exchange  Model  <AEM> , 
developed  by  Martin-Mariettta  Corporation,  is  the  most  widely 
used  strategic  analysis  model.  It  was  initially  built  in  1964 
as  a  cost  effectiveness  model  for  use  in  advanced  ICBM 
studies.  But,  in  response  to  user  requirements,  it  has 
evolved- into  a  large-  expected  value  model  used  to  study  the 
structure  of  strategic  forces.  As  of  the  last  revision, 
documentation  published  in  March,  1982,  the  AEM  requires  1500k 
of  computer  storage. 

AEM  enjoys  tremendous  flexibility  in  its  allocation  rules. 
It  permits  analysis  of  different  types  and  levels  of  exchanges 
and  is  constructed  to  accept  a  variety  of  strike  objectives. 
The  analyst  also  has  the  ability  to  structure  a  variety  of 
employment  constraints  input  as  hedges.  These  hedges  fall 
into  three  separate  categories  <Ref  3j10)s 

1)  A  specified  amount  of  value  destroyed  on  a  particular 
target  set  must  be  accomplished  by  a  particular 
weapon  set . 

2)  A  specified  amount  of  weapons  must  be  used  against  a 
particular  target  class. 
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3>  A  target  set  can  not  be  attacked  by  more  than  a 
certain  amount  of  a  particular  set  of  weapons. 

The  model  employs  a  singleshot  kill  probability  against 
point  targets  as 

<Ri  i/CEP,  ) 

SSPK  -  I  -  .3  A 


where  Rij  1  s  the  lethal  radius  of  weapon  i  when  it  is 
allocated  against  target  j,  CEP^  is  the  circular  error 
probable  of  weapon  i.  For  area  targets,  AEM  adjusts  the  point 
target  SSPK  as  follows! 


SSPRarea  a  SSPk||(1-(1-§)B)+(1-|) (|) (1/B) 

where  X,  B,  and  C  are  shape  parameters  resulting  from  a 
regression  by  RAND  (Ref  20x33).  Area  and  terminal  defenses 
are  also  incorporated  in  the  model. 

The  optimal  allocation  was  attained  by  a  goal -oriented 
Integer  Linear  Program  using  the  Lagrangian  method.  However, 

because  users  wanted  to  trade  off  integer  solution  assurance 

1 

| 

for  greater  computational  efficiency,  the  model  has  been! 

i 

converted  to  a  mul ti -objective  model  solved  by  a  linear  ! 
programming  <LP>  method  of  allocation.  \ 

B SHADE.  The  Strategic  Weapon  Allocation  Damage 
Expectation  (SHADE)  model  was  developed  as  an  "in-house" 
project  by  Headquarters  SAC/XPS  in  the  mid-l?70's  and  was 
designed  to  have  a  fast  run  time  (Ref  12s35) .  Consequently, 
it  aggregates  and  clusters  much  of  the  targeting  data  so  that 
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data  manipulation  is  kept  to  a  minimum.  In  its  current  -form, 
there  are  two  separate  programs,  BSUDS  and  BALLOT. 

The  Single  Uni-form  Data  System  <SUDS>  is  a  user-friendly 
program  that  allows  an  analyst  to  enter  data  from  a  card  deck 
<Ref  1).  It  takes  all  of  the  pre-spec  if ied  force  structures 
and  allocation  rules  and  converts  them  into  a  form  which  is 
capable  of  being  used  in  the  allocation  phase.  The  model  uses 
the  WTK  target  vulnerability  system  to  determine  its  SSPKs. 

The  allocation  phase  uses  a  computer  program  called 
BALLOT,  a  version  of  an  earlier  program  called  ALLOT  <Ref  23). 
it  uses  linear  programming  techniques  to  determine  the 
optimal  laydown  of  a  nuclear  arsenal  against  an  aggregated 
target  base,  where  optimal  is  defined  as  attaining  the 
maximum  DE  of  the  entire  target  complex.  It  is  a  one-si  Jed 
exchange  where  the  optimal  weapon-to-target  assignment  can  be 
found  according  to  one  of  the  following  objectives: 

1)  The  total  value  destroyed  is  to  be  a  maximum. 

2)  The  total  number  of  weapons  employed  is  to  be  a 
minimum. 

3>  The  total  megatonage  employed  is  to  be  a  minimum. 

4)  The  total  EMT  (Equivalent  Megatonnage)  employed  is 
to  be  a  minimum. 

5)  The  total  CMP  (Countermilitary  Potential)  is  to  be  a 
minimum. 

BALLOT  has  the  ability  to  accept  manual  1 y-inser ted 
constraining  equations  into  the  LP  (Ref  23:8)  and,  as  was 
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mentioned  before,  it  has  the  ability  to  declare  certain 
target  classes  as  being  time  urgent. 

Hodson .  This  paper  describes  an  approach  to  weapon 
allocation.  The  authors  point  out  that  the  general  problem 
of  weapon  allocation  is  that  "of  assigning  weapons  to  targets 
in  such  a  way  as  to  achieve  a  desired  level  of  destruction". 
The  authors  go  on  to  say  that  "when  different  combinations  of 
weapon  types  and  numbers  are  applied  to  the  same  target 
system  with  the  same  targeting  philosophy,  it  is  possible  to 
measure  the  effect  of  adding,  retaining,  or  deleting  certain 
weapon  systems"  <Ref  17i2> .  The  authors  use  linear  pro¬ 
gramming  to  determine  an  optimal  allocation.  If  Sjj 
represents  the  single  shot  probability  of  survival  <SSPS>  of 
a  type  j  target  assigned  a  type  i  weapon,  the  expected 
survival  probability  of  that  target  would  be 


where  is  the  number  of  type  i  weapons  assigned  to  the 
type  J  target.  Hodson  linearizes  this  survival  function  so 
that  LP  can  be  used.  By  using  the  monotonic  characteristics 
of  the  logarithmic  function, 

logcTT  <S,  .  i;*  >3 

i  13 


becomes 


ZC  <  1  og  SjjJXjj]. 


This  model  assumes  that  SSPS's  are  known  and  weapon  and 


target  characteristics  are  not  discussed.  ^Iso,  since  no 
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attempt  is  made  to  keep  X^  an  integer,  this  allocation 
method  is  only  an  approximation.  However,  the  authors 
present  analysis  that  concludes  the  error  from  this 
approximation  is  small. 

Hambsoanss .  In  an  Air  Force  Institute  of  Technology 
thesis,  Michael  C.  Mambsganss  developed  a  framework  that 
demonstrated  that  a  realistic  and  flexible  model  could  be 
built  and  still  be  manageable  in  terms  of  limited  computer 
resources.  The  model  could  .  duress  both  weapon  performance 
and  cost  issues. 

f 

This  framework  was  developed  by  extracting  features  o-f 

various  nuclear  exchange  models  in  existence.  The  formulae 

i 

to  compute  the  single-shot  probabilities  of  survival  <SSPS> 
for  various  weapon/target  combinations  was  acquired  from  AEM. 

i 

To  incorporate  fast  run  time,  a  linear  programming  technique 
was  used. |  To  linearize  the  damage  function 


where  is  the  single-shot  probability  of  survival  for 
target  j  when  attacked  by  weapon  i,  and  Xjj  is  the  amount  of 
i  weapons  allocated  to  target  J,  Hodson's  logarithmic 
function  was  used. 

Mambsganss  introduced  the  idea  of  linear  goal  programming 
to  allow  consideration  of  multiple  objectives.  The  goal 
programming  methods  used  were  weighted  goal  programming  for 
the  weapon  allocation  and  compromise  programming  when 
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investigating  cost  issues. 

The  objective  -function  -for  the  allocation  problem  is 

Minimize  Z  M  ^  M^ 

where  Wj  is  the  weight  or  value  associated  with  goal  jt  and  Mj 
is  the  variable  representing  the  underachievement  o-f  goal  j 
<Re-f  27:79).  The  model  Mambsganss  developed  required 
considerable  amount  o-f  manual  input  o-f  data. 

This  model  and  the  three  models  previously  discussed  were 
the  primary  sources  o-f  information  and  techniques  -for  the 
authors.  The  -following  models,  although  not  used  as 
extensively  in  this  research,  are  included  -for  the  reader's 
bene-f  it. 

SROTTE.  This  is  a  two-strike  model  in  which  Red  employs 
all  its  weapons  against  both  Blue's  force  and  value  targets, 
and  Blue  retaliates  with  its  surviving  arsenal  against  Red's 
value  targets.  The  model  can  handle  four  weapon  types  for 
both  Red  and  Blue.  SSPKs  for  each  side  are  expressed  as 
functions  of  reliability,  penetration  probability,  and 
single-shot  probability  of  kill  with  the  damage  function 
assumed  to  have  an  exponential  nature.  The  model  optimizes 
the  damage  difference  between  Blue's  and  Red's  strikes  using 
Fauld's  branch  and  bound  algorithm  as  a  piece-wise  linear 
approximation  to  the  damage  difference  (Ref  15) . 

COPE  58,  The  is  a  widely  used  aggregated  effectiveness 
model  developed  by  Lambda  Corporation  which  can  handle  a 
mixture  of  weapon  types,  permits  up  to  three  strikes  and  is 
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restricted  to  48  target  types.  The  model  considers  area  and 
terminal  defenses  and  contains  submodels  to  determine  missile 
and  bomber  penetration  probabilities  as  well  as  weapon/target 
allocation.  Because  of  the  complexity  of  the  damage  function 
and  the  nonlinearity  of  portions  of  the  constraint  set,  a 
Lagrangian  solution  method  is  used  <Ref  13). 

DAY.  This  model  splits  the  general  allocation  problem 
into  two  problems,  email  and  large.  The  small  problem  is  to 
allocate  weapons  within  target  complexes,  while  the  large 
problem  coordinates  the  allocations  and  assigns  an  optimal 
desired  ground  zero  <DG2)  and  height  of  burst  (HOB)  to  each 
weapon  in  the  weapon  stockpile.  The  model  is  complex  and 
requires  the  solution  of  a  nonlinear  programming  problem  for 
every  target  complex.  DAY  assigns  values  for  each  target  in 
the  complex  (Ref  11). 

SOLAND.  This  model  examines  discrete  ABM  levels  so  as  to 
minimize  any  damage  inflicted  from  an  optimal  attack  by  an 
opponent.  Both  point  and  area  defenses  are  considered  and 
are  assumed  to  be  in  the  exhaustive  mode,  which  implies  that 
no  damage  can  be  inflicted  on  a  target  until  the  defenses 
have  been  exhausted.  The  attacker  is  assumed  to  know  the  ABM 
levels  and  a  min-max  approach  is  used,  which  is  to  minimize 
the  weapons  necessary  to  exhaust  the  defenses  and  to  inflict 
as  much  damage  as  possible.  There  are  bounds  which  specify 
the  number  of  missiles  that  can  be  allocated  against  any 
particular  target  (Ref  23). 
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MANNE .  This  paper  assumes  that  all  weapons  are  homo¬ 
geneous.  The  model  uses  X* ^  as  the  probabi 1 i ty  with  which 
weapon  i  is  allocated  against  target  j  and  is  the 
conditional  probability  of  kill  for  weapon  i  against  the 
target  j.  The  model  minimizes 

£  V,  IT  (1-P,  ,  X-  .) 

3  J  {  iJ 

where  V.  is  the  worth  of  target  j.  Manne  formulates  a  linear 
J 

approximation  of  the  above  equation  and  solves  the  allocation 
using  a  transportation  algorithm  (Ref  19). 

ARMS.  This  paper  considers  both  a  cost  optimization 
model  and  a  single-sided  allocation  model.  Damage  functions 
are  determined  via  analytical  fatality  curves  taken  from  the 
Code  30  model.  The  cost  optimization  model  contains  the 
allocation  model  as  a  submodel  and  considers  only  annual 
costs  and  production  costs.  Development  costs  are  incor¬ 
porated  implicitly  with  these  costs  by  assuming  linear  build¬ 
up  of  weapons  over  a  specified  time  period.  Both  models 
consider  100  variables  and  are  solved  by  the  sequential- 
unconstrained-minimization  technique  (SISTF)  (Ref  27»  19)  . 

BRACKEN.  This  is  a  convex  programming  model  which  opti¬ 
mizes  an  SLBM  attack  on  bomber  bases.  The  model  is  split 
into  two  stages.  The  first  stage  allocates  the  submarines  to 
a  set  of  feasible  launch  areas  and  the  second  stage  is  the 
actual  allocation  of  SLBMs  to  bomber  bases  to  maximize 
damage.  Not  all  missiles  are  allowed  to  launch  simultan¬ 
eously,  and  some  bombers  are  allowed  to  scramble  safely  as  a 
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function  of  th*  missile's  flight  tint*  and  rang*.  A 
standard,  non I  inoar  programming  rout  in*  was  *mploy*d  to  yield 
th*  global  maximum  <R*f  7). 

MIERCORT.  This  is  a  two-sided  nuclear  *xchang*  mod*] 
which  allocates  homogeneous  missiles  against  cities  that  are 
defended  by  area  and  point  defenses.  The  model  operates  in 
the  exhaustive  mode,  allocating  missiles  to  the  defenses 
prior  to  attacking  th*  cities.  Th*  allocation  is  achieved 
through  a  nonlinear  integer  progamming  algorithm  using  a 
branch  and  bound  technique  (Ref  21) . 

DAVID.  This  model  allocates  a  non homogeneous  weapon  base 
against  an  undefended  target  complex.  Th*  PSSK  formula  is  a 
function  of  target  hardness  in  PSI  and  the  yield  and  CEP  of 
the  weapon.  Th*  damage  function  is  linearized  and  through  a 
method  of  separable  convex  functions  can  be  solved  as  a 
transportation  problem.  By  using  a  numerical  example,  th* 
authors  show  that  round-off  errors  are  small  when  th*  number 
of  weapons  is  much  greater  than  th*  amount  of  targets,  and 
th*  round-off  error  is  large  when  the  weapons  are  roughly 
equivalent  or  less  than  th*  numbqr  of  targets  (Ref  9) . 

PERKINS.  This  model  develops,  optimal  defensive  strat- 
*  1 

egies  for  an  anti-missile  defense  and  can  be  used  for  optimal 

I 

offensive  allocation  strategies.  The  author  contends  that 
the  success  of  defending  units  is  jproportional  to  th*  ratio 
of  defending  units  to  attacking  units.  A  12-step  procedure 
is  developed  to  determine  optimal  offensive  and  defensive 
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strategies.  The  method  is  -fundamentally  based  on  marginal 
rates  of  effectiveness  which  are  represented  by  partial 
derivatives  of  the  damage  function.  Allocations  are  made 
until  diminishing  returns  are  realized  (Ref  22> . 

MATHESON .  This  is  a  game-theoretic  approach  to  missile 
allocation  in  terms  of  costs  or  resources.  There  are  four 
cases  considered: 

Case  1.  Defense  has  no  available  resources. 

Case  2.  Defense  has  less  than  enough  resources  to 
equalize  offense  system  costs. 

Case  3.  Defense  has  exactly  enough  resources  to 
equalize  offense  system  costs. 

Case  4.  Defense  has  more  than  enough  resources  to 
equalize  offense  system  costs. 

Employing  game-theory  min-max  principles  enables  the 
model  to  generate  explicit  expressions  for  both  offensive  and 
defensive  acquisition  strategies  in  terms  of  unit  costs.  Unit 
cost  coefficients  for  two  offensive  and  two  defensive  systems 
are  considered  in  the  allocation. 

Summary 

The  models  examined  range  from  very  detailed  to  aggre¬ 
gated,  and  from  scenario  specific  to  generalized.  Most 
models  contain  the  same  essential  components  and  differ  only 
in  the  assumptions  associated  with  these  components.  The 
most  "realistic"  models  express  parameters  such  as  proba¬ 
bility  of  kill  as  functions  of  target  and  weapon  input  data. 


It  is  the  opinion  o-f  the  authors  of  this  thesis  that  the 
most  useful  models  are  those  that  permit  multiple  choices  of 
objectives  such  as  A EM  and  BSMADE. 
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Ill  Model  Characteristics 


This  chapter  presents  -four  characteristics  o-f  BRIK.  The 
-first  section  lists  both  model  restrictions  and  scenario 
limitations  and  the  second  section  presents  all  of  the  input 
requirements.  The  third  section  explains  how  single  shot 
probability  of  survival  <SSPS)  is  computed  and  the  final  part 
discusses  the  damage  function. 

Limi tat i ons 

R.D.  Specht  quotes  E.  S.  Quade  as  saying,  "All  of  the 
assumptions  of  a  model  must  be  made  expl  ici t. . .that  his  (the 
modeler's)  errors  will  be  more  evident"  (Ref  26:219).  While 
assumptions  may  not  necessarily  reflect  errors,  they  do  rep¬ 
resent  the  modeler's  perception  of  the  "real  world."  Con¬ 
sequently,  it  is  important  for  any  user  to  be  aware  of  judge¬ 
ments  and  limitations  built  into  any  program.  Therefore, 
BRIK's  limitations  are  explicitly  brought  forth,  not  to  ex¬ 
pose  "error",  but  to  provide  a  further  understanding  of  the 
basis  for  the  model's  construction.  This  section  is  divided 
into  two  parts.  The  first  part  lists  the  restrictions  inher¬ 
ent  in  the  model,  and  the  second  part  describes  limitations 
that  must  be  imposed  on  any  scenarios.  While  these  limita¬ 
tions  are  a  direct  result  of  the  model's  restrictions,  they 
are  included  in  a  separate  subsection  to  better  express  the 
shortfalls  a  user  might  find  unacceptable. 

Model  Restrictions.  The  development  and  potential  appli¬ 
cation  of  BRIK  are  based  on  the  following  items: 
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1)  Aggregation.  SRIK  is  a  fully  aggregated  nuclear 
exchange  model.  All  of  the  weapons  and  targets  used  for  any 
allocation  must  be  represented  by  classes.  Ulithin  each 
class,  any  subsystem  must  have  identical  characteristics. 

For  example,  all  items  in  a  target  class  called  BRIDGE  must 
have  the  same  input  parameters  of  size,  vulnerability,  and 
value. 

2)  No  individual  targeting  strategy.  The  output  con¬ 
sists  of  the  number  of  weapons  from  each  weapon  class  that 
are  used  to  attack  a  particular  target  class.  If  more  than 
one  weapon  class  is  allocated  to  a  single  target  class,  SRIK 
does  not  provide  the  individual  breakdown  concerning  which 
weapons  go  to  which  targets. 

3)  Non- integer  solutions.  No  attempt  in  made  to  Keep 
the  weapon  numbers  in  whole  units.  For  example,  the  model 
could  suggest  assigning  24. 6  weapons  of  type  i  to  target 
class  J . 

4)  No  footprinting  restr i ct i ons.  In  reality,  multiple 
independently  targeted  re-entry  vehicles  <MIRV>  from  a  single 
missile  are  restricted  in  the  degree  of  dispersion  of  their 
targets.  BRIK  does  not  account  for  this. 

3)  Prompt  damage  only.  The  only  nuclear  effects  con¬ 
sidered  when  computing  SSPS  are  damage  due  to  overpressure , 
dynamic  pressure,  blast  duration,  and  cratering.  Prompt  and 
residual  radiation  from  thermal  effects,  neutrons,  and  gamma 
rays  is  not  considered. 
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6)  No  collateral  damage.  Each  target  in  each  class  is 


assumed  to  be  separate.  For  example,  if  two  targets  are 
close  to  each  other,  destruction  of  one  will  have  no  effect 
on  the  other. 

7)  No  nuclear  fratricide.  There  exists  the  possibility 
that  nuclear  expl osi ons  wi 1 1  adversely  affect  the  detonation 
of  subsequent  warheads.  This  is  called  fratricide.  BRIK 
does  not  take  fratricide  into  account. 

Scenario  Limitations.  Any  research  design  concerning  the 
study  of  nuclear  forces  must  consider  scenarios.  A  force 
structure  that  looks  good  under  one  plan  may  be  totally 
ineffective  in  another.  For  example,  a  force  composed  en¬ 
tirely  of  very  accurate  land-based  missiles  may  look  extreme¬ 
ly  impressive  if  the  scenario  calls  for  those  missiles  to 
attack  known  targets.  However,  if  the  current  policy  dic¬ 
tates  they  will  only  be  used  for  retaliation,  the  force  would 
be  very  ineffective  if  they  could  be  destroyed  in  a  first 
strike.  Since  scenario  is  so  important,  the  limitations 
imposed  on  any  attack  scheme  are  explicitly  stated. 

1>  The  target  and  weapon  base  must  be  completely  known. 
The  number,  diameter,  value,  and  some  expression  of  vulnera- 

l 

bility  must  be  known  for  each  target.  For  each  weapon,  the 
user  must  know  the  number  of  delivery  vehicles,  the  number  of 
warheads  Iper  weapon,  circular  error  probable  <CEP>,  the  yield 
of  each  warhead,  a  reliability  figure  expressing  the 


23 


V.  :  i  .  /  '  .  ;  '  •  ,  V  . 


■  •.  / 


probability  that  the  weapon  will  successfully  arrive  at  its 
destination  and  detonate,  and  daily  and  generated  alert  rates. 

2)  Any  exchange  can  only  be  represented  as  a  sequence  of 
attacks.  Before  Blue  can  retaliate  to  a  Red  strike,  it  must 
be  assumed  that  the  Red  strike  has  been  completed. 

3)  Time  is  not  explicitly  modeled.  The  model  does  not 
account  for  the  fact  that  some  weapons  will  detonate  earlier 
than  others.  However,  it  is  possible  to  restrict  seme  wea¬ 
pons  from  attacking  certain  targets.  This  permits  considera¬ 
tion  of  time  urgent  targets. 

4)  Point  defenses  are  not  considered.  The  advantages 
gained  from  a  defensive  system  placed  only  in  one  geographic 
location  cannot  be  studied  with  BRIK.  However,  the  reliabil¬ 
ity  percentage  of  an  individual  weapon  class  can  be  used  to 
study  the  effect  a  particular  area  defensive  system  has  on 
the  overall  problem.  For  example,  if  it  is  determined  that  a 
particular  defensive  system  degrades  a  weapon  class  perfor¬ 
mance  by  20X,  this  percentage  can  be  included  in  the  rel la¬ 
bility  factor,  REL,  for  that  particular  weapon.  However,  this 
factor  will  now  apply  whether  the  weapon's  attack  is  in  a 
heavily  defended  area  or  on  an  undefended  target. 

5)  Imprecisely  located  targets  cannot  be  studied.  One 
of  the  assumptions  of  BRIK  is  the  complete  description  of  the 
target  and  weapon  bases.  UJhile  coordinates  of  each  target 
are  not  an  input  requirement,  it  is  assumed  that  the  location 
of  each  target  is  known. 


BRIK  does  not 


6)  Command  and  Control  is  not  considered, 
account  for  the  effects  from  loss  of  any  positive  control 
elements. 

This  concludes  the  limitations  the  model  places  on  the 
user.  Next  the  input  requirements  are  presented. 

Input  Requirements 

To  perform  an  allocation,  BRIK  requires  complete  descrip¬ 
tion  of  the  weapon  and  target  base.  Table  I  lists  the  input 
requirements  for  the  weapon  complex.  Table  II-A  gives  the 
inpi’t  requirements  for  each  target  class  and  Table  II-B  lists 
the  additional  input  requirements  if  BRIK  is  asked  to  com¬ 
pute  a  value  associated  with  force  targets.  Using  this  data, 
BRIK  internally  computes  an  SSPS  for  each  weapon/target  com¬ 
bination.  The  methods  used  for  the  SSPS  computations  and  the 
appropriate  formulas  are  described  in  the  next  section. 

Single  Shot  Probability  of  Survival 

Most  nuclear  exchange  models,  including  AEM,  consider 
only  prompt  damage  mechanisms  to  compute  SSPS.  Even  though 
this  underestimates  the  total  destruction  resulting  from  a 
nuclear  detonation,  the  vogue  of  using  only  prompt  damage 
will  be  continued  with  this  effort. 

This  model  has  the  capability  to  consider  four  prompt 
damage  effects  when  computing  SSPS.  These  effects  are  over¬ 
pressure,  dynamic  pressure,  blast  wave  duration,  and  crater¬ 
ing.  Overpressure  is  the  transient  pressure,  usually  expres¬ 
sed  in  PSI ,  exceeding  the  ambient  pressure,  manifested  in  the 
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TABLE  I 

Input  Requirements  for  the  Weapon  Complex 


wpname< i > 

Name  of  each  weapon  class  i 

numwpn< i ) 

Number  of  delivery  vehicles  in  each 
class 

numwrh< i ) 

Number  of  warheads  per  delivery 
veh i cl e 

rel < i ) 

Probability  weapon  class  i  will 
reach  its  target  and  detonate 

wpncep<i) 

Circular  error  probable  <CEP)  of 
weapon  i  in  nautical  miles 

wpnyl d< i ) 

Yield  in  megatons  of  each  weapon  i 

wpndda( i ) 

Daily  alert  rate  of  weapon  class  i 

wpnga( i ) 

Generated  alert  rate  of  weapon 
class  i 

TABLE  1 1 -A 

Input 

Requirements  for  the  Target  Complex 

tgtnam< j ) 

Name  of  each  target  class  j 

numtgt< j ) 

Number  of  targets  in  each  class  j 

tgtpsi <J> 

Vul nerabi 1 i ty  of  each  target  class  j 
expressed  as  a  VNTK  or  PSI  number 

tgtdi a< j  > 

Diameter  in  nautical  miles  that 
encompasses  95X  of  a  target  in  class  j 

tgtcat< j ) 

Target  category.  F,  V,  or  M  for  force, 
val ue  or  mi  1 i tary . 

ntgtpr< j  > 

Preemptive  priority  of  target  class  j 

tgtval < j ) 

Number  reflecting  the  value  of  each 

' 

target  in  class  j 

TABLE  II-B 

Input  Requirements  to  Compute  Force  Target  Value 


tgtrel < j ) 

Reliability  of  weapons  at  target 
class  j 

tgtcep< j) 

CEP  of  weapons  at  target  class  j 

tgtyl d< j  > 

Yield  per  warhead  in  megatons  of  each 
weapon  at  target  class  J 

ntgtwh< j ) 

Number  of  warheads  located  at  target 
class  j 
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shock  wave  -from  an  explosion.  Dynamic  pressure  is  the  air 
pressure  which  results  -from  the  mass  air  flow  (or  wind) 
behind  the  shock  front  of  a  blast  wave.  Blast  wave  duration 
is  the  time  either  the  overpressure  or  dynamic  pressure  is 
present,  and  cratering  is  the  pit,  depression,  or  cavity 
formed  in  the  surface  of  the  earth  by  a  surface  or  under¬ 
ground  explosion  (Ref  14:632-637).  A  description  of  the 
methods  available  to  consider  these  damage  effects  and  the 
formulas  that  appear  in  BRIK  are  presented  next. 

Vulnerabi 1 i ty.  The  vulnerability  of  a  target  class  can 
be  expressed  either  as  a  VNTK  figure  or  a  PS1  number.  De¬ 
pending  on  the  input  from  the  user,  BRIK  computes  the  SSPS  of 
a  weapon/target  combination  from  one  of  two  methods.  The 
first  method  is  called  the  Physical  Uulnerabi 1 i ty  System.  In 
this  system,  a  target's  susceptibility  to  blast  damage  is 
indicated  by  a  combination  of  numbers  and  letters.  The 
vulnerability  number,  WTK,  consists  of  a  two-digit  number 
reflecting  the  target  hardness  relative  to  a  specified  damage 
level,  followed  by  a  letter  indicating  predominant  sensitiv¬ 
ity  to  overpressure  (P>  or  dynamic  pressure  (Q),  then  finally 
a  K  factor.  The  K  factor  allows  for  hardness  adjustments  to 
be  made  to  account  for  the  effects  of  variations  in  blast 
wave  duration  due  to  different  weapon  yields  (Ref  3:34).  The 
following  example  indicates  the  importance  of  considering 
blast  wave  duration.  If  the  overpressure  duration  is  short 
(1  to  3  milliseconds),  overpressures  between  390  PSI  and  470 
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PSI  are  required  to  achieve  a  SOX  mortality  rate  in  humans. 
However,  animal  tests  have  shown  that,  if  the  overpressure 
duration  is  long  <80  to  100  milliseconds),  40  to  50  PSI  will 
produce  that  same  50X  mortality  rate  (Ref  8:7-92).  Because 
buildings  are  also  susceptible  to  blast  duration,  the  LNTK 
system  has  been  included  in  this  moael.  However,  VNTK  num¬ 
bers  may  not  always  be  available;  therefore,  the  less 
accurate  method  of  designating  an  estimated  PSI  level  has  also 
been  included.  This  entails  estimating  the  peak  overpressure 
in  pounds  per  square  inch  (PSI)  that  will  provide  the  desired 
level  of  damage.  A  discussion  concerning  appropriate  PSI  | 

i 

numbers  can  be  found  in  Reference  27,  page  31. 

VNTK.  The  ALLOT  computer  program,  currently  a  part  of 
the  SAC  BLUE  SUADE  model ,  uses  the  VNTK  system  to  compute  i 

SSPS.  A  listing  of  BLUE  SUADE  was  acquired  and  DPDX,  the  j 

i 

subroutine  that  computes  probability  of  kill  (PK),  was  ex-j 

I 

tracted.  After  a  few  minor  modifications,  that  subroutine' 

I 

was  placed  in  BRIK  as  Function  VTK.  It  is  known  that  LtCol 

j 

Donald  J.  Berg  wrote  Allot  (Ref  23:100);  however,  little  is 
known  about  the  regression  techniques  used  to  develop  the 
formulas  in  DPDX.  It  is  certain,  though,  that  the  subroutine 
is  being  used  properly  (see  Verification  in  Chap  VI).  First, 
this  section  will  discuss  the  input  requirements  for  DPDX  and 
how  it  is  being  used  in  BRIK.  Then,  the  formulas  and  the 
data  table  in  the  subroutine  will  be  listed. 
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The  parameter  statement  associated  with  DFDX  in  ALLOT 

/ 

contains  ten  elements.  The  first  three-- LN,  T,  and  K1  — 
are  ';he  three  parts  to  the  LfrfTK  number  input  by  the  user. 

The  next  three  elements  concern  weapon  and  target  character¬ 
istics.  They  are  95%  of  the  radius  of  a  circle  that  encom¬ 
passes  the  target,  weapon  yield,  and  weapon  CEP.  The  seventh 
element,  H,  is  always  either  zero  or  one.  It  denotes  whether 
the  weapon  detonation  was  an  air  burst  or  a  ground  burst. 

The  last  three  elements  concern  numbers  computed  by  the 
subroutine.  The  first  one  is  probability  of  kill  and  the 
last  two  are  used  in  a  sensitivity  printout  in  the  output 
section  of  the  program.  lit  revising  the  subroutine,  the  two 
elements  associated  with  the  sensitivity  output  were  removed. 
Also,  the  equations  that  computed  those  figures  were  taken 
out.  Next,  a  loop  was  included  inside  the  subroutine  to 
compute  PK  assuming  the  attacking  weapon  detonated  as  both  an 
airburst  and  a  ground  burst.  This  negated  the  requirement  to 
designate  whether  or  not  a  weapon  detonated  in  the  air  or  on 
the  ground.  To  compute  SSPS,  BRIK  uses  the  highest  PK  com¬ 
puted  from  these  two  cases.  Following  are  the  formulas  and 
the  data  tables  used  in  Function  VTK. 

PK  »  EXPC-EXPCb<l ,T)+U*<b<2,T>+U*<b<3,T>+U*b<4,T>>>>> 
where 

T  =»  1  if  target  is  susceptible  to  over-pressure 

T  **  2  if  target  is  susceptible  to  dynamic  pressure 

U  -  S/<XD*C> 
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XD  a  l/Y^1^ 

C  -  /cep2  +  .231*  R2 

R  «  radius  of  a  circle  in  -feet  that  encompasses  95X  of 
the  target 

S  *  EXPCa< 1 ,H,T)+V#<a<2,H,T>+V*< a<3,H,T)+V*<a<4,H,T> 
+V<a<5,H,T>+V*<a<<S,H,T)+V*a<?,H,T>>>>>>> 

H  *  0  or  1  representing  airburst  or  ground  burst 

V  »  VN+d<l ,K,T>+XD*<d<2,K,T)+XD*<d<3,K,T>+XD*d<4,K,T>>> 

K  »  K  Factor 

Y  »  Yield  in  Kilotons 

CEP  *  Circular  Error  Probable. 

It  should  be  noted  that  there  are  limits  in  the  subroutine. 

If  U  is  greater  than  10,  PK  is  automatically  set  to  1.  Also, 
if  U  if  less  than  0.1,  PK  is  assigned  0. 

The  final  step  in  acquiring  SSPS  is  to  include  weapon 
reliability.  The  following  formula  is  used  for  each  weapon 
target  combination: 


SSPS..  =  1  -  REL.*PK 

where  SSPSjj  represents  the  single  shot  probability  of  survi¬ 
val  of  a  target  in  class  j  assigned  a  weapon  from  class  i, 
and  REL^  is  the  reliability  of  weapon  class  i. 

There  are  three  arrays  in  the  above  f ormul as ' that  obtain 
their  data  from  tables.  Array  a(7,2,2)  is  a  three-dimen¬ 
sional  array  that  can  be  thought  of  as  having  two  sets  of  two 
rows  of  numbers,  each  row  having  seven  elements.  Those  23 
numbers  are  in  Table  III.  Array  b(4,2)  has  two  rows  of  four 
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numbers  each.  Its  elements  ere  in  Table  IV.  Finally,  array 
d<4,10,2>  has  two  sets  of  ten  rows,  each  row  having  four 
numbers.  These  numbers  are  listed  in  Table  V. 

This  concludes  the  description  of  Function  VTK.  BRIK 
automatically  distinguishes  between  a  WTK  number  and  PSI . 

If  a  PSI  number  is  input,  Function  PK  is  used  to  compute 
SSPS.  The  formulas  in  that  function  are  presented  next. 

PSI .  To  compute  SSPS  given  a  PSI  number,  the  model 
starts  by  computing  a  lethal  radius  <LR).  If  PSI  is  greater 
than  10,  the  lethal  radius  resulting  from  a  weapon  of  class  i 
attacking  a  target  of  class  J  is  expressed  as  follows  <Ref 
27*49)  i 

LRi;j  *  2.8Yi(1/3)(PSI;j-7.37)'*352.  (la) 

If  PSI  is  less  than  or  equal  to  10  the  lethal  radius  is  (Ref 
10*214) 

LRij  a  (6.8l*Y^^3^  )/psij  *62  <ib> 

where  Y^  equals  yield  in  megatons.  For  cratering,  it  is 
assumed  that  the  target  is  destroyed  if  it  is  inside  the  lip 
crest  of  a  crater  formed  in  soft  rock.  The  formula  to  com¬ 
pute  this  LR  is  derived  in  the  following  manner.  The  apparent 
crater  radius  for  a  one-kiloton  surface  burst  in  dry  soft 
rock  is  41  feet.  Scaling  the  yield  to  the  .3  power  and 
computing  the  lip  crest  radius, 

Rij  a  l‘25*6l*Y1*3 

where  Y^  equals  yield  in  kilotons  and  R^j  equals  the  crater 
lip  crest  radius  in  feet  (Ref  14*253).  In  BRIK,  the  formula 
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TABLE  III 


Data  for  Array  a<7,2,2> 


a(x,x,l) 

8. 214, -.1118,  5.265a-4,  2. 162#-3,-*6.638e-7,  7. 132e-9, -3.064e-l  1 
8.783,-. 1353,,  2.355*-3,-2.086e-4,  9.901*-6,-l .872e-7,  1.227e-09 
a<x,x,2) 

8.313,-. 1033,-7. 908»-4, -9. 039e-5,  1 .458e-5,-5.220e-7,  3.726e-09 
8.789,-. 1120 , -6.658a-5,-5.803e-4,  5.853e-5,-l .905»-6,  2.056e-08 


TABLE  IV 


Data  for  Array  b(4,2) 

1.66904,  -2.17442,  .260926,  -.0752178 

1.65946,  -2.15466,  .295833,  -.0484718 


TABLE  V 


Data  for  Array  d(4,10,2> 


d(x ,x  ,  1 ) 

0.00000, 

0.00000, 

0.00000, 

0.00000 

•• 

0.377874 

,  1.56916, 

0.0013762, 

0.0063101 

- 

1 .22391 , 

3.32978, 

-0.0020688, 

-0.0469539 

— 

1  .937, 

5.35163, 

-0.0477323, 

-0.134546 

- 

2.80456, 

7.74241 , 

-0.223298, 

-0.343665 

— 

3.81168, 

10.7222, 

-0.84178, 

-0.513504 

- 

5.05081 , 

14.6336, 

-2.37473, 

-0.476883 

— 

6.65796, 

20.1955, 

-5.98136, 

0.313062 

— 

8.92061 , 

28.9801 , 

-14.2174, 

3.11471 

- 

12.7263, 

45.9934, 

-35.3644, 

12.2164 

d<x ,x ,2> 

0.00000, 

0.00000, 

0.00000, 

0.00000 

— 

0.288917 

,  0.79886, 

-0.0399065, 

0.0011525 

— 

0.611921 

,  1.72876, 

-0.18691 , 

0.0095419 

- 

0.978187 

,  2.83679, 

-0.511558, 

0.0514480 

- 

1 .40112, 

4.19642, 

-1.13330, 

0.18118 

— 

1 .90063, 

5.91 5*8, 

-2.23005, 

0.484749 

— 

2.50907, 

8.17913, 

-4.11606, 

1.12029 

— 

3.28374, 

11.3172, 

-7.35613, 

2.37622 

— 

4.34296, 

16.0486, 

-13.2112, 

4.90666 

“ 

6.01000, 

24.4300, 

-25.3672, 

10.643 
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used  expresses  radius  in  nautical  miles  and  yield  in  mega¬ 
tons.  Therefore, 

*  is§f^<vi000)‘3 

or, 

LHjj  «  .125*Y*3  <2> 

Cratering  is  independent  o-f  PSI .  Therefore,  this  formula 
places  a  ceiling  on  the  protection  gained  from  super-harden¬ 
ing.  Figure  1  plots  the  relationship  between  weapon  yield 
and  the  PSI  ceiling.  For  PSI  values  below  this  ceiling,  BRIK 
uses  formula  set  1  to  compute  LR.  For  PSI  values  above  this 


ceiling,  the  lethal  radius  is  computed  by  formula  2. 

Vsing  one  of  the  above  lethal  ridi i  formula,  the  circular 

'  ! 

normal  distribution  is  used  to  compete  the  single  shot  proba¬ 
bility  of  Kill  <SSPK)  for  a  point  target  in  class  j  allocate. 


a  weapon  from  class  i.  The  formula  Is  as  follows! 

SSPK1;j  a  l-.5(LRij/<[EPi)2 


where  CEP^  equals  the  circular  error  probable  of  weapon  i 
LR^  equals  the  larger  of  the  two  Itthal  radii  computed 


and 


above . 


This  formulation  assumes  what  is  commonly  called  a 
"cookie-cutter"  damage  function!  if  the  warhead  lands  at  or 
within  the  lethal  radius,  the  target  is  completely  destroyed; 
if  the  warhead  lands  outside  the  lethal  radius,  the  target  is 
completely  undamaged.  However,  in  reality,  this  does  not 
happen.  Because  of  variability  in  target  hardness,  war¬ 
head  effects,  and  the  hardness  of  different  target 
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_ _  YIELD  IN  MEGATONS _ _ 

Figurt  1.  Boundary  Between  LR  Formulas  of  Cratering 
and  PSI 

components,  there  is  some  region  around  a  target  where  a 
weapon  can  be  detonated  and  the  target  is  neither  completely 
destroyed  nor  completely  unharmed.  To  represent  this  region, 
a  model  could  use  a  log-normal  damage  function.  A  technique, 
sometimes  referred  to  as  a  quadrature  scheme,  is  used  to 
compute  SSPK.  Given  CEP  and  LR,  points  are  established  at 
Known  distances  from  a  target.  The  probability  that  the  weapon 
hits  each  point  is  multiplied  by  the  probability  that  the  weapon 
destroys  the  target,  assuming  the1  weapon  hit  each  point.  The 
average  of  these  products  provides  the  probability  of  Kill. 

This  technique  is  not  used  here  because  the  added  complexity 
of  computing  SSPK  usirg  a  log-normal  function  results  in  an 
insignificant  change  in  probability  of  survival — about  2.2 
percent  (Ref  4x16). 

If  the  area  of  the  target  class  J  is  large  in  comparison 
to  a  weapon's  lethal  diameter,  a  different  interpretation  of 
damage  is  needed.  For  a  point  target  in  class  j,  SSPK^j  is 

39 


.rat 


the  probability  a  weapon  i  will  detonate  with  the  target 
located  within  the  weapon's  lethal  radius.  However,  for  an 
area  target  J,  SSPK^j  can  be  interpreted  as  the  percentage  of 
the  target  destroyed  if  allocated  one  warhead  of  type  i.  To 
account  for  this  change  of  interpretation,  an  adjustment  is 
made  to  the  SSPK^^  computed  for  a  point  target.  The  fol¬ 
lowing  formula  was  found  in  Reference  27,  page  41,  credited 
to  Reference  20: 

poo  »  SSpKij[§a-(i-§)Bmi-§)(§)(1/B)] 


where 


X  »  <Lyo> 

L  ■  Lethal  Diameter 
D  ■  Target  Diameter 


*  c 

'1C2  y  3 

0<y<l 

1  ♦  040jyy°« 

1  <y<  °* 

C7C8yy°9 

0<y<2 

°10'0ll012yy°13 

2<y<~ 

y  -  LR/CEP. 


The  values  of  cl  to  cl3  are  as  follows: 


cl 

■ 

7.43993 

c3 

■ 

.0930393 

c9  • 

-.29298 

c2 

- 

.207333 

c4 

- 

3.83323 

clO  - 

2.73109 

c3 

- 

2.4242 

c7 

- 

1 .489 

ell  - 

1.12109 

c4 

• 

3.89243 

c8 

. 

1.09218 

c!2  - 

.79845 

.439752 


X 


The  final  step  to  acquire  SSPS  includes  weapon  system 


reliability.  Therefore, 

SSPS..  ■ 


1-REL. *SSPK. . 

1  i  J 


where  SSPS..  is  the  single  shot  probability  of  survival  of  a 
^  J 

target  in  class  j  assigned  a  weapon  from  class  i,  and  REL^  is 
the  reliability  of  weapon  class  i. 

This  concludes  the  presentation  of  BRIK's  computational 
techniques  for  probability  of  survival.  How  the  single  shot 
probability  of  survival  for  a  tsrget  in  class  j  assigned  a 
weapon  in  class  i  is  used  in  the  damage  function  is  discussed 
next . 

Damage  Function 

The  damage  function  in  BRIK  works  with  expected  value 

concepts.  Specifically,  the  measure  of  merit  involves  the 

expected  percentage  of  targets  in  a  class  that  survives  an 

allocation.  For  example,  let  S. .  denote  the  survival  probabi- 

*  J 

1 i ty  of  a  class  J  target  if  allocated  a  single  type  i  weapon. 
If  more  than  one  type  i  weapon  was  assigned  to  target  j,  the 
survival  probability  would  become 

siiXiJ 

where  X..  represents  the  number  of  class  i  weapons  assigned  to 

**■  J 

target  j.  If  there  were  several  different  types  of  weapons 
that  could  be  assigned  to  this  single  target,  the  survival 
probabi 1 i ty  could  be  expressed  as 

Tfs.  ,Xij 
i  d 

where  S.  •  represents  the  survival  probabi 1 i ty  of  target  j  if 

X  J 


allocated  a  type  i  weapon.  For  example,  if  target  j  has  an 
SSPS  of  .4  if  assigned  a  type  1  weapon  and  .3  if  assigned  a 
type  2  weapon,  the  probability  that  target  j  would  survive  if 

assigned  one  weapon  of  each  type  is 

2 

T  S.Aj  *  (.4)1(.3)1  s  .12 
i=l 

Next,  assume  that  there  are  two  more  targets  in  class  j  for  a 

total  of  three  targets.  Also,  assume  that  the  second  target 

was  allocated  one  type  1  weapon  and  that  the  third  target  was 

allocated  one  type  2  weapon.  The  probabi 1 i ty  of  survival  for 

the  second  target  is 

T  S±1Xij  =  (.4)1(.3)°  =  A 
i=l 

and  the  probability  of  survival  of  the  third  target  is 

TfS^ij  -  (.4)°(.3)1  ■  .3. 

i*l 

The  average  probability  of  survival  of  the  entire  class  is 

(l/3)(.12-.4+.3)  -  .273. 

Expressed  mathematically,  the  survival  function  for  class  J 

**  a.  m 

(1/a.)  t  (  ITs^ilc) 

3  k*l  i-1  13 

where  a^  represents  the  number  of  targets  in  class  J  and  m 
equals  the  number  of  weapon  classes.  This  expression  for  the 
average  probability  of  survival  is  the  sum  of  products  and  tis 

difficult  to  deal  with  in  terms  of  solvability.  Therefore,^ 

1 

the  following  approximation  is  used  in  BRIK.  The  average  j 
probability  of  survival  of  the  entire  class  j  can  be  aproxi-\ 


mated  by 


when  Xj  ^  represents  the  average  number  of  weapons  allocated 
against  each  target  in  class  j.  In  the  current  example,  a 
total  of  two  weapons  from  each  weapon  class  was  assigned  to 
the  three  targets  in  class  j.  Using  the  approximation  to 
compute  the  average  probability  of  survival  yields 

\  S.  .Xij  =  (.4)(2/3)(. 3)<2/3)  s  .243 
i=l  1J 

This  compares  with  .273  obtained  from  using  formula  <C  1  >  .  If 
Xij  *op  *  are  *nte9ersi  there  is  no  error,  and  the  larg¬ 

est  error  found  occurs  when  X.*  is  equal  to  .5,  1.5,  2.5, 

J 

,etc.  Both  Reference  27  and  Reference  17  agree  that  the 
magnitude  of  the  error  in  this  approximation  is  "believed  to 
be  small;"  This  belief  combined  with  the  increase  in  solva¬ 
bility  convinced  the  thesis  authors  to  use  it  in  8RIK. 

In  conclusion,  the  damage  function  employed  in  BMK  is  an 
approximation  that  expresses  the  probability  of  survival  of  a 
target  class  when  allocated  different  weapons.  How  that 
function  is  used  in  the  linear  constraints  of  a  goal  program¬ 
ming  formulation  is  discussed  in  Chapter  IV. 
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IV  Constraint  Development 


BRIK  uses  a  linear  goal  programming  routine  to  allocate 
weapons  to  targets.  Goal  programming  requires  two  things,  an 
objective  -function  and  a  constraint  set.  A  discussion  o-f  ;,oal 
programming  and  a  presentation  of  the  objective  functions  in 
BRIK  will  be  found  in  Chapter  V,  while  this  chapter  derives 
all  of  the  constraints  used  by  BRTK  in  the  linear  programming 
matrix.  Regardless  of  the  problem,  the  constraint  set  is 
identical.  This  set  consists  of  one  constraint  for  each 
target  class,  one  constraint  for  each  weapon  class,  one 
constraint  for  the  extreme  noal ,  and  one  constraint  for  each 
hedging  option  designated  by  the  user.  As  currently  dimen¬ 
sioned,  BRIK  can  handle  a  maximum  of  61  constraints. 

Target  Constraint 

Let  Sjj  represent  the  single  shot  probability  of  survival 

of  a  target  in  class  j  allocated  one  weapon  of  type  i.  If  X, . 

**•  J 

represents  the  number  of  weapons  i  allocated  to  each  target 
in  class  j,  the  survival  of  target  j  can  be  represented  by 

SijXi^ 

If  different  types  of  weapons  i  are  allocated  to  target  class 
j ,  survival  of  each  target  becomes 

In  a  goal  programming  problem,  a  constraint  containing  the 
survival  expression  and  a  damage  expectancy  <DE)  goal  can  be 
written  as  follows: 

*  1-DEj  •  pvj 
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where  PV.  equals  < 1-DE  •> ,  the  survival  goal  of  a  target  in 
J  J 

class  j;  n j  equals  the  negative  deviation  from  the  survival 

goal;  and  p.  equals  the  positive  deviation  from  the  survival 
J 

goal.  If  this  constraint  was  used  in  a  non-linear  goal  pro¬ 
gramming  problem,  minimizing  p.  would  drive  the  survival 

J 

expression  to  the  goal  PV j . 

It  is  possible  to  linearize  this  constraint  to  permit  its 


use  in  a  standard  linear  programming  model.  By  moving  the 
goal  variables  n^  and  p^  to  the  right-hand  side, 

fi/ij  •  pvrnj+pr 

If  PV-  is  divided  out  of  each  term  in  the  right-hand  side, 
J  —  y  -n4+p. 


j  s  pY1+-kf>- 

,  if  the  natural  logarithm  of  each  side  is  taken,  the 


following  is  attained! 


ln(ITs.  .xij)  ■  In  PV.(1+ 


-n<+ P; 


Because  the  logarithm  of  products  is  the  sum  of  logarithms, 

|  Zln(SijXij)  -  lnPVj+ln(l+~ 

and 

Xln(  S .  jXi  j )  -ln(  If-  J  y-  J )  =  lnPVj . 
i  j 

The  survival  function  will  always  be  greater  than  zero  and  less 
than  or  equal  to  one.  This  is  also  true  of  a  goal  PV ^ . 

Since  the  logarithm  of  any  number  between  zero  and  one  is 
negative,  it  is  necessary  to  multiply  this  constraint  by 


negative  one  to  obtain  a  positive  right-hand  side: 

-Zln(S.  .Xij)fln(l+-^rJ-)  =  -InPV. . 

It  is  also  possible  to  bring  the  exponent  X*  *  outside  the 

*  J 
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parenthesis  and  obtain 

lnSi  j  >xi j+ln<  1+~p^’)  *  -^y 

In  a  goal  programming  problem,  n.  represents  the  negative 

J 

deviation  from  a  goal  and  p.  represents  the  positive  devia- 

J 

tion  from  a  gou'l .  If  a  goal  is  not  met  exactly,  it  is 
impossible  to  simultaneously  deviate  from  that  goal  in  both  a 
positive  and  negative  direction.  Therefore  only  one  of  the 
variables,  n^  or  pj  ,  can  be  non-zero  at  any  given  time. 

Thus,  either  nj  or  Pj  must  always  be  zero.  Since  1n(l>  »  0, 
it  is  possible  to  split  the  second  term  in  the  above  equation 
and  obtain 


This  model  is  an  aggregated  nuclear  exchange  model . 
Therefore,  let  S^j  represent  the  survival  probabi 1 i ty  of  all 
class  j  targets  allocated  class  i  weapons.  Since  X* *  equals 
the  average  number  of  weapons  i  allocated  to  each  target  in 
class  J,  it  is  possible  to  perform  a  variable  transformation 
of  Xjj  ,  such  that 

xi3  *  *ij*J 


where  Yjj  represents  the  number  of  class  i  weapons 
allocated  to  target  class  j  and  N.  equals  the  number  of 

U 

targets  in  class  j.  Using  this  transf ormation  greatly  reduces 
round-off  error  resulting  from  non-integer  programming;  how¬ 
ever,  it  introduces  the  approximation  errors  mentioned  in 


4  6 


Chapter  III.  This  transformation  yields  the  following: 

-I(lnSi;j)(Yi;3/N.)f-ln(l+p^-)-(-ln(l+^i))  =  -lnPV^ 

1  J  J 


-I(lnSij)Yia+Njln(l+^)-(-Njln(l+^f))  =  -NjlnPV. 
l  J  J 


Since  0  <=*  n^  <*  , 


H.ln(l+^sp) 

3  FVj 


is  always  non-positive.  And,  since  0  <=  p^  <=  1  -  PVj  , 

p.  J  J 

'  J 

is  always  non-negative.  Therefore,  these  two  expressions  can 
represent  the  dev  i  at  i  onaK  war  i  abl  es  in  the  final  form  of  the 
target  constraint.  The  final  form  of  the  target  constraint 
used  in  this  effort  is 


-£(lnsij)llfdj--d/  =  -N.lnPVj 


where  dj  equals 


and  d*  equals 


Hjln(l+p7i) 

J 


-Njlnd+p^f-) . 

As  mentioned  above,  minimizing  pj  would  drive  the 

survival  expression  to  the  goal  PUs  .  Since  p*  appears  in  the 

J  j 

term  that  represents  the  negative  deviational  variable,  ck~, 

J 

minimizing  d*”  would  drive  this  constraint  to  the  right-hand 
J 

si-de  goal.  Likewise,  minimizing  dj+  would  effectively  place 
an  upper  bound  on  the  constraint,  preventing  any  over¬ 
achievement  of  the  designated  goal.  How  these  deviational 
variables  are  used  in  the  different  objective  functions  in 
BRIK  will  be  discussed  in  the  next  chapter. 
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This  is  not  the  first  time  a  non-linear  damage  function 
has  been  linearized  with  logarithms.  Both  Hodson  and 
Wambsganss  did  so.  Hodson 's  formulation  had  only  a  single 
deviational  variable,  sometimes  referred  to  as  minimizing  the 
maximum  deviation  (Ref  16s 178) ,  and  his  goals  were  ratios  of 
DE.  For  example,  if  there  were  three  target  classes  and  it 
was  desired  to  have  an  allocation  that  would  give  DE  in  the 
ratio  of  3:5:4,  his  PVj's  were  set  as  follows: 

PV1  »  In  3  PV2  =  In  5  PU3  =  In  4. 

This  worked  satisf actor i 1 y  for  him  (Ref  17:14).  However, 
just  as  the  constraints  are  being  used  here,  Mambsganss  used 
them  to  express  DE  goals  in  terms  of  a  percentage  of  survival 
for  each  class  j  and  each  constraint  had  a  different  devia¬ 
tional  variable  (Ref  27:71-74) .  This  sometimes  produces  some 
undesirable  side  effects.  The  cause  of  those  side  effects 
and  possible  solutions  will  be  discussed  in  Chapter  VII. 
Weapon  Constraints 

In  the  target  constraints,  the  decision  variable  Y- • 

J 

represents  the  number  of  type  i  weapons  assigned  to  all  tar¬ 
gets  in  class  j.  The  weapon  constraints  sum  all  of  the 
weapons  In  class  i  assigned  to  each  target  class  j.  This  sum 
is  set  equal  to  the  number  of  warheads  available  in  each 
weapon  class  i.  This  is  depicted  by  the  following: 


fid  ■  Ai 

J 


where  A^  equals  the  number  of  warheads  available  in  class  i 
To  maintain  the  goal  programming  formulation,  deviational 
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variables,  d^~  and  d^  ,  are  added.  This  gives  the  -final 
-form  of  the  weapon  constraints: 

^ij^i  “di*  =  Ai* 

The  presence  o-f  these  deviations]  variables  provides  a  tremen¬ 
dous  amount  o-f  -flexibility.  Both  variables  are  used  at 
various  times  in  BRIK's  objective  functions.  Minimizing  d^+ 
prevents  the  use  o i  more  weapons  than  are  available,  while 
minimizing  d^”  forces  the  use  of  as  many  weapons  as  possible. 
As  seen  in  the  next  section  on  hedging  constraints,  complete¬ 
ly  different  interpretations  of  the  same  constraint  result 
from  switching  the  deviations!  variables  that  appear  in  the 
objective  function. 


Hedging  options  permit  a  user  to  designate  a  variety  of 
allocation  rules.  In  its  current  form,  the  model  has  seven 
options  to  choose  from.  As  currently  dimensioned,  up  to  20 
hedging  constraints  can  be  added.  In  this  section,  the  seven 
hedging  options  and  the  form  of  their  constraints  are  listed. 

The  following  is  a  list  of  the  seven  types  of  hedging 
options  available: 

1)  Enforce  a  minimum  level  of  damage  on  a  particular 

target  class.  \ 

\ 

2)  Enforce  a  minimum  Uevel  of  damage  on  a  particular 

V 

target  class  using  a  specific  set  o-f  weapons. 


3)  Enforce  an  upper  level  of  damage  on  a  particular 


target  class  resulting  from  a  specific  set  of  weapons 


4)  Restrict  the  number  o-f  a  class  of  weapons  which  can 
be  allocated  to  a  target  class. 

5)  Build  your  own  constraint. 

6)  Enforce  a  minimum  level  o-f  damage  on  a  particular  set 
o-f  target  classes. 

7)  Restrict  the  number  o-f  weapons  which  can  be  allocated 
to  each  target  in  a  particular  target  class. 

There  is  a  constraint  associated  with  each  o-f  the  above 
list  o-f  options.  For  option  one,  the  constraint  is  identical 
to  a  targeting  constraint.  Usually  the  DE  entered  here  would 
be  lower  than  the  normal  achievement  goal . 

Min  {d“l 
st  (subject  to) 

-plnSi^)Yi;j+d"-d+  *  -N^lnPVj 

where  the  variables  are  as  previously  defined.  For  hedging 
type  two,  the  constraint  differs  from  type  1  in  that  not  all 
weapon  types  must  be  in  the  constraint.  The  user  can  desig¬ 
nate  a  subset  k  of  weapons  to  be  included  in  this  hedge. 

Min  jd"j 
st 

*  -NjlnPY 

Type  3  hedge  can  be  used  to  restrict  the  damage  resulting 
from  a  designated  set  of  weapons.  There  is  no  difference 
between  this  constraint  and  that  of  hedge  2.  The  difference 
comes  from  minimizing  the  positive  deviational  variable  in¬ 
stead  of  the  negative  deviational  variable. 
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-  Z  (InS. . )Y. .fd”-d+  *  -N.lnPV. . 
ifk  J  j 

With  any  particular  target  class,  a  DE  goal  -for  the 

targeting  constraint  would  normally  be  higher  than  the  DE 
goal  for  this  type  3  hedge.  The  purpose  of  this  constraint 
is  to  further  restrict  the  damage  from  some  subset  of  wea¬ 
pons.  It  would  not  be  needed  if  the  DE  goal  cn  this  con¬ 
straint  was  set  higher  than  the  overall  goal  for  the  class. 

Type  4  hedge  sets  the  decision  variable  Yjj  equal  to  the 
maximum  number  of  type  i  weapons  that  can  be  assigned  to 
target  class  j.  Deviations!  variables  are  included  to  pre¬ 
serve  the  G.P.  formulation. 

Min{d+/ 

st 


Yij+d*’~d+  S  RHS* 

Type  5  hedge  consists  of  customizing  a  constraint.  For 
example,  a  user  may  have  an  unlimited  number  of  two  different 
types  of  warheads  but  a  limited  number  of  del ivery  systems. 

A  constraint  could  be  built  that  restricts  the  total  number 
of  the  two  warheads  delivered.  For  example,  assume  there  are 
already  two  weapon  constraints  in  the  problem  such  that 

fu*'-*1  ’  Ai 

and 

faj*'- d*  s  A2 

where  A^  and  A£  represent  large  warhead  availabilities,  but 


the  delivery  vehicles  -for  type  one  and  type  two  warheads 


restrict  the  total  number  o-f  deliverable  warheads  to  A^  . 

The  constraint  would  be  as  -follows: 
min  (d+| 
st 

f  u+f2.rd‘-d'  *  a3- 

Minimizing  d+  would  prevent  A^  -from  being  exceeded.  BRIK 
gives  the  user  the  ability  to  build  any  constraint.  All  that 
is  needed  is  i,  j,  the  coe-f-f icient  o-f  each  required  term, 
and  the  appropriate  right-hand  side. 

Type  6  hedge  permits  designating  a  minimum  level  o-f 
damage  across  a  set  k  o-f  targets. 

Min  (d~l 
st 

-  r  (KlnS,  JY,  ,)+d“-d+  -  -(ZN.)lnPV,. 
jfki  1J  13  j  3  3 

Type  7  hedge  restricts  the  number  of  weapons  which  can  be 

allocated -to  each  target  in  a  particular  target  class.  The 

constraint  sums  the  variabl  >s  th  .t  repres  .it  the  weapons  that 

can  be  assigned  to  some  class  J  end  se  :s  that  sum  equal  to 

some  multiple  M  o-f  the  number  of  targets  in  class  j. 

Min  (d+J 

st 

^Yij+d"~d+  =  my 

This  concludes  the  presentation  of  all  seven  hedging 
options  included  in  the  model.  A  great  variety  of  allocation 
rules  are  thus  permitted  allowing  tremendous  flexibility. 
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Extrema  Goal 


A  lowest  priority  goal  is  built  into  the  constraint  set 
to  give  the  model  the  ability  to  drive  the  allocation  to  a 
unique  solution.  I f  there  are  weapons  remaining  -following  an 
•  1  location,  there  are  an  in-finite  number  o-f  allocations 
available  that  meet  the  attained  DE  goals.  To  insure 
that  the  allocation  is  unique,  a  lowest  priority  goal,  *n 
extreme  goal,  is  included.  The  constraint  can  take  one 
o-f  two  different  forms, 

minfd/j 

**  £(£coefiYi j )+de"-de+  *  0  (1) 

where  coef^  is  internally  computed  based  on  a  choice  of  the 
user,  or 

min{de~] 

**  U  L  YiJ+d  "-d  +  =  TA.  (2) 

i  jek  e  e  i  1 

where  k  is  the  set  of  targets  chosen  to  receive  the  remain¬ 
ing  weapons. 

The  user  has  five  choices  of  extreme  goals.  They  are 

1)  Minimize  warheads  used.  _  . _ _ 

2)  Minimize  the  amount  of  megatonnage  used. 

3)  Minimize  the  countermi 1 i tary  potential  of  the  force. 

4>  Minimize  the  equivalent  megatonnage  of  the  force. 

9)  Use  as  much  of  the  remaining  arsenal  as  possible. 

Constraint  one  is  used  for  the  first  four  choices  and 
constraint  two  is  used  for  choice  five.  If  minimizing 
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warheads  is  chosen,  coef^  equals  1  in  all  cases.  If  minimi¬ 
zing  megatonnage  is  chosen,  coef^  is  equal  to  the  megatonnage 
of  weapon  class  i.  I-f  the  user  wishes  to  minimize  the  counter¬ 
military  potential  o-f  the  allocated  -force,  coe-f^  is  a  -func¬ 
tion  o-f  yield.  I-f  y^<  0.2  megatons  then 

c°efi  =  yi*8/(6*CEPi)2 

otherwise. 

coef^  *  y^/(6*CEP^)2  <Re-f  23s 2)  . 

Finally,  i-f  it  is  chosen  to  minimize  the  equivalent  mega- 
tonnage,  the  -following  is  used.  I-f  y^<  1  megaton  then 

coef^  *  y^(2/3) 

otherwise 

coef^  =  <Re-f  23:2). 

To  use  as  many  o-f  the  remaining  warheads  as  possible,  the 
user  must  select  the  set  k  o-f  targets.  Currently*  that  set  k 
o-f  targets  must  be  either  the  -force  targets,  other  military 
targets,  value  targets,  or  any  remaining  targets  where  the  DE 
goals  have  not  been  met.  After  selecting  the  appropriate 
target  set,  constraint  two  is  used  -for  the  extreme  goal. 

This  concludes  the  derivation  o-f  all  the  constraints  that 
can  appear  in  a  problem  solved  by  BRIK.  Some  set  of  these 
constraints  combined  with  an  objective  function  constitute  a 
linear  goal  programming  problem.  Goal  programming,  the 
objective  functions,  the  constraint  sets,  and  the  solution 
algorithm  will  be  presented  in  the  next  chapter. 
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U  Model 

8RIK  is  an  interactive,  goal  programming  mode!  -for  nu¬ 
clear  exchange  problems.  The  first  section  in  this  chapter 
reviews  goal  programming,  providing  an  introduction  to  goal 
programming  for  those  readers  unfamiliar  with  the  topic.  It 
also  offers  a  brief  review  for  those  readers  fluent  in  goal 
programming  techniques.  Following  the  goal  programming  des¬ 
cription  are  the  mathematical  formulations  of  BRIK.  This 
section  includes  the  various  objective  functions  and  the 
constraint  sets  that  form  the  goal  programming  matrix.  This 
chapter  concludes  with  a  description  of  PAGP,  the  preemptive 
goal  programming  algorithm  used  in  BRIK.  As  mentioned  in  the 
preface,  PAGP  was  acquired  from  the  Association  for  Computing 
Machinery  <ACM) .  It  is  copyrighted,  therefore  its  commer¬ 
cial  .use  is  prohibited. 

Gcal  Programming 

Use  of  linear  programming  <LP)  is  restricted  to  a 
single  overriding  objective,  such  as  maximizing  total  profit 
or  minimizing  total  cost.  However,  this  is  not  always  real¬ 
istic.  Sometimes  it  is  desirable  te  conduct  studies  that 
focus  on  a  variety  of  other  objectives,  e.g.,  to  maintain 
stable  profits,  increase  one's  share  of  the  market,  diversify 
products,  or  improve  worker  morale.  An  example  of  these 
statements  that  concerns  many  nuclear  exchange  models  is  the 
case  where  it  is  desirable  to  attain  damage  on  a  certain 
group  of  targets  but  also  to  minimize  the  weapons  used  in  the 
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Allocation.  With  a  standard  LP,  there  are  two  ways  to 
formulate  this  problem.  One,  constraints  could  be  made  re¬ 
flecting  the  desired  damage  expectancy  <DE>  of  the  respective 
target  classes  with  the  objective  of  minimizing  the  amount  of 
weapons  used.  Or,  two,  the  objective  could  maximize  DE,  but, 
to  limit  the  amount  of  weapons  used,  constraints  could  be 
built  which  represent  the  amount  of  weapons  available.  In 
the  first  case,  there  is  the  possibility  of  having  an  infea¬ 
sible  problem.  In  the  second  case,  an  answer  will  be  ob¬ 
tained,  but  it  may  not  be  the  best  answer  because  of  the 
changing  of  the  weapon  constraints  to  limit  warhead  usage. 
This  obstacle  of  multiple  objectives  that  are  incommensurate 
and  often  incompatible  can  be  solved  by  providing  a  method  of 
striving  toward  several  such  objectives  simultaneously.  That 
method  is  called  goal  programming  <S.P.) . 

The  basic  idea  is  to  establish  a  specific  numerical  goal 
for  each  of  the  objectives,  formulate  an  objective  function 
for  each  objective  and  then  seek  a  solution  that  minimizes 


the  weighted  sum  of  deviations  of  these  objective  functions 
from  their  respective  goals  <Ref  14:172).  Mathematically, 
the  goal  programming  model  reduces  to  the  following: 

Min  z  -  r(VikV*Vi5,'di') 


£  — d ^  *  g^  i  ■  1, * • • ,p 

plus  including  any  other  linear  programming  constraints  on 
the  X..  The  variable  g.  represents  each  objective's  goal  and 


Wik*  and  Wis“  provide  the  relative  weight*  of  the  deviationa! 
variables  in  the  kth  or  sth  ranking.  The  use  of  weights 
permits  ranking  of  the  variables  at  each  priority  assuming 
each  variable  is  in  terms  of  a  single  uniform  measure. 

Pk  and  Pg  represent  preemptive-pr iori ty  factors.  Any 
goal  at  preemptive  priority  k  (designated  by  P^>  will  always 
be  preferred  to  (i.e.,  preempt)  any  goal  at  a  lower  priority 
k+l,...,K,  regardless  of  any  scalar  multiplier  W^k+  or  W^g” 
associated  with  these  lower  priorities.  For  example,  consi¬ 
der  the  decision  procedure  in  purchasing  a  home.  The  buyer's 
first  priority  may  be  to  consider  only  a  home  that  is  within 
a  20-mile  radius  from  his  or  her  place  of  work.  All  other 
homes  are  excluded  from  consideration.  Next,  the  buyer  may 
desire  to  limit  the  purchase  price  to  under  *100,000.  Thus, 
even  though  there  may  be  homes  under  *100,000  21  miles  from 
work,  they  are  excluded  (or  preempted)  from  consideration  by 
the  first  priority.  Thus,  the  preemptive-priori ty  concept  is 
used  as  an  iterative  screening  process  (Ref  18:380). 

The  symbols  d^  and  dj*  represent  deviational  variables. 

The  difference  between  what  is  accomplished  and  what  is 
aspired  to  is  the  deviation  from  a  goal .  Having  both  a 
negative  and  a  positive  variable  in  the  formulation  permits 
both  under-  as  well  as  overach i evemen t  of  a  goal.  Sometimes 
it  is  useful  to  include  these  variables,  occasionally  called 
auxiliary  variables,  in  standard  simplex  problems  just  to 
simplify  the  model.  Auxiliary  variables  are  simply  variables 
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that  are  added  to  a  model  -for  convenience  in  addition  to  the 
original  decision  variables  o-f  the  problem.  "In  most  cases, 
this  requires  constructing  an  equation  -for  each  auxiliary 
variable  that  de-fines  this  variable  in  terms  o-f  the  other 
variables  in  the  model.  This  definition  then  is  incorporated 
into  the  model  by  adding  this  equation  as  another  equality 
constraint  in  the  model.  Including  these  equations" .. .al - 
lows  incorporation  o-f  "auxiliary  variables  into  the  objective 
-function  in  a  linear  programming  format*  (Ref  16:176-177). 

In  conclusion,  the  strength  of  goal  programming  comes 
from  its  flexibility.  Many  different  aspects  of  a  constraint 
set  can  be  examined  by  simply  changing  the  priorities  of  the 
objective  terms.  This  will  become  apparent  in  the  following 
section  inhere  only  the  objective  function  is  changed  to  meet 
different  problem  requirements. 

Mathematical  Formulation 

The  feasibility  of  a  model  that  had  a  linearized  damage 
function  with  auxiliary  variables  in  the  target  constraints 
had  already  been  demonstrated  in  an  earlier  effort  by 


Hambsganss  (Ref  27).  However,  to  give  his  model  flexibility, 
his  thesis  presented  many  different  formulations.  To  perform 
an  allocation  to  meet  damage  goals,  auxiliary  variables  from 
the  target  constraints  appeared  in  a  single  objective  func¬ 
tion.  To  complete  the  model,  hard  constraints  were  added  to 
represent  weapon  availability  and  any  hedging  options. 
Mambsganss  cautioned  that  too  many  hedges  may  create  a  prob- 
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1  cm  with  no  -feasible  solutions  (Ref  27s 77).  In  another 
-formulation,  to  compute  a  possible  reserve  force,  hard  con¬ 
straints  were  made  out  of  both  the  target  and  weapon  con¬ 
straints.  The  terms  in  the  objective  function  consisted  of 
all  of  the  decision  variables.  To  preclude  an  allocation 
consisting  of  entirely  one  weapon,  or  the  greater  portion  of 
a  weapon  class,  Wambsganss  suggested  reducing  the  weapon 
availability  prior  to  performing  the  allocation  (Ref  27:80). 
Again,  the  possibility  of  an  infeasible  probtem  exists. 

These  problems  of  infeasibility  plus  the  programming  complex¬ 
ity  resulting  from  occasionally  adding  auxiliary  variables 
motivated  a  search  for  a  different  formulation.  The  formula- 

i 

tion  eventually  agreed  upon  used  goal  programming  and  appears 

! 

♦ 

in  Figure  2. 

Three  advantages  of  goal  programming  are  immediately 

! 

apparent.  One,  the  form  of  the  constraint  set  does  not 
change.  Regardless  of  the  problem  being  analyzed,  every 


constraint  contains  deviational  variables.  The 


advantages 


gained  from  this  standardization  were  realized  during  the 
flowcharting  of  BRIK.  Knowledge  of  the  objective  function  is 
not  needed  to  keep  track  of  which  columns  hold  which  decision 
variables.  The  second  advantage  of  goal  programming  came 
from  the  objective  function.  All  of  the  different  alloca¬ 
tions  can  be  performed  with  just  the  deviational  variables. 

It  is  now  possible  to  minimize  the  weapons  used  in  a  problem 
without  having  the  decision  variables  appearing  in  the 
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Objective  Functions 

1) 
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Any  additional  hedging  constraints 

Figure  2.  Mathematical  Formulations  of  BRIK 


objective  function.  This  also  simplified  programming. 
Finally,  the  third  advantage  gained  from  6.P.  is  flexibility. 
By  just  varying  the  priorities  of  the  deviational  variables, 
many  different  questions  can  be  answered  concerning  the  for¬ 
mulation.  This  flexibility  will  be  better  understood  given  a 
description  of  each  of  the  three  objective  functions. 

The  first  objective  function  is  designed  to  take  the 
available  arsenal  and  allocate  it  to  meet  the  DE  goals  as 
best  possible.  It  has  four  parts  and  can  use  as  many  as  ten 
preemptive  priorities.  Part  one  bounds  the  problem  with  the 
available  arsenal.  At  priority  one,  PI,  the  positive  devia¬ 
tional  variables  in  the  weapon  constraints  are  minimized. 

This  will  prevent  any  lower  goal  from  using  more  weapons  than 
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are  available.  In  the  second  part,  all  of  the  user-def ined 
hedging  options  are  met  at  priority  two,  P2.  With  G.P.,  i-f  a 
particular  goal  is  not  completely  satisfied,  one  of  the  two 
deviational  variables  will  have  some  positive  value.  This 
fact  makes  it  easy  to  determine  if  there  is  an  insufficient 
arsenal  to  meet  the  hedges.  Should  there  be  an  insufficient 
arsenal,  some  deviational  variables  in  the  hedging  con¬ 
straints  will  have  a  positive  value  in  the  solution.  In  part 
three  of  objective  function  one,  weapons  are  allocated  to  the 
target  classes  in  order  of  priority.  BRIK  permits  division 
of  the  target  classes  into  as  many  as  seven  preemptive  prio¬ 
rities.  This  is  done  by  choice  of  the  user.  If  it  is  chosen 
not  to  separate  the  classes  into  priorities,  all  of  the 
target  goals  will  appear  at  priority  three,  P3. 

At  this  point,  it  is  convenient  to  discuss  the  weights 
put  on  the  deviational  variables  for  the  targeting  goals.  As 
mentioned  earlier,  the  use  of  weights  permits  ranking  of  the 
variables  at  each  priority.  This  is  normally  accomplished  by 
multiplying  the  deviational  variable  in  the  objective  func¬ 
tion  by  some  value  associated  with  a  particular  goal.  How¬ 
ever,  this  process  has  been  complicated  by  the  linearization 
of  the  constraints  (see  Chapter  IV> .  Before  the  lineariza¬ 
tion,  a  target  constraint  had  the  following  form: 

pl/i^nj-Pj  ■  PVr 

As  mentioned  earlier,  minimizing  p*  would  drive  the  con- 

straint  to  the  survival  goal  of  PV.t  .  After  linearizing  this 

J 
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constraint!  it  was  shown  that 

tTjlnd*^) 

can  play  the  role  of  the  negative  deviational  variable  d.”. 

,  J 

Thereforey  to  drive  the  linearized  constraint  to  its  goal,  it 


is  necessary  to  now  minimize  d  •  ,  because  this  term  holds  the 

J 

variable  p^.  If  the  constraints  had  not  been  linearized,  the 
targets  could  be  weighted  within  each  priority  by  multiplying 
the  deviational  variable  Pj  by  a  number  representing  target 
class  value  •  However,  the  deviational  variable  in 

each  constraint  is  not  p^,  it  is 

V  * 

If  this  equality  was  solved  for  pj ,  the  result  would  be  the 
following: 

pj  *  pVj(e(dr/^J)-l). 

Multiplying  both  sides  by  yields 

vjNjpj  *  vjWe(dj  <i> 

Therefore,  since  dj  appears  in  each  target  constraint,  to 
properly  weight  the  deviational  variables,  the  right  side  of 
equation  one  should  be  in  the  objective  function  with  the 


goal  of  minimizing  d*  .  This  is  not  solvable  with  LP! 

J 

Therefore,  the  following  approximation  was  used  in  BRIK. 

V/nj 

gets  smal 1 , 

(e^r^V-D  *  d^/Nj 

yiel ding 


V.N.p.  *  V.PV.d 

J  J  J  J  V  J 


As 
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Thus,  the  number  used  to  weight  the  negative  deviations]  var¬ 
iables  in  the  targe’t  constraint  in  objective  function  one  is 

vi  ■  YY 

The  last  part  of  objective  one  and  the  lowest  priority 
goal  is  the  extreme  goal.  The  purpose  of  this  goal  is  to 
minimize  some  aspect  of  the  allocation.  Because  of  the' 
preemptive  features  of  G.P.,  this  goal  may  or  may  not  be 
considered.  If  some  of  the  DE  goals  cannot  be  met  due  to  an 
insufficient  arsenal,  the  extreme  goal  will  not  be  con¬ 
sidered.  However,  if  all  of  the  goals  are  met,  there  will  be 
weapons  left  over.  If  that  is  the  case,  because  of  the  use 
of  non-integer  programming,  virtually  an  infinite  combination 
of  allocations  could  be  reported.  Thus,  the  extreme  goal 
serves  the  role  of  driving  the  allocation  to  a  unique  solu¬ 
tion.  This  first  objective  function  performs  an  allocation 
restricted  to  weapon  availability.  Thus,  depending  on  the 
size  of  the  arsenal,  the  goals  may  or  may  not  be  met.  This 
contrasts  with  the  second  objective  function  where  the  OE 
goals  will  be  met. 

The  second  objective  function  gives  the  user  the  ability 
to  build  an  arsenal  to  meet  all  DE  goals.  This  is  accom¬ 
plished  with  five  preemptive  priorities.  At  PI,  the  positive 
deviational  variable  in  each  target  constraint  is  minimized, 
guaranteeing  that  the  targeting  goals  will  not  be  exceeded. 

At  P2,  all  of  the  hedging  constraints  are  met  and  at  P3  all 
of  the  targeting  goals  are  met.  At  this  point,  the  user  has 
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two  options  concerning  the  weapon  base.  The  weapon  classes 
can  be  input  where  the  availability  is  zero  or  some  number 
can  be  input  that  represents  the  "minimum*  numbers  of  each 
weapon  type  that  should  be  used  in  the  allocation.  Since 
weapon  availability  is  unbounded,  there  is  no  limit  to  the 
number  of  weapons  from  each  class  that  can  be  allocated.  So, 
regardless  of  how  many  weapons  the  user  inputs  for  the  mini¬ 
mum  arsenal,  all  DE  goals  will  be  met.  However,  should  some 
minimum  arsenal  be  entered,  this  objective  function  will 
attempt  to  use  all  of  those  weapons  prior  to  "building" 
any  additional  warheads.  {This  is  accomplished  at  P4  by 

i 

minimizing  the  negative  deviational  variables  on  the  weapon 
constraints.  There  is  a  side-effect  for  this  priority  of 

which  the  user  must  be  aware.  Should  the  minimum  arsenal  be 

j 

sufficient  to  meet  all  of  the  DE  goals,  P4  will  have  the 
effect  of  maximizing  the  number  of  warheads  used  and  the 
following  priority,  PS,  will  have  no  effect.  Finally,  just 

like  ths  lowest  priority  goal  in  the  first  objective  func- 

{ 

tion,  PS  insures  some  minimum  aspect  of  the  used  arsenal. 

The  final  abjective  function,  type  three,  has  four  parts 
and  requires  four  preemptive  priorities.  Unlike  the  first 
two  objectives  which  bound  the  problem  with  either  weapon 
availability  or  DE  goals,  this  third  objective  function 
bounds  the  problem  with  both.  At  Pi,  the  positive  devia¬ 
tional  variables  of  the  weapon  constraints  are  minimized. 
Also,  at  P2,  the  positive  deviational  variables  of  the  target 
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constraints  are  minimized.  This  guarantees  that  neither  the 
weapon  availabilities  nor  the  designated  goals  will  be  ex¬ 
ceeded.  At  P3,  all  oF  the  hedging  options  will  be  met  within 
the  limits  set  by  PI  and  P2.  Finally,  the  extreme  goal  will 
be  considered  at  P4.  IF  this  objective  Function  is  chosen, 
there  are  two  major  diFFerences  From  the  First  two  objective 
Functions  oF  which  the  user  must  be  aware.  First,  the  DE 
goals  are  not  used  to  drive  the  allocation.  They  are  used 
only  as  upper  bounds  to  the  problem.  Second,  the  extreme 
goal  has  a  FiFth  option-  that  is  not  present  in  the  other  two 
objectives.  That  option  permits  dumping  any  remaining  arse¬ 
nal  on  targets  whose  goals  are  riot  yet  met.  Finally,  the 
type  three  objective  Function  relies  heavily  on  the  hedging 
constraints  to  drive  the  allocation.  ThereFore,  it  is  esti¬ 
mated  that  it  would  take  longer  to  set  up  a  problem  using  the 
type  three  Function  as  opposed  to  the  type  one  or  the  type 
two  objective. 

In  conclusion,  Q.P.  provides  a  tremendous  amount  oF  Flex¬ 
ibility.  So  much  in  Fact,  it  could  be  argued  that  restricting 
the  model  to  only  three  objective  Functions  restricts  some  oF 
the  gain  From  G.P.  It  could  have  been  possible  to  build  BRIK 
such  that  the  user  could  stipulate  what  variables  From  what 
constraints  go  into  what  priorities.  That  indeed  would  have 
permitted  the  user  to  take  Full  advantage  oF  the  mathematics. 
However,  it  was  decided  not  to  allow  the  ability  to  manually 
juggle  objective  terms  because  oF  our  stated  purpose  oF 
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developing  a  user-friendly  model.  A  user  does  not  require  an 
in-depth  knowledge  of  G.P.  to  use  BRIK.  To  select  one  of  the 
above  objective  functions,  ill  that  is  needed  is  the  ability 
to  choose  between  1,  2,  or  3.  Everything  is  then  built 
internally. 

O.P.  has  been  described  and  the  formulations  have  been 
presented.  The  last  part  of  this  chapter  concerns  the  algo¬ 
rithm  used  to  solve  the  problem  just  described.  That  algo¬ 
rithm  is  called  a  Partitioning  Algorithm  for  Goal  Programming 
Problems  <PAGP) . 

PAGP 

BRIK's  evolution  into  a  G.P.  model  occurred  simultaneous¬ 
ly  with  the  search  for  an  independent  allocation  routine.  To 
increase  the  model's  portability,  it  was  concluded  early  in 
the  project  that  BRIK  would  not  be  fettered  by  any  local 
library  package.  At  first,  time  was  spent  looking  for  an 
LP  package  that  could  be  used  as  a  subroutine.  But  as  the 
current  model  evolved,  the  search  expanded  to  include  goal 
programming  routines.  Mhile  it  was  not  beyond  the  scope  of 
this  effort  to  write  a  G.P.  subroutine,  finding  an  already 
existing  package  helped  simplify  the  coding  problem.  It  was 
decided  to  use  PAGP  because,  besides  meeting  all  of  the 
requirements  for  portability,  the  program  had  coding  effi¬ 
ciencies  of  which  BRIK  could  take  advantage.  These  efficien¬ 
cies  include: 

1)  Partitioning  of  the  goal  constraints. 
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2)  Variable  elimination. 

3)  Special  termination  rules. 

This  section  will  describe  each  of  these  three  coding 
"tricks"  included  in  PAGP  and  explain  why  they  are  an  advan¬ 
tage  to  BRIK. 

The  -first  advantage  o-f  PAGP  is  the  partitioning  proce¬ 
dure.  The  ordinal  priority  -factors  in  the  goal  programming 
objective  -function  are  used  to  partition  the  goal  constraints 
o-f  the  problem.  "This  is  accomplished  -first  by  observing 
that  -for  any  goal  constraint  i,  one  and  only  one  f  three 
things  may  occur: 

1)  only  d^~  appears  in  the  objective  function, 

2)  only  d^b  appears  in  the  objective  function,  or 

3>  both  d^“  and  d^+  appear  in  the  objective  function. 

In  case  < 1) ,  the  partition  would  assign  goal  constraint  i  to 
the  priority  factor  associated  with  d^”f  in  case  <2> ,  con¬ 
straint  i  would  be  assigned  to  the  priority  factor  associated 
with  dj* j  while  in  case  <3> ,  the  partition  would  determine  the 
higher  order  priority  factor  <in  terms  of  the  ordinal  rank¬ 
ing)  associated  with  either  d^~  or  d^+  and  constraint  i  would 
be  assigned  to  that  priority"  <Ref  2:379).  Only  the  con¬ 
straints  that  had  deviational  variables  in  PI  are  in  the 
initial  problem.  Next,  the  only  new  constraints  added  at  P2 
are  those  constraints  that  had  deviational  variables  at  P2 
and  not  at  PI.  This  continues  until  ei ther  all  of  the  con¬ 
straints  in  the  problem  are  added  or  the  stopping  rule  is 
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invoked.  BRIK  takes  advantage  of  partitioning  because  dif¬ 
ferent  sets  of  constraints  are  brought  in  at  different  prior¬ 
ities.  Consider  objective  one.  Only  the  weapon  constraints 
are  in  the  problem  at  Pi,  the  hedging  constraints  are  brought 
in  at  P2,  the  target  constraints  are  added  at  P3,  and  finally 
the  extreme  goal  is  input  at  P4.  PAGP  is  actually  solving  a 
series  of  smaller  subproblems  in  order  to  find  a  solution  to 
the  original  problem.  This  increases  efficiency. 

The  second  "trick"  in  PAGP  is  the  elimination  procedure. 
Not  only  does  it  help  BRIK,  it  improves  efficiency  in 
any  problem  solved  by  PAGP.  The  motivation  behind  the  eli¬ 
mination  procedure  comes  from  the  theory  of  L.P.  If  Z  is  the 
optimal  value  of  an  L.P.  problem  and  <Zj-Cj>>0  for  some  non- 
basic  Xj,  Xj  cannot  enter  the  basis  to  form  an  alternative 
optimal  solution.  A  corollary  to  this  statement  that  applies 
to  G.P.  concerns  the  optimal  solution  to  a  subproblem  S^. 

Any  nonbasic  variable  which  has  at  least  one  positive  rela¬ 
tive  cost  <z.-c.>0>  can  be  eliminated  from  entering  the  basis 
J  J 

in  subproblem  S  ,  where  y  ■  k+l,...,P  <Ref  2x380).  This 
¥ 

elimination  makes  the  program  more  efficient  because  fewer 
nonbasic  elements  are  considered  to  enter  the  basis  at  each 
pivot.  PAGP  accomplishes  this  by  maintaining  an  indicator 

l 

row,  INO.  While  this  indicator  row  increases  the  efficiency 
of  PAGP,  it  should  be  noted  that  PAGP  still  pivots 
columns.  The  program  could  be  further  improved  if 
to  pivot  the  eliminated  columns. 
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The  third  "trick"  used  in  PAGP  is  a  special  stopping 
rule.  This  is  accomplished  by  using  the  vector  IND.  If  it 
indicates  no  nonbasic  elements  can  enter  the  basis,  the 
routine  stops,  regardless  of  how  many  lower  priorities  are 
left  to  be  considered.  •  BRIK's  extreme  goal  takes  advantage 
of  this  last  efficiency  item.  In  objective  one,  the  DE  goals 
are  met  as  best  possible  at  P3.  At  P4,  if  all  of  the  DE 
goals  are  attained,  some  aspect  of  the  allocation  is  mini¬ 
mized.  However,  if  some  DE  goals  were  not  met,  PAGP's  stop¬ 
ping  rule  prevents  the  constraint  at  PA  from  ever  being  added 
to  the  G.P.  matrix.  This  results  in  added  efficiency. 
Conclusion 

This  concludes  the  mathematical  presentation  of  BRIK.  It 
is  a  goal  programming  model  that  has  its  own  independent 
allocation  routine.  The  User  Guide  is  in  Appendix  A  and  a 
listing  appears  in  Appendix  E.  The  following  chapter  con¬ 
cerns  verification  and  demonstration. 
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This  chapter  will  discuss  verification  of  BRIK,  then 
demonstrate  some  of  its  capability.  As  defined  by  the  Random 
House  Dictionary  of  the  English  Language,  to  verify  means  to 
‘ascertain  correctness  of,  as  by  examination,  research,  or 
comparison."  This  section  examines  the  steps  taken  to  verify 
the  correctness  of  this  model.  The  demonstration  section 
will  start  with  a  sample  problem  that  compares  BRIK  solutions 
with  output  from  AEM.  Then,  sample  results  used  to  illus¬ 
trate  BRIK's  objective  functions  and  hedging  options  will'-be 
displayed. 


BRIK  has  3723  lines  of  code  divided  among  49  subroutines. 
This  section  will  describe  the  effort  undertaken  to  ascertain 
its  correctness.  Also,  this  section  indirectly  provides  a 


more  detailed  description  of  the  physical  composition  of  the 
model  by  reviewing  the  functions  of  the  subroutines  and  the 
use  of  the  data  files  that  are  built.  While  most  small 
algorithms  can  be  verified  with  hand  calculations,  verifica¬ 
tion  of  a  program  of  this  size  must  use  some  type  of  systema¬ 
tic  approach.  The  approach  used  divided  the  model  into  six 
sections! 


1)  Editor 


2)  SSPS  Computation 

3)  Matrix  Formulation 

4>  Objective  Development 


5)  PAGP  (Assignment  Algorithm) 

6)  Output  ' 

The  steps  taken  to  insure  verification  of  each  section  will 
be  presented  separately. 

Edi tor .  The  editor  has  two  parts,  the  target  editor  and 
the  weapon  editor.  Except  for  the  fact  that  one  is  for 
target  data  and  the  other  is  for  weapon  data,  they  are  essen¬ 
tially  the  same.  The  verification  process  consisted  of  exer¬ 
cising  the  program  to  determine  if  both  sections  were  acting 
as  expected.  For  example,  both  methods  of  data  entry,  either 
through  files  or  interactively,  were  exercised  until  it  was 
visually  determined  that  information  was  going  into  the  pro¬ 
per  storage  locations.  This  visual  capability  was  made  pos¬ 
sible  by  the  subroutines  TGTEDT  and  WPNEDT.  Each  of  these 
subroutines  permit  a  viewing  of  the  data  on  the  monitor. 

After  it  was  determined  that  weapon  and  target  data  was 
being  entered  properly,  effort  was  concentrated  on  the  trans¬ 
portation  between  subroutines  called  by  TGTEDT  and  WPNEDT. 
When  outlining  the  flowchart  for  this  model,  it  was  deter¬ 
mined  that  it  was  easiest  to  control  transportat i on  if  the 
user  was  returned  to  a  main  menu  after  executing  a  parti¬ 
cular  option.  This  technique  proved  invaluable,  especially 
since  the  first  option  or  each  menu  permitted  the  user  to 
view  the  data  base.  Verification  again  consisted  of  exerci¬ 
sing  the  program.  After  performing  the  different  menu 
options,  a  visual  inspection  of  the  data  was  made  to  insure 
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that  the  desired  results  were  attained.  Eventually,  it  was 
possible  to  move  throughout  both  edi tors  and  always  return  to 
their  respective  main  menu.  That  eventual  return  coupled 
with  a  correct  data  display  insured  validity  of  the  editors. 

SSPS.  When  BRIK  was  in  its  earliest  stages  of  develop¬ 
ment,  it  was  visualized  to  have  two  separate  programs.  The 
-first  program  would  contain  the  editor  and  the  matrix  generator 
and  the  second  program  would  have  the  allocation  routine  and 
the  output.  The  intention  was  to 'parallel  the  BLUE  SWADE  model 
at  SAC  Headquarters.  To  -facilitate  these  two  programs,  files 
were  built  by  the  first  program  that  would  be  used  as  input 
for  the  second  program.  One  of  those  files  is  called  SPARSE. 

SPARSE  is  a  j  by  i  matrix  containing  the  single  shot 
probability  of  survival  <SSPS)  of  a  target  in  class  j  if 
assigned  a  warhead  from  a  weapon  in  class  i.  There  are  two 
different  sources  for  each  SSPS  number  calculated.  If  the 
user  inputs  the  target  vul nerabi 1 i ty  in  terms  of  PSI ,  proba¬ 
bility  of  kill  <PK)  is  calculated  by  function  PK.  If  i  VNTK 
number  is  input  to  express  target  vul nerabi 1 i ty ,  PK  is  com¬ 
puted  by  Function  VTK.  After  a  PK  is  computed,  SSPS  is  at¬ 
tained  by  subtracting  the  product  of  PK  and  weapon  reliabil- 
ity  from  one.  Each  of  the  two  functions,  PK  and  VTK,  were 
verified  separately. 

Function  PK  was  verified  by  following  a  series  of  inputs 
with  hand  calculations.  These  calculations  matched  the  en¬ 
tries  of  SPARSE.  Therefore,  PK  was  assumed  to  be  correct. 

72 


X  -X 

HI  /•' 


7V;V 


-X- 

J- 


A 


This  relatively  simple  process  contrasts  greatly  with  the 
steps  taken  to  verify  Function  VTK. 

VTK  is  a  function  extracted  directly  from  BALLOT.  Since 
BALLOT  has  enjoyed  extensive  use  at  SAC  Headquarters,  it  was 
assumed  that  the  output  would  be  acceptable  if  the  function 
was  being  used  properly.  First,  it  was  necessary  to  determine 
the  units  of  the  input  variables.  This  was  gleaned  from  the 
program  listing.  In  BALLOT,  CEP  is  in  feet,  Yield  is  in 
kilotons,  and  R-95  is  in  nautical'  miles.  Second,  sample 
calculations  were  made  with  the  function  to  insure  reason¬ 
ableness.  When  the  authors  were  satisfied  that  the  function 
was  being  properly  used,  tables  similar  to  Figures  3-A  and  3-B 
were  computed  for  both  T  factors,  P  and  Q. 

The  intuitive  consistency  of  these  tables  is  sufficient 
assurance  of  the  verification  of  Function  VTK.  In  the 
center  row,  as  yield  increases,  the  PK  increases.  In  the 
center  column,  as  CEP  increases,  PK  decreases.  Both  of  these 
make  sense.  Also,  notice  that  the  PK  increases  as  K  in¬ 
creases.  The  K  factor  allows  for  hardness  adjustments  to  be 
made  to  account  for  the  effects  of  variations  in  blast  wave 
duration  due  to  different  weapon  yields  (Ref  5:34>.  Since 
targets  with  a  high  K  factor  are  more  susceptible  to  damage 
from  shock  wave  duration  than  targets  with  a  low  K  factor,  it 
makes  sense  that  as  the  K  factor  increases,  so  does  PK. 

To  further  verify  Function  VTK,  a  comparison  was  made 
between  PK' s  computed  by  BRIK  and  some  PK's  computed  by  AEM. 
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Figure  3-A.  Various  Probabi 1 i ty  of  Kill 
<R  95  =>  0  and  T  =  P) 


The  first  sample  problem  presented  in  the  demonstration 
section  of  this  chapter  contains  five  weapon  classes  and  ten 
target  classes.  To  solve  this  problem,  it  was  necessary  for 
both  programs  to  compute  50  different  PK's,  one  for  each 
weapon/target  combination.  A  visual  comparison  of  those  two 
sets  of  numbers  revealed  that  Function  VTK's  numbers  matched 
AEM's  probabi 1 i ty  of  kill  computations  to  the  third,  and 
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Figure  3-B.  Various  Probabi 1 i ty  of  Kill 
<R  95  -  0  and  T  -  Q> 


sometimes  the  fourth  decimal  place.  This  further  verifies 
that  the  borrowed  subroutine  is  being  used  correctly. 

In  conclusion,  verification  of  the  SSPS  computations  for 

I 

PSI  inputs  came  from  hand  calculations  and  verification  of 
the  UNTK  system  calculations  came  from  a  comparison  with  an 
established  nuclear  exchange  model. 

Matrix  Formulation.  AIJ  is  a  file  that  contains  the 


matrix  of  the  current  problem.  This  file  was  built  primarily 


to  provide  the  output  section  right-hand-side  values.  Its 
secondary  purpose  was  to  simplify  the  verification  process. 
Small  problems  were  run  using  all  of  the  hedges  and  each  of 
the  five  types  of  extreme  goals.  The  AIJ  files  generated  for 
each  run  were  examined  for  verification. 

There  were  two  examinations  the  files  had  to  pass,  the 
numbers  in  each  constraint  had  to  be  in  the  proper  columns 
and  they  had  to  be  accurate.  A  visual  inspection  verified 
that  the  numbers  were  in  the  right  columns  for  the  small 
sample  problems.  Since  the  numbers  are  properly  located  for 
small  problems,  they  will  be  properly  positioned  for  large 
problems  because  the  algorithm  which  is  used  for  column 
placement — 

£1  ♦  j  +  <i-l>*NTGTS3 

where  i  is  the  weapon  class  number,  j  is  the  target  class 
number,  and  NTGTS  is  the  number  of  target  classes  in  the 
current  problem — is  dependent  on  problem  size.  After 
verifying  that  the  numbers  were  in  their  correct  positio.  s, 
it  was  necessary  to  check  their  correctness.  The  coef¬ 
ficients  in  the  target  constraints  are  the  negative  loga¬ 
rithms  of  the  SSPS's  found  in  SPARSE.  It  was  a  simple  matter 
to  verify  those  numbers  having  both  SPARSE  and  AIJ  available. 
Right  hand  sides,  EMT,  and  CMP  figures  were  all  checked  with 
hand  calculations. 

Objective  Development.  PAGPIN  is  a  file  used  by  the  goal 
progamming  package.  It  contains  all  of  the  information  PAGP 
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requires  to  solve  the  goal  programming  problem.  Since  all  of 
the  data  describing  the  objective  functions  is  in  this  file, 
verification  consisted  of  visually  checking  its  accuracy.  It 
must  be  noted  that  the  information  contained  in  PAGPIN  is  the 
authors'  i nterpretat i on  of  PAGP's  requirements.  This  inter¬ 
pretation  resulted  from  reading  the  program  listing  acquired 
from  ACM.  The  verification  of  that  interpretation  will  be 
discussed  in  the  next  section. 

PAGP.  Two  things  were  required  to  verify  PAGP,  verifica¬ 
tion  of  the  algorithm  and  verification  of  the  interpretation 
of  its  requirements.  It  was  decided  that  this  could  be 
accompl i shed  by  using  it  to  run  text  book  examples.  Despite 
PAGP  being  proofed  prior  to  publication  and  the  test  examples 
coming  from  a  text  book,  it  is  presumptuous  to  think  that 
both  could  be  without  error.  However,  because  they  are 
independent,  the  chance  of  achieving  the  same  answer  given 
that  one  or  both  are  in  error  is  considered  negligible. 

Thus,  the  examples  on  pages  394,  404,  and  408  from  Ignizio 
<Ref  18)  were  run.  PAGP  computed  the  same  answers  reported 
in  the  text.  Therefore,  since  the  output  compared  favorably, 
three  statements  could  be  made! 

1)  The  program  PAGP  was  operating  properly. 

2)  The  text  book  examples  were  correct. 

3>  Our  interpretation  of  the  needs  of  PAGP  was  accurate. 

Output .  The  output  routine  uses  the  results  from  the 
allocation  and  computes  the  damage  achievement.  This 
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achievement  is  displayed  on  the  monitor  screen  along  with  the 
target  classes  attacked,  the  weapon  systems  that  attacked 
them,  the  number  of  weapons  assigned  to  each  target  class,  and 
the  expected  number  o-f  targets  remaining  in  each  target  class. 
Since  the  output  section  uses  the  number  o-f  weapons  allocated 
to  compute  damage  achievement,  veri-fying  the  accuracy  o-f  the 
achievement  would  insure  verification  of  the  output  section. 

Verification  of  the  achievement  is  as  follows.  At  this 
point,  PA6P  is  working  properly.  Therefore,  if  the  problem  is 
formulated  correctly  and  there  are  sufficient  weapons  in  the 
inventory  to  meet  the  OE  goals,  the  reported  damage  achieve¬ 
ment  from  an  allocation  should  equal  the  goals.  Example 
problems  were  run  and  the  achievement  did  equal  the  goals. 
Therefore,  the  output  section  is  verified. 

This  concludes  the  formal  process  of  verifying  BRIK.  It 
is  allocating  weapons  to  targets  exactly  as  it  is  being  asked. 
Demonstration 

This  demonstration  section  consists  of  two  parts.  In  the 
first  part,  a  test  problem  is  presented  where  a  solution 
generated  by  AEM  is  compared  to  solutions  from  BRIK.  The 
second  part  runs  through  a  series  of  small  examples  to  show 
the  capability  and  versatility  of  BRIK. 

Test  Problem.  BRIK  allocates  weapons  to  targets  only  if 
the  DE  goals  are  known.  Therefore,  to  compare  BRIK  with  the 
Arsenal  Exchange  Model,  the  following  problem  was  designed.  A 
target  base  and  a  weapon  base  were  supplied  by  the  authors  to 
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USAF/SA.  These  two  bases  were  run  on  AEM  with  the  objective 
to  maximize  DE.  There  were  no  special  allocation  rules  or 
hedges  invoked  in  the  problem.  After  receiving  AEM's 
allocation,  the  OE  attained  on  each  target  class  became  the  DE 
goals  -for  BRIK.  This  made  it  possible  to  -form  a  comparison 
between  the  two  models.  The  -following  is  an  in-depth 
description  of  the  entire  process  along  with  the  results. 

The  problem  is  to  allocate  a  weapon  base,  consisting  of 
five  classes,  to  a  target  base  consisting  of  ten  classes 
(Table  VI)  in  order  to  maximize  damage  expectancy  <DE> . 


Table  VI 
Test  Problem 
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NAME 

NUMB  VUL 

DIM 

TYPE 

PRI0 

VAL 

REL 

CEP 

YLD  UHDS 

civil 

140  24Q0 

.56 

B 

1 

1.00 

.00 

.00 

.00 

0 

local 

215  13P1 

.49 

B 

1 

1.00 

.00 

.00 

.00 

0 

c3i 

450  35P7 

.00 

B 

1 

1.00 

.00 

.00 

.00 

0 

icb» 

1000  52P8 

.00 

f 

1 

1.00 

.00 

.00 

.00 

0 

lcc 

200  39P0 

.00 

f 

1 

1.00 

.00 

,00 

.00 

0 

subpts 

20  22P1 

.36 

f 

1 

1.00 

.00 

.00 

,00 

0 

irbis 

150  11P0 

.00 

f 

1 

1.00 

.00 

.00 

.00 

0 

afbase 

100  10Q1 

.79 

f 

1 

1.00 

.00 

,00 

.00 

0 

stores 

430  31P6 

.00 

m 

1 

1.00 

,00 

,00 

.00  _ 

0 _ 

facil 

wr 

520  23Q0 

.31 

A 

1 

1.00 

.00 

,00 

.00 

0 

sJ 

NAME  NUMWPNS  WHD/WPN 

REL 

CEP 

YLD 

DrtYALRT 

GENALRT 

MMII 

450 

1 

.85 

.200 

2.00 

.85 

1,00 

MMIII1 

250 

3 

.85 

.150 

.17 

.85 

1.00 

Poseid 

104 

10 

.85 

.240 

.05 

.85 

1,00 

B52grv 

106 

4 

.85 

.600 

.35 

.33 

1.00 

FB111 

60 

6 

.85 

.200 

.20 

.33 

1,00 
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The  -first  half  of  Table  VI  consists  of  the  number  of  target 
classes  in  the  base,  the  header  row,  and  ten  rows  that  fully 
describe  each  target  class.  Column  one,  NAME,  gives  the  name 
of  each  class,  and  column  two,  NUMB,  gives  the  number  of 
targets  in  each  class.  Column  three,  VUL,  is  the  target 
vulnerability  index.  In  this  problem  WTK  numbers  are  used. 
Column  four,  DIA,  gives  the  diameter  of  the  smallest  circle 
that  encompasses  95%  of  the  target.  Columns  five  and  six, 
TYPE  and  PRIO,  indicate  the  type  of  target  class  and  each 
class'  priority.  In  this  problem  all  target  classes  are 
priority  1.  Finally,  column  seven,  UAL,  shows  the  value  of 
each  target  in  the  class.  All  the  targets  in  this  problem 
have  a  value  of  one.  The  last  four  columns,  REL,  CEP,  YLD, 
and  WHDS  are  used  to  calculate  'force  target”  values  and  are 
not  used  in  this  problem.  Therefore,  the  input  numbers  are 
zero. 

The  next  section  of  Table  VI  contains  the  weapon  base. 

The  first  two  rows  are  similar  to  the  target  base.  They 
contain  the  number  of  weapon  classes  in  the  base  and  the 
header  row.  For  this  problem,  there  are  five  weapon  classes. 
Each  weapon  class  name  and  the  number  of  weapon  systems 
available  in  each  class  are  listed  in  columns  one  and  two, 
NAME  and  NUMWPNS.  The  third  column,  WHD/WPN ,  indicates  the 
number  of  warheads  each  weapon  system  contains.  The  fourth 
column,  REL,  is  the  probability  that  each  warhead  will  arrive 
and  detonate,  while  the  fifth  column  is  the  weapon  system  CEP 

80 


in  nautical  miles.  The  sixth  column,  YLD,  is  the  yield  in 
megatons  of  each  warhead.  The  final  two  columns,  DAYALRT  and 
GENALRT,  contain  the  percentage  of  weapon  systems  on  alert. 

AEM's  solution  is  given  in  Table  UI I .  The  target  class 
name  is  listed  first.  The  next  two  columns,  value  and  type, 
are  self-explanatory.  The  fourth  column,  no.,  is  the  number 
of  targets  in  the  class.  Column  six  is  the  number  of  targets 
that  survived  the  allocation.  Columns  seven  and  eight  are  not 
used.  Column  nine  provides  the  expected  percentage  of  kill  in 
each  target  class.  If  more  than  one  weapon  type  attacked  the 
target  class,  the  total  expected  kill  percentage  is  in  the  row 
marked  "sum".  Columns  10,  11,  and  12  give  the  allocation 
information.  For  example,  with  target  class  "civil*  in  the 
top  row,  20  targets  (column  10)  were  each  attacked  by  1 
warhead  (column  11)  of  weapon  type  fblll  (column  12). 

Although  AEM's  total  value  destroyed  is  not  given  in  the 
table,  AEM's  score  was  1770. 

Since  BRIK  allocates  weapons  to  targets  to  meet  DE  goals, 
the  DE  achieved  by  AEM  for  each  target  class  was  used  as  input 
for  BRIK.  Objective  function  one  and  extreme  goal  one  (meet 
goals  with  the  fewest  amount  of  weapons)  were  used  in  the 
allocation.  The  output  is  listed  in  Tab’s  'vi  1 1 .  The  first 
two  columns  list  the  target  class  name  and  the  weapon  types 
allocated  to  it.  The  third  column  is  the  number  of  weapons  of 
the  weapon  class  allocated  to  the  target  class.  The  fourth 
column  is  the  DE  goal.  The  fifth  column  is  the  percentage  of 
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the  target  class  destroyed  in  the  allocation.  It  is  not  the 
percentage  o-f  the  goal  met!  The  final  column  is  the  number  of 
targets  remaining  after  the  allocation. 

There  are  many  similarities  between  the  two  allocations. 
The  total  value  destroyed  was  the  same  and  many  of  the  target 
classes  received  the  same  type  and  amount  of  weapon  types. 
However,  there  are  two  discrepancies.  BRIK  indicated  that 
there  were  unused  weapons.  Also,  four  target  classes  were  not 
assigned  the  same  type  or  number  of  weapons.  Because  of  these 
discrepancies,  several  iterative  runs  were  made,  gradually 
increasing  the  DE  goals  until  the  remaining  weapons  were  used. 
The  final  result  is  listed  in  Table  IX. 

The  most  notable  exception  between  BRIK's  final  run  in 
Table  IX  and  AEM's  solution  is  the  increase  in  value 
destroyed,  1788.  This  is  due  to  BRIK  reporting  a  slightly 
higher  damage  achievement  than  AEM  for  some  of  the  target 
classes.  For  example,  BRIK  reports  a  DE  of  .75  for  ”C3i" 
while  AEM  reported  .733.  This  higher  DE  occurs  approx imatel y 
four  times  in  BRIK's  answer  and  is  caused  by  the  damage 
function  approximation  described  in  Chapter  III. 

Even  though  BRIK's  output  does  not  match  AEM's  solution 
exactly,  they  are  very  close.  Using  AEM' s  solution  as  the  DE 
goals,  a  difference  of  only  6X  in  the  amount  of  weapon  usage 
occurred.  When  BRIK's  DE  goals  were  gradually  increased  until 
all  of  the  weapons  were  used,  a  difference  of  1%  in  the  value 
destroyed  was  attained.  In  an  attempt  to  determine  why  the 
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differences  occurred,  the  single  shot  probabilities  of 
survival  as  computed  by  BRIK  were  compared  with  those  of  AEM. 
As  mentioned  in  the  verification  section,  these  numbers 
compared  favorably.  It  was  determined  that  the  differences  in 
destroyed  value  were  attributed  to  the  damage  function 
approximation  described  in  Chapter  III.  This  approximation 
‘error*  will  be  further  described  in  Chapter  VII. 

BRIK  Capabilities.  BRIK  is  not  a  single  use  model. 
Instead,  due  to  the  various  objective  functions  and 
constraining  rules  that  can  be  combined  for  use  in  any 
problem,  BRIK  has  the  versatility  to  be  used  in  many  types  of 
nuclear  exchange  analyses.  This  section  will  demonstrate  the 
use  and  interaction  of  BRIK's  objective  functions  and  hedges 
in  the  allocation  process. 

Examples  will  be  used  to  clarify  the  use  of  the  various 
rules  that  can  be  employed  in  BRIK.  These  examples  will  be 
small  and  simple,  so  that  the  effect  of  each  function  or 
constraint  is  not  hidden  in  the  at  location.  A  target  base 
<Figure  4-A>  of  two  classes  will  be  attacked  by  a  weapon  base 
{Figure  4-B>  of  two  classes  for  all  of  this  section's 


examples.  The  OE  goals  for  both  target  classes  w 


ill  be  . 
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with  enough  weapons  available  to  meet  these  OE  goals  unless 

i 

specified  otherwise.  The  extreme  goal  used  in  alii  problems 
was  to  minimize  warheads  used.  The  first  examples  will 
demonstrate  the  three  objective  functions.  The  remaining 
examples  will  describe  the  effect  of  each  hedge. 
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NAME 

NUMu 

VUL 

DIA 

TYPE 

PRI0 

0AL 

REL 

CEP 

YLD  UHDS 

civil 

10 

2508 

.49 

m 

1 

1.00 

,00 

,00 

.00  0 

local 

10 

1001 

.49 

n 

1 

1.00 

.00 

.00 

.00  0 

Figure  4-A.  Example  Target  Base 


2 

NAME 

NUMWPNS 

UHD/WPN 

REL  CEP  YLD  DAYALRT 

GENALRT 

titan 

49 

1 

,80  .900  4.00  .85 

1.00 

MMII 

45 

1 

,85  .300  1,00  .85 

1.00 

Figure  4-B.  Example  Weapon  Base 

Objective  Functions:  There  are  three  objective  -functions 
in  BRIK  from  which  to  choose.  They  include  the  following: 

1.  OE  goals  are  met  as  best  possible  using  only  the 
available  arsenal. 

2.  Force  all  DE  goals  to  be  met,  regardless  of  the  weapon 
avai 1 abi 1 i ty. 

3.  Bounds  the  problem  with  both  weapon  availability  and  DE 
goal  levels.  This  function  requires  hedging  constraints 
or  a  number  5  extreme  goal — forcing  BRIK  to  use  all 
remaining  weapons — to  run  an  allocation. 
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The  following  sample  output  <F-gure  5)  is  from  an 


allocation  with  a  type  one  objective  function.  Both  target  DE 
goals  were  met  using  33.4  MMIIs.  The  remaining  MMIIs  and 
Titans  were  not  allocated  and  were  placed  in  reserve. 


TOTAL  VALUE 

DESTROYED  UAS  19.0 

♦ 

TGTNAH 

WPNCLASS 

NUMBER 

ASSIGNED 

GOAL 

NUMBER 

ACHIEVEMENT  REMAINING 

civil 

MMII 

17,45 

.95 

.95  .5 

local 

MMII 

15.79 

.95 

*95  ,5 

WPNAME 

NUMBER 

USED 

REMAIN 

titan 

49.0 

.0 

49.0 

MMII 

45.0 

33.4 

11.4 

Figure  S.  Objective  One  Example  Output 

The  second  objective  function  forces  all  DE  goals  to  be 
met  and,  if  necessary,  creates  weapons  to  meet  the  DE  goals. 
The  function  will  either  build  the  weapon  classes  and  numbers 
needed  from  scratch  or  i t  can  build  additional  weapons  in 
addition  to  those  already  included  in  the  weapon  base.  Only 
weapon  classes  already  in  the  weapon  base  will  be  used  in  the 
allocation.  The  following  example  (Figure  6>  is  run  with  both 
weapon  classes  having  0  weapons.  In  order  to  meet  its  goals, 
BRIK  had  to  create  33.4  MMIIs,  which  is  the  same.’  amount  needed 
in  the  first  objective  example.  The  weapon  usage  output 
indicated  that  -33.4  MMIIs  were  remaining.  This  is 
interpreted  as  adding  33.4  MMIIs  to  the  existing  arsenal. 


If  there  are  already  weapons  in  the  inventory,  BRIK  will 
use  those  weapons  already  existing  before  adding  more  to  the 
weapon  base.  For  example,  for  the  solution  reported  in  Figure 
7,  instead  of  8  weapons  in  both  classes,  the  Titan  class 
initially  had  ten  weapons  and  the  MMII  class  had  8.  The 
allocation  used  all  the  18  Titans  and  added  25  Mills  to  meet 
the  goals. 

The  type  two  objective  function  was  specifically  designed 
to  use  whatever  minimum  weapons  were  made  available  pridr  to 
'building*  more.  If  the  user  provides  a  sufficient  arsenal  to 
meet  the  DC  goals  and  uses  the  type  2  objective  function,  BRIK 
wilt  allocate  the  maximum  number  of  weapons  possible  to  meet 
those  goals. 


TOTAL  VALUE 

DESTROYED 

WAS  19,0 

♦  . 

TGTNAM 

civil 

local 

WPNCLASS 

MMII 

MMII 

NUMBER 

ASSIGNED 

17.65 

15.79 

GOAL 

.95 

.95 

NUMBER 

ACHIEVEMENT  REMAINING 
.95  .5 

.95  .5 

WPNAME 

titan 

MMII 

NUMBER 

.0 

.0 

USED 

.0 

33.4 

REMAIN 

.0 

-33,4 

Figure  6.  Objective  Two  Example  Output  -  8  weapons 
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TOTAL  VALUE  DESTROYED  WAS  19.0. 


NUMBER  NUMBER 


TGTNAM 

civil 

local 

WPNCLASS 
MMI  I 
ti  tan 

MMI  I 

ASSIGNED 

17.65 

10.00 

7.31 

GOAL 

.95 

.95 

ACHIEVEMENT 

.95 

.95 

REMAINING 

.5 

.5 

UPNAME 

titan 

MMI  I 

NUMBER 

10.0 

.0 

USED 

10.0 

25.0 

REMAIN 

.0 

-25.0 

Figure  7.  Objective  Two  Example  Output  -  10  Titans 

The  next  example  is  an  allocation  using  a  type  three 
objective  -function  and  no  hedging  or  type  5  extreme  goal 
constraints.  The  allocation  places  an  upper  limit  on  the  DE 
that  can  be  accomplished  and  the  weapons  that  can  be  used. 

The  type  1  extreme  goal  used  minimizes  the  weapons. 

Therefore,  since  there  is  no  minimum  OE  level  specified  and 
the  extreme  goal  forces  a  minimization  of  weapons,  none  of  the 
weapons  are  used  <Figure  8).  If  a  hodge  is  added  to  force  a 
.5  OE  on  the  target  class  "civil",  BRIK  will  allocate  4.1  MMI I 
weapons  to  achieve  the  DE  level  {Figure  9) .  It  should  be 
noted  that  BRIK  used  4.1  weapons  to  destroy  5  targets;  this  is 
due  to  the  same  damage  function  assumption  which  was  evident 
in  the  AEM  comparison  in  the  previous  section.  This  pr obi em 
is  accentuated  in  small  problems,  and  is  not  as  evident  in 
larger,  more  complex  problems. 


TOTAL  VALUE 

DESTROYED 

WAS  .0 

♦ 

TGTNAH 

civil 

lOCQl 

WPNCLASS 

NUMBER 

ASSIGNED 

GOAL  ACHIEVEMENT 

.95  .00 

.95  .00 

NUMBER 

REMAINING 

10,0 

10.0 

WPNAME 

titan 

MMII 

NUMBER 

49.0 

45.0 

USED 

.0 

.0 

REHAIN 

49.0 

45.0 

Figure  8. 

Objective  Three  Example  Output  -  No 

Hedges 

TOTAL  VALUE 

DESTROYED  WAS  5.0. 

TGTNAH 

civil 

local 

WPNCLASS 

MMII 

NUMBER 

ASSIGNED 

4.08 

GOAL  ACHIEVEMENT 

.95  .50 

'95  .00 

NUMBER 

REMAINING 

5.0 

10.0 

WPNAME 

titan 

MMII 

NUMBER 

49.0 

45.0 

USED 

.0 

4.1 

REMAIN 

49,0 

40,9 

Figure  9.  Objective  Three  Example  Output  -  ('edge 

Another  way  to  -force  a  type  three  objective  -function  to 
allocate,  is  to  include  a  type  5  extreme  goal  to  -force  BRIK  to 
use  all  of  the  remaining  weapons.  To  demonstrate  this  method, 
both  weapon  classes  were  reduced  to  10  weapons  each  and  the 
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results  i.*e  displayed 

BRIK  to  allocate  all  : 

in  Figure 

20  weapons. 

10. 

The  extreme  goal  forced 

TOTAL  VALU 

E  DESTROYED 

WAS  13.7. 

NUMBER 

NUMBER; 

TGTNAM 

URN CLASS 

ASSIGNED 

GOAL 

achievement  REMAINING 

civil 

1 1  tran 

10.00  ■ 

.95 

.37  1.3 

MMII 

6.35 

local 

MMII 

3-65 

.  95 

>50  5.0 

WPNANE 

NUMBER 

USED  ? 

EM  A  IN 

titan 

10.0 

10 . 0 

.0 

MMII 

10.0 

10.0 

A 

Figure  10.  Objective  Three  Example  Output  -  Type  5  Extreme 
Coal 

Hedges:  The  -following  seven  constraining  hedges  are 
available  in  BRIK: 

1.  En-force  a  minimum  level  o-f  damage  on  a  particular  target 
class. 

2.  En-force  a  minimum  level  o-f  damage  on  a  particular  target 
class  with  a  speci-fic  subset  o-f  weapon  classes. 

3.  Allow  only  a  certain  amount  o-f  damage  <50  on  a  particular 
target  class  resulting  -from  a  specific  set  of  weapons. 

4.  Restrict  the  number  of  weapons  from  a  weapon  class  that  can 
be  allocated  to  a  specific  target  class. 

5.  Allows  the  construction  of  a  constraint  that  is  not 
indigenous  to  B.1IK. 
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6.  Cn-force  a  minimum  level  of  damage  on  a  particular 
aggregated  set  of  target  classes. 

7.  Restrict  the  amount  of  total  weapons  that  can  be  allocated 
to  each  target  in  a  specific  target  class. 

If  there  is  a  certain  minimum  level  of  damage  that  must  be 
attained  on  a  target  class,  a  type  one  nedging  constraint  can 
be  used.  This  hedge  was  used  in  Figure  8  to  force  the  third 
objective  function  to  allocate.  If  it  was  desired  that  the  DE 
level  be  accomplished  using  Titans  instead  of  MMIIs,  a  type 
two  hedge  could  be  used.  This  would  allow  the  user  to 
specifically  designate  the  titan  class  as  the  only  weapon 
class  to  be  allocated  to  target  class  "civil*  to  attain  the 
specified  OE  level  (Figure  11). 


TOTAL  VALUE 

DESTROYED 

WAS  5.0. 

• 

TGTNAM 

civil 

local 

V 

WPNCLASS 

titan 

NUMBER 

ASSIGNED 

7.27 

GOAL 

.95 

,95 

NUMBER 

ACHIEVEMENT  REMAINING 
,50  5.0 

.00  10.0 

WPNAME 

titan 

MMII 

NUMBER 

49.0 

45,0 

! 

USED 

7.3 

.0 

REMAIN 

41.7 

45.0 

Figure  11.  Example  of  Type  Two  Hedging 


To  demonstrate  a  type  3  hedge*  consider  the  -following. 
Figure  4  showed  output  when  an  objective  function  one  was  used 
to  attain  .95  DE  levels.  If  it  is  desired  to  restrict  the 
Mills  from  inflicting  greater  than  50/C  of  the  DE  on  the 
"civil"  class*  a  type  three  hedging  constraint  can  be  used. 
This  hedge  allowed  BRIK  to  <use  only  4.08  Mills  to  attain  a  .5 
DE  of  the  target  class  (Figure  12).  Titans  are  used  to  meet 
the  rest  of  the  D’£  goal  . 


TOTAL  VALUE 

DESTROYED 

WAS  19.0 

♦ 

TGTNAM 

NUMBER 

NUMBER 

WPNCLASS 

ASSIGNED 

GOAL 

ACHIEVEMENT 

REMAINING 

civil 

titan 

24.16 

.95 

.95 

»5 

MMII 

4.08 

local 

MMII 

15.79 

.95 

.95 

.5 

WPNAME 

HUMBER 

USED 

REMAIN 

titan 

47.0 

24.2 

24.8 

MMII 

45.0 

19.9 

25.1 

Figure  12.  Example  of  Type  Three  Hedgir-.q  -  Single  Weapon 
Class 

If  the  weapon  subset  included  both  "titan”  and  "Mill”  classes, 
BRIK  will  only  attain  a  DE  level  of  .5  on  the  "civil"  class 
(Figure  13) . 
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TOTAL  VALUE  DESTROYED  WAS 


14.5* 


NUMBER 


TGTNAM 

WPNCLASS 

ASSIGNED 

civil 

MMII 

4.08 

local 

MMII 

15.79 

WPNAME 

NUMBER 

USED 

titan 

49.0 

.0 

MMII 

45.0 

19.9 

NUMBER 

GOAL 

ACHIEVEMENT  REMAINING 

.95 

a  50  5*0 

.95 

*95  *5 

REMAIN 

49.0 

25.1 

Figure  13.  Type  Three  Hedging  -  Both  Weapon  Classes 


If  it  is  desired  to  limit  the  number  of  a  weapon  class  to 
be  used  against  a  specific  target  class,  a  type  four  hedging 
constraint  can  be  added.  Figure  14  lists  the  output  for  a 
type  one  objective  function.  A  type  four  hedge  was  used  to 
prevent  more  than  5  MMIIs  from  being  used  in  attaining  the 
.93  DE  goal  against  target  class  "local". 


TOTAL  VALUE 

DESTROYED 

WAS  19.0 

♦ 

NUMBER 

NUMBER 

TGTNAM 

WPNCLASS 

ASSIGNED 

GOAL 

ACHIEVEMENT  REMAINING 

civil 

MMII 

17.65 

,95 

►  95  .5 

local 

titan 

12,72 

,95 

.95  ,5 

MMII 

5,00 

WPNAME 

NUMBER 

USED 

REMAIN 

titan 

49.0 

12,7 

36.3 

MMII 

45.0 

22.6 

22,4 

Figure  14.  Example  of  a  Type  Four  Hedge 
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Since  BRIK's  set  of  hedges  is  not  inclusive  of  all  the 
possible  allocation  ru’es  that  might  be  needed  for  a  problem, 
the  ability  to  input  a  user-built  constraint  is  included  in 
BRIK.  This  is  a  type  5  hedging  constraint.  One  possible  use 
of  this  hedge  is  to  account  for  a  common  carrier.  If  a 
weapon  system  contains  two  different  types  of  warheads  for 
use  in  a  specific  ratio,  such  as  a  B-52  carrying  gravity 
bombs  and  SRAMs,  this  function  can  be  used  to  build  a 
constraint  that  only  allows  use  of  weapon  systems  in  a 
specific  ratio.  For  example,  if  it  is  desired  to  use  the 
same  exact  amount  of  titans  as  Mills,  a  hedge  can  be  built  to 
force  the  same  use  of  both  weapon  classes.  The  output  from 
an  allocation  which  used  a  type  one  objective  function  and  a 
user-built  hedge  is  shown  in  Figure  15. 

TOTAL  VALUE  DESTROYED  WAS  19.0.  | 


TGTNAH 

civil 

local 

WPNCLASS 

MMII 

titan 

NUMBER 

ASSIGNED 

17.65 

18.10 

GOAL 

.95 

.95 

NUMBER 

ACHIEVEMENT  REMAINING 
.95  ,5 

.95  .5 

WPNAME 

titan 

HMII 

NUMBER 

49.0 

45.0 

USED 

18.1 

18.1 

REMAIN 

30.9 

26.9 

Figure  15.  Example  of  a  User-Built  Hedge  -  Type  5 
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In  the  screen  output,  BRIK  does  not  report  any  individual 
weapon  allocation  less  than  .5  ,  but  it  does  add  up  all  o-f 
the  weapons  allocated  to  determine  the  usage.  This  describes 
why  there  is  a  dexcrepancy  between  the  number  of  weapons 
allocated  and  the  number  of  weapons  that  are  reported  used. 

An  inspection  of  "alloc*  indicated  that  .45  MMIIs  were 
allocated  against  the  target  class  "local". 

Hedge  6  can  be  used  to  create  a  super  target  class  -from 
several  smaller  ones.  There  might  be  times  when  it  is 
desired  to  attain  a  DE  level  against  an  entire  category  of 
classes  and  not  be  particularly  interested  in  the  DE  on 
speci-fic  target  classes,  especially  if  the  target  class 
characteristics  vary  too  much  to  allow  aggregation  into  a 
single  class.  For  this,  a  type  six  hedge  is  used.  It  is 
recommended  that  this  hedge  be  used  only  with  the  type  3 
objective  function,  because  this  function  places  an  upper 
bound  on  the  DE  that  can  be  attained  on  each  class.  Without 
this  upper  bound,  BRIK  could  allocate  all  of  its  weapons  to 
the  softest  target  in  the  set.  The  results  in  Figure  16  are 
from  an  example  that  used  a  type  three  objective  function  and 
a  type  six  hedge  of  .5  DE  on  a  combination  of  both  classesT 
There  are  10  targets  in  each  class  for  a  total  of  20  targets. 
However,  BRIK  allocated  7.31  MMIIs  to  destroy  7.5  "local" 
targets,  which  was  below  the  desired  10  target  destruction 
goal.  This  is  due  to  the  damage  function  approximation  and 
will  be  discussed  in  Chapter  VII. 


97 
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DESTROYED 

WAS  7.5 

♦ 

TGTNAH 
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GOAL 

ACHIEVEMENT  REMAINING 
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.95 

.00  10.0 
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.75  2.5 
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titan 
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.0 

49.0 
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43.0 

7.3 

37.7 

Figure  16.  Example  of  a  Type  Six  Hedge 

A  type  seven  hedge  restricts  the  amount  of  weapons  that 
can  be  used  against  each  target  in  a  target  class.  Figure  17 
shows  example  output  from  an  allocation  that  used  a  type  one 
objective  function  and  two  type  7  constraints.  The 
constraints  forced  a  maximum  of  one  weapon  per  target  in  each 
target  class.  There  were  ten  targets  in  each  class, 
therefore,  only  ten  weapons  were  a! 1 ocated  to  each  class. 


TOTAL  VALUE  DESTROYED  WAS  16.7 ♦ 
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Figure  17. 

Example  of  a  Type  Seven  Hedge 
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Thia  concludes  the  demonstration.  BRIK  incorporates 
several  objective  functions  and  many  allocation  rules  into  an 
extremely  versatile  nuclear  exchange  model.  The  next  chapter 
will  conclude  the  report  with  some  observations  concerning 
BRIK's  behavior  and  will  suggest  areas  'or  further  study. 
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VI I  Concl usi on 


BRIK  is  unlike  any  other  nuclear  exchange  model  because 
it  uses  preemptive  goal  programming  <G.P.)  to  meet  all  of  the 
allocation  rules  and  to  assign  weapons  to  targets.  Because 
of  G.P.,  BRIK  is  extremely  versatile.  This  final  chapter  is 
divided  into  three  sections  which  discuss  this  versatility. 
First,  a  summary  of  the  modeling  process  will  be  presented  in 
order  to  provide  a  perspective  of  the  underlying  process 
which  led  to  the  development  of  BRIK.  The  second  section 
will  provide  a  discussion  of  some  potential  applications  of 
BRIK.  Finally,  the  third  section  will  present  observations 
concerning  some  of  BRIK's  quirks. 

Summary 

This  project  started  as  a  continuation  of  a  thesis  defen¬ 
ded  in  December  1982  <Ref  27) .  Michael  C.  Mambsganss  demon-  ■ 
strated  that  auxiliary  variables  in  linear  programming  con¬ 
straints  could  be  used  to  allocate  weapons  to  targets  in 
order  to  meet  user-defined  damage  expectancy  goals. 

Wambsganss  employed  a  damage  function  approximation  whtre  his 
model  allocated  non-integer  weapons  to  each  target  in  a 
target  class.  Ne*t,  to  increase  solvability,  he  linearized 
the  damage  function  with  logarithms,  permitting  the  use  of  LP. 
Finally,  he  included  auxiliary  variables  in  each  target  con¬ 
straint.  With  these  auxiliary  variables  appearing  in  a  sin¬ 
gle  objective  function,  Hambsganss  demonstrated  that  an  allo¬ 
cation  could  be  made  to  meet  user-defined  DE  goals.  From 
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this  preliminary  work,  BRIK's  evolution  into  an  elaborate, 
user-friendl y  model  is  a  noteworthy  advancement. 

Flexibility  was  the  -first  requirement  placed  on  BRIK.  A 
model  is  -flexible  if  it  can  handle  many  different  types  of 
allocations.  In  a  nuclear  exchange  model,  this  is  best 
handled  by  permitting  numerous  allocation  rules.  Various 
models  were  reviewed  to  determine  what  rules  and  options  are 
currently  in  use.  From  the  list  compiled  by  this  review, 
BRIK's  hedging  options  were  developed.  The  model  also  in¬ 
cludes  the  ability  to  prevent  a  designated  weapon  class  from 
attacking  any  specific  target  class.  These  options  give  BRIK 
the  flexibility  that  was  originally  desired. 

The  next  requirement  placed  on  BRIK  was  the  ability  to 
compute  single  shot  probability  of  survival  <SSPS> . 
Mambsganss'  model,  computed  SSPS  given  the  vul  nertbi  1  i  ty  of  a 
target  class  expressed  as  peak  overpressure — a  PSI  number. 
Because  this  system  is  inaccurate,  the  authors  decided  that 
the  usefulness  of  BRIK  would  be  increased  if  the  model  in¬ 
cluded  the  nuclear  targeting  vu 1 nerabi 1 i ty  <VN)  system  used 
by  the  Defense  Intelligence  Agency.  This  was  accomplished  by 
using  Subroutine  DPDX  from  the  SAC  BLl»J  SWADE  model  .  Because 
VN  numbers  may  not  always  be  available,  BRIK  also  retained 
the  ability  to  compute  SSPS  from  vulnerability  expressed  as 
PSI  numbers. 

It  was  also  felt  that  BRIK  should  be  user-f r i endl y . 

Large,  complicated  models  sometimes  require  many  weeks  of 
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study  before  they  can  be  adequately  used.  To  reduce  this 
learning  time  -for  BRIK,  the  program  was  written  to  supply 
prompts  which  guide  a  user  through  each  section  of  a  problem. 
BRIK  then  uses  the  provided  information  to  perform  all  of  the 
mathematics  internally.  For  example,  to  build  a  hedge  which 
enforces  a  minimum  amount  of  damage  on  a  designated  target 
class  resulting  from  a  particular  set  of  weapon  classes,  a 
user  only  has  to  perform  a  few  simple  steps.  First,  the  user 
is  asked  to  choose  the  desired  type  of  hedge  from  a  menu. 
Next,  he/she  is  prompted  to  supply  the  number  of  the  desig¬ 
nated  target  class  and  the  minimum  DE.  Finally,  the  user  is 
asked  to  input  the  weapon  class  numbers  that  should  be  used 
to  attain  that  DE.  BRIK  now  uses  all  of  this  information  to 
internally  build  the  necessary  constraint. 

Another  requirement  imposed  on  BRIK  was  transpor tabi 1 i ty . 
Since  use  of  a  local  library  package  might  inhibit  the  use¬ 
fulness  of  thiv.  model,  the  literature  once  again  was  searched 
for  an  appropriate  allocation  routine.  At  first,  the  search 
was  for  linear  programming  packages,  but,  as  the  design  of 
BRIK  progressed,  goal  programming  routines  were  also  sought. 
PAGP,  the  Partitioning  Algorithm  for  Goal  Programming,  was 
chosen  because  it  offered  several  efficiencies  of  which  BRIK 
could  take  advantage.  These  included  partitioning  of  the 
constraints,  an  indicator  array  which  reduced  the  number  of 
variables  that  needed  to  be  considered  for  entry  into  the 
basis,  and  special  stopping  rules.  The  routine  is  copy- 
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righted,  therefore  it  cannot  be  used  -for  commercial  qain. 
However,  the  Association  -for  Computing  Machinery  gave  permis¬ 
sion  to  include  it  in  this  model. 

The  use  of  G.P.  permitted  BRIK  to  evolve  into  a  multipur¬ 
pose  model.  BRIK  has  the  capability  to  build  different 
objective  functions  allowing  vastly  different  questions  to  be 
asked  about  the  same  LP  matrix.  This  capability  will  be 
further  discussed  in  the  next  section. 

ii§ £ 

Models  are  nothing  but  "black  boxes"  <Ref  24:53).  One 
should  always  remember  that  BRIK  is  only  a  model  and  behaves 
according  to  the  algorithms  written  into  it.  Given  the  same 

input,  it  will  give  the  same  output;  but  two  different 

! 

analysts  may  come  up  with  different  answers  to  the  same 

i 

question.  This  is  due  to  the  fact  that  allocations  are  made 

within  a  set  of  constraints  which  should  reflect  strategy  and 

! 

operational  considerations.  These  constraints  can  be  input  as 

i 

I 

allocation  rules.  Because  of  BRIK's  flexibility,  a  virtually 

i 

endless  variety  of  allocation  rules  are  available.  Also, 
because  of  the  different  objective  functions,  various  kinds 
of  analyses  can  be  performed. 

BRIK  can  be  used  to  analyze  build-down  strategies.  In 
the  spring  of  1983,  US  Senator  William  Cohen  <D-Maine>  sug¬ 
gested  a  strategy  to  stabilize  the  current  nuclear  forces. 

He  presented  the  idea  of  removing  a  number  of  weapons  from 
the  existing  arsenal  in  exchange  for  new  weapons  added  to  the 
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force.  BRIK  is  a  tool  that  can  aid  analysis  of  various 
build-down  strategies.  The  following  is  a  description  of  how 
the  model  could  be  used  for  this  purpose.  It  is  first  neces¬ 
sary  to  determine  the  DE  capability  of  the  existing  nuclear 
arsenal.  Next,  a  problem  can  be  run  using  the  current 
forces,  a  current  target  base,  and  the  DE  goals  that  the 
force  can  attain.  Then,  the  sensitivity  section  of  BRIK  is 
used  to  increase  weapon  availability  on  some  select  system. 

If  the  first  objective  function — meet  goals  with  some  minimum 
aspect  of  the  available  arsenal  —  is  used,  some  weapons  will 
go  unused  when  the  problem  is  rerun.  If  the  unused  weapons 
are  different  than  the  weapons  added,  the  unused  weapons 
could  become  candidates  for  build-down. 

BRIK  could  also  be  used  to  determine  future  force  struc¬ 
ture.  The  second  objective  function  can  be  used  to  "build" 
the  number  of  weapons  needed  to  meet  user-defined  DE  goals. 

To  choose  between  alternative  systems,  the  following  steps 
could  be  performed.  A  weapon  base  that  includes  the  alter¬ 
native  weapon  classes  should  be  entered  along  with  a  target 
base  and  appropriate  DE  goals.  After  inputting  all  of  the 
appropriate  allocation  rules,  the  problem  would  be  run.  The 
output  would  show  which  systems  were  "built,"  thus  giving  the 
analyst  indication  as  to  which  system  to  purchase. 

A  third  type  of  analysis  that  BRIK  is  well  suited  to 
perform  is  a  sensitivity  analysis.  Because  of  the  interac¬ 
tive  qualities  of  this  model,  an  analyst  can  sit  at  a  termi- 
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swers  to  these  questions  are  displayed  on  the  screen  enabling 
the  analyst  to  immediately  see  the  results  of  the  changes 
made  to  the  problem.  This  sensitivity  capability  is  a  step 
toward  building  the  tools  of  strategic  force  analysis  around 
the  machines  used,  making  the  machines  extensions  of  human 
thought. 

BRIK  has  a  tremendous  amount  of  potential  to  aid  analy¬ 
sis  in  many  different  areas  of  nuclear  exchange  as  well  as  a 
versatility  heretofore  unknown.  However,  it  is  not  perfect. 
The  following  section  describes  two  areas  a  user  may  find  as 
shortfalls  to  the  model. 

Observations 

There  are  two  aspects  of  BRIK's  behavior  which  may  concern 
potential  users  of  this  model.  The  first  area  deals  with 
BRIK's  current  lack  of  capability  to  maximize  damage  expect¬ 
ancy  while  the  second  area  concerns  the  damage  function 
approximation.  After  presenting  observations  depicting 
BRIK's  behavior  in  these  two  areas,  a  discussion  of  possible 
corrections  follow. 

The  first  area  of  concern  is  BRIK's  lack  of  capability  to 
maximize  DC.  Given  a  weapon  base  and  a  target  base,  most 
nuclear  exchange  models  will  maximize  damage  expectancy 
<Ref  24:53).  To  see  if  BRIK  has  this  capability,  the  test 
problem  used  to  compare  BRIK  with  the  Arsenal  Exchange  Model 
<AEM)  was  run  with  DE  goals  of  .98  for  each  target  class. 

The  results  are  in  Table  X.  Comparing  this  table  with  the 
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solution  -from  AEM  in  Table  VII,  it  can  be  seen  that  BRIK's 
allocation  is  not  close.  Also,  the  total  value  destroyed  by 
BRIK,  1367.2,  is  considerably  lower  than  the  damage  achieve¬ 
ment  by  AEM  of  1770.  Therefore,  in  its  current  form,  BRIK 
does  not  maximize  damage  expectancy. 

The  second  area  of  concern  is  the  "optimistic"  DE  from 
the  damage  function  approximation  described  in  Chapter  III. 
Figure  18  shows  the  results  of  a  small  problem  where  an 
achievement  of  .5  was  attained  on  target  class  "civil."  From 
vthe  figure  it  can  be  seen  that  4.08  weapons  were  assigned  to 
the  class.  Since  the  class  originally  had  10  targets  and 
there  are  S  targets  remaining,  5  targets  were  destroyed  by 
4.08  weapons!  These  areas  are  mentioned  because  they  may  be 
of  concern  to  potential  users  of  BRIK  in  its  current  form. 
However,  this  does  not  imply  that  both  limitations  need  be 
permanent. 

It  is  believed  that  BRIK's  current  lack  of  capability  to 
maximize  DE  can  be  solved.  One  potential  technique  is  an 
iterative  scheme  that  systematical  1 y  changes  DE  goals  on  the 
target  constraints.  The  iterations  should  concentrate  on 
allocating  weapons  to  targets  that  have  not  been  assigned  any 
weapons  before  allocating  weapons  to  targets  that  already 
have  weapons  asigned.  This  could  possibly  solve  the  first 
area  of  concern.  Giving  BRIK  the  capability  to  maximize  DE 
seems  realistic  and  would  be  a  good  addition  to  the  model. 


Teat  Problem  BRIK  Solution t  High  DE  Goals 
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T3TAL  VALUE 

DESTROYED 

WAS  5.0. 

TGTMAM 

civil 

local 

WPNCLASS 

MMII 

NUMBER 

ASSIGNED 

4.08 

GOAL  • 
.95 
.95 

ACHIEVEMENT 

.50 

.00 

NUMBER 

REMAINING 

5.0 

10.0 

WPNAME 

titan 

MMII 

NUMBER 

49.0 

45.0 

USED 

.0 

4.1 

REMAIN 

49.0 

40.9 

Figure  18.  Sample  Problem 

I 

However,  this  deficiency  does  not  negate  the  advantages 
gained  from  the  availability  of  this  tool.  It  matches  the 

i 

methodology  of  sequential  achievement  of  damage  goals,  as 

outlined  in  Reference  24. 

| 

The  damage  function  approximation  in  BRIK  results  from 

I 

assigning  non-integer!  weapons  to  each  target  in  a  class.  It 

I 

is  used  internally  i n ;  two  different  places.  First,  it  is 
used  to  build  all  of  the  target  constraints  and  the  hedging 
constraints  which  deal  with  minimum  or  intermediate  DE  goals. 
Second,  the  damage  function  approx imati on  is  used  in  the 
output  section  to  compute  the  achievement  on  each  target 
class.  To  eliminate  the  damage  function  approximation,  the 
model  could  be  completely  reformulated  where  the  decision 
variable  is  some  number  of  targets  that  receive  integer 
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weapons.  From  the  authors"  points  o f  view,  this  would  be  an 
extremely  difficult  task.  It  is  possible,  though,  to  take  the 
allocation  and  recompute  the  damage  achievement.  For  exam¬ 
ple,  in  Figure  18,  since  4.08  weapons  were  assigned  to  target 
class  "civil,”  the  achievement  could  be  computed  by  assuming 
4.08  targets  were  each  allocated  one  weapon  and  5.92  targets 
were  allocated  no  weapons.  Recomputing  the  DE  in  this  manner 
would  allow  for  a  more  accurate  representation  of  the  achieve¬ 
ment,  however,  it  would  also  result  in  some  inconsistency  in 
the  model.  Again,  using  Figure  18  as  an  example,  4.08  MMII 
warheads  were  assigned  to  target  class  "civil*  by  using  a 
hedge  that  asked  for  a  DE  goal  of  .5.  The  current  output 
reflects  that  this  goal  was  met  and  there  were  remaining 
weapons.  If  the  achievement  was  recomputed  as  described 
above  where  4.08  targets  received  1  weapon  and  5.92  targets 
received  no  weapons,  the  achievement  would  have  been  reported 
as  .33,  compared  with  the  desired  goal  of  .5.  By  making 
conparisons  between  BRIK's  non-integer  allocations  with  inte¬ 
ger  solutions,  this  difference  in  achievement  is  the  largest 
one  that  was  found.  While  the  current  damage  function 
approximation  makes  weapons  appear  optimistic,  the  difference 
is  consistent  within  the  model  as  it  is  currently  formulated. 
Also,  reported  DE  is  optimum  if  fractions  of  weapons  could  in 
fact  be  allocated.  Therefore,  BRIK's  use  as  an  analysis  tool 
is  not  compromised  by  the  presence  of  this  damage  function 
approximation . 
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BRIK  started  with  an  LP  -formulation  and  a  single  objec¬ 
tive  -function.  It  evolved  into  a  model  that  contained  a 
unique  application  o-f  G.P.,  -from  which  a  tremendous  amount  o-f 
capability  is  enjoyed.  The  deviational  variables  and  how 
they  are  handled  in  the  different  objective  functions  is  the 
key  to  BRIK's  flexibility  and  potential  for  multipurpose 
application.  BRIK  permits  the  user  to  determine  its  applica¬ 
tion  and  to  examine  numerous  objectives  and  strategies 
providing  a  virtually  acustom  designed"  model  tailored  to  the 
individual  user's  requirements. 
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APPENDIX  A 


Ust  Quid* 

This  guide  will  -familiarize  the  user  with  many  o-f  the 
-features  o-f  BRIK.  It  will  step  through  BRIK  in  the  order  o-f 
appearance  o-f  the  various  model  -functions.  Examples  will  be 
supplied  when  necessary  for  clarity. 

BRIK  was  developed  to  be  "user  -friendly".  The  standard  o-f 
"user  -friendly"  is  taken  to  mean  that  anyone  with  the 
requisite  data  and  limited  background  in  either  computer  or 
mathematical  programming  could  use  the  model  -for  an 
allocation.  Since  BRIK  displays  prompts  and  menus  on  the 
screen  to  guide  the  user  through  the  model's  -functions,  the 
user  needs  access  to  a  computer  equipped  with  a  monitor. 

Many  o-f  the  prompts  which  appear  on  the  screen  are 
questions  that  require  answers.  Realizing  that  a  user  could 
input  an  illogical  response  or  an  incorrect  piece  o-f  data,  the 
model  incorporates  many  internal  checks  to  guard  against  the 
possibility  of  an  unintentional  mistake.  An  example  of  a 
check  can  be  found  at  any  yes  or  no  <y/n>  question.  If  the 
user  gives  an  illogical  response,  the  cursor  drops  a  line  and 
waits  for  a  correct  response,  y  or  n.  Also,  the  model 
discourages  the  user  from  inputting  inappropr iate  data.  For 
example,  probabilities  are  accepted  only  if  they  are  between  0 
and  1.  The  procedure  is  the  same  as  the  yes/no  question;  the 
cursor  will  wait  for  the  proper  input. 

On  some  menus,  the  user  must  enter  a  0  to  terminate  that 
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section  o-f  the  model.  It  has  been  discovered  that  the  VAX 
computer  interprets  any  non-numeric  character  as  a  0  when  it 
is  expecting  an  integer.  It  is  unknown  if  any  other  computer 
has  the  same  logic.  Therefore,  the  user  should  be  cautioned 
that  if  any  character  value  is  entered  for  an  integer,  that 
particular  model  section  can  terminate  prematurely. 

BRIK  can  be  used  for  two  types  of  analysis.  The  first 
type  is  used  for  any  new  problem  where  the  user  must  input 
data,  select  the  type  of  objective  function,  and  enter  the 
allocation  rules  for  the  problem.  The  second  type  enables  the 
user  to  conduct  a.  sensitivity  analysis  on  some  of  the  input 
parameters.  The  sensitivity  section  assumes  that  an  initial 
type  1  analysis  has  been  accomplished. 

This  user  guide  will  be  broken  up  into  two  main  sections. 
The  first  section  will  discuss  the  interaction  required  to 
conduct  a  new  analysis.  The  second  section  will  describe  the 
capability  of  the  sensitivity  section. 

NEW  ANALYSIS 

This  section  will  describe  the  four  major  sub-areas  of  a 
new  analysis!  infut  requirements,  selection  of  an  objective 
function,  allocation  constraining  rules,  and  BRIK's  output 
data. 

Input  Requirements.  BRIK  allows  the  user  to  input  data 
interactively  or  through  user-designated  files.  The  data 
input  is  divided  into  three  subsections!  target  data,  weapon 
data,  and  damage  expectancy  <DE>  goals.  This  section  will 
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describe  the  type  of  data  needed  by  BRIK  and  the  types  of 
operations  and  editing  that  are  available. 

Oata  is  entered  in  the  -following  order:  target  base, 
weapon  base,  and  DE  data.  Since  BRIK  -follows  the  same  basic 
operation  o-f  input  -for  all  three  data  types,  the  -following 
explanation  applies  to  all  three  types.  The  -first  decision 
the  user  must  make  is  whether  the  data  is  to  be  entered 
interactively  or  -from  a  file.  If  there  has  not  been  any 
previous  analysis  which  created  a  usable  data  base  stored  in  a 
user -design a ted  file,  the  data  must  be  entered  interactively. 
If  the  data  is  entered  interactively,  BRIK  will  guide  the  user 
with  a  series  of  questions.  These  questions  insure  that  the 
required  data  is  input  in  the  correct  fashion. 

Some  input  data  is  entered  as  a  group  and  there  is  no 
individual  check  for  data  correctness.  Therefore,  at  the  end 
of  all  the  data  input  sections,  any  illogical  or  incorrect 
data  will  be  converted  into  a  form  that  BRIK  can  use.  These 
default  values,  except  for  probabilities  and  percentages,  will 
be  indicated  in  the  input  data  lists.  Probabilities  and 
percentages  are  not  allowed  to  be  greater  than  1  or  less  than 
0.  If  data  is  input  greater  than  1,  BRIK  will  decrease  the 
value  to  1;  and  if  a  number  less  than  0  is  input,  it  is 
corrected  to  0.  It  is  also  possible  to  enter  some  data  that 
is  physically  too  large  to  display  on  the  screen  due  to  format 
limitations.  Even  though  the  data  will  be  read  correctly,  all 
screen  output  for  this  data  will  be  a  series  of  stars.  These 
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format  limitations  will  be  listed  in  the  data  lists. 

It  is  possible  to  build  a  data  file  independently  of  BRIK, 
but  it  is  easier  to  input  the  data  interactively  and  allow  the 
model  to  build  an  external  file.  If  the  data  base  is  input 
from  an  external  file,  the  user  must  insure  the  file  exists 
prior  to  operating  BRIK.  If  a  data  file  is  built,  its  name 
must  have  less  than  seven  characters.  Since  the  file  must  be 
formatted,  it  is  easier  to  allow  BRIK  to  create  the  file. 

The  rest  of  this  section  describes  the  input  requirements 
for  the  three  types  of  data  bases  used  in  BRIK.  Each  section 
will  list  and  briefly  describe  the  required  parameters  and 
discuss  the  types  of  editing  that  can  be  done  with  BRIK. 

Target  Input:  The  following  is  a  list  of  required  target 
input  data: 

Class  name  —  Maximum  of  6  characters 

Number  in  the  class  —  Integer,  format  limit  «  9999 

Target  vulnerability  —  VNTK  number  or  peak  overpressure 

needed  to  achieve  a  desired  level 
of  damage.  If  the  VNTK  number  is 
incorrectly  input,  the  program 
will  terminate  after  the  DE  goals 
are  input. 

Target  diameter  —  The  diameter  of  a  circle  in 

nautical  miles  <NM)  that 
encompasses  95 V.  of  the  target. 

If  the  target  is  a  point  target 
the  diameter  input  is  0.  The 
default  is  0  if  a  negative 
number  is  input. 

Target  category  —  FORCE  =  f 

VALUE  =  v,  default  value 
MILITARY  *  m 
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Target  class  priority  —  Integer  1  to  7,  de-fault  is  1  if  a 

value  less  than  1  is  input  and  7  if 
a  value  greater  than  7  is  input. 

Value  —  The  worth  of  a  target  to  the 

attacker.  The  default  is  1. 

The  program  also  has  an  algorithm  to  calculate  values  for 
force  type  targets.  If  the  user  exercises  this  option,  the 
following  additional  target  data  is  needed: 

Target  reliability  —  Probabi 1 i ty  that  the  weapon  located 

at  force  target  will  be  successful 
against  the  attacker. 

CEP  —  Nautical  Miles  ,  defaults  to 

.80031 

Yield  per  warhead  —  Megaton  Yield,  format  limit  is 

999.99 

Warheads  per  target  .  —  The  number  of  warheads  located  at 

each  target.  If  the  target  is  a 
submarine  base,  a  possibility  of 
many  warheads  exists.  Format  limit 
is  99 

If  the  force  value  algorithm  has  been  selected,  the  user  must 
enter  values  for  any  value  or  military  target  classes  so  that 
the  relative  values  between  the  classes  is  represented. 

Target  Editor:  The  first  item  the  user  sees  after 
inputting  the  target  data  is  a  menu  for  editing  the  data 
CFigure  A-l) .  This  menu  allows  the  user  to  either  view  or 
change  any  target  input  parameter  and  to  add,  replace,  or 
delete  an  entire  target  class.  Again,  the  editor  sections  ask 
sel f-expl anatory  questions.  It  should  be  noted  that  after 
making  any  change,  the  user  should  view  the  data  before 
continuing  in  the  program.  This  will  insure  the  new  data  set 
is  correct. 
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If  data  is  to  be  input  via  an  external  -file,  the  user  can 
manually  edit  the  data  file  before  running  BRIK  (Figure  A-2) . 
The  user  should  be  cautioned  not  to  displace  any  data  by  even 
one  space,  as  the  format  BRIK  uses  to  read  the  data  is  an 
exact  one.  Also,  if  the  user  adds  or  deletes  a  class,  the 
first  number  in  the  file  must  be  changed.  This  number 
designates  how  many  classes  are  in  the  file. 


MENU 

Target  Class  Category 

1) 

Check/View  Data 

2) 

Change  Parameters 

3) 

Add  Class 

4) 

Replace  Class 

5) 

Delete  Class 

Enter  your  choice.  Type  0  if  you  do  not  wont  any  of 

these  options 

- 1 - 

or  are  finished  with  your  changes. 

FIGURE  A- I .  Target  Editor  Menu 
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NAME 

NUMB 

VUL 

BIA 

TYPE 

PRI0 

VAL 

REL 

CEP 

YLD 

MHOS 

ldrshp 

200 

20a  5 

,20 

V 

1 

1,00 

.00 

.00 

.00 

0 

iebra 

500 

45oO 

.00 

f 

2 

1.00 

.00 

.00 

.00 

0 

Elbm 

30 

30q0 

.30 

f 

■3 

i.OO 

.00 

,00 

.00 

0 

pal 

500 

30 

1.00 

V 

3 

1,00 

.00 

.00 

.00 

0 

afbase 

50 

35 

1.00 

m 

3 

1.00 

.00 

.00 

.00 

0 

FIGURE  A-2.  Example  of  Target  Base  External  File 

Weapon  Input i  The  following  list  contains  required  weapon 

input  parameters.  As  in  the  target  input  section,  data  can  be 

input  interactively  or  by  a  user -design a ted  file. 

Class  name  — ■  Maximum  of  6  characters 

Number  in  the  class  —  Integer,  format  limit  is  99999 

Warheads  per  weapon  —  The  amount  of  warheads  per  weapon, 

carrier,  format  limit  is  999 

— ■  Probabi 1 i ty  that  the  weapon  will  be 
successful . 

—  Nautical  Miles  <NM> ,  defaults  to 

.88881 

—  Megaton  Yield,  format  limit  is 
999.99 

—  Percentage  of  the  weapon  class 
useable  for  the  allocation 

Generated  alert  rate  —  Percentage  of  the  weapon  class 

useable  for  the  allocation 


Reliability 

CEP 

Weapon  yield 
Daily  alert  rate 
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The  daily  and  generated  alert  rate  data  can  be  used  in  a 
counterstrike  scenario  where  the  daily  alert  rate  can 
.  represent  the  pessimistic  percentage  of  surviving  weapons  and 
the  generated  alert  rate  can  represent  the  optimistic 
percentage  of  surviving  weapons. 

Weapon  Editor:  The  weapon  editor  is  similar  to  the  target 
editor  in  the  types  o-f  available  options  (Figure  A-3>  .  It 
also  has  sel f-expl anatory  questions  to  help  the  user  through 
the  section.  The  weapon  data  -file  is  similar  in  construction 
to  the  target  data  file  (Figure  A-4)  .  It  can  be  manually 
edited,  but  again  caution  is  advised. 


MENU 


Weapon  Class  Category 

1)  Check /View  Data 

2)  Change  Paraneters 

3)  Add  Class 

4)  Replace  Class 

5)  Delete  Class 


Enter  your  choice.  Type  0  if  you  do  not  want  any  of 
these  options  or  are  finished  with  your  changes. 


FIGURE  A-3.  Weapon  Class  Editor 


NAME 

NUMWPNS 

UHD/UFN 

REL 

icbm 

1000 

3 

oer 

♦  /  tJ 

slbi* 

20 

60 

.30 

bomber 

50 

9 

.60 

CEP 

YLD 

DAYALP.T 

0ENALRT 

.130 

.30 

.90 

1.00 

.300 

.40 

.90 

1.00 

.'250 

.50 

.80 

i.00 

Figure  A-4.  Example  of  Weapon  Base  External  File 

Damage  Expectancy  (DE) :  Damage  Expectancy  is  the  percen¬ 
tage  of  destruction  the  user  wants  to  inflict  on  a  particular 
target  class.  DE  can  be  input  from  a  file  or  interactively. 

If  there  have  been  any  target  class  additions  or  deletions  in 
the  target  editor,  the  model  will  only  accept  DE  input 
interactively.  The  user  can  enter  DEs  by  either  target 
category  (m,  v,  or  f)  or  by  individual  class  (Figure  A-5)  .  If 
the  DEs  are  entered  by  category,  all  the  target  classes  which 
are  classified  in  that  category  will  have  the  same  DE.  If  a  DE 
value  greater  than  or  equal  to  1  or  less  than  0  is  input  from 
a  file,  BRIK  will  either  reduce  it  to  . 999  or  raise  it  to  0, 
respectively.  After  DE  values  have  been  entered,  the  user  can 
check  the  DE  values  and  make  corrections.  Again,  BRIK  guides 
the  user  with  questions  and  prompts. 


Do  you  wont  to  enter  your  DE  data  by* 

1)  category  <v=Value,  i»=Military,  f=Force> 

2)  individual  closs. 

Select  one. 


Figure  A-5.  Example  of  DE  Menu 

Once  all  the  DE  values  have  been  input  and  edited,  BRIK 
will  compute  the  single  shot  probabilities  of  survival  <SSPS> . 
This  is  indicated  by  a  ‘PLEASE  STANDBY"  on  the  monitor  screen. 
The  single  shot  probabilities  of  survival  are  not  listed  on 
the  screen  and  cannot  be  externally  changed.  They  are  listed 
in  the  file  called  "SPARSE". 

When  the  calculations  are  complete,  BRIK  will  asK  the  user 
what  type  of  alert  rate  his  weapon  base  is  on.  This  will 
determine  what  percentage  of  the  target  base  to  use  for  the 
at  location. 

This  concludes  the  entry  of  all  required  parameters.  The 
weapon  availability  and  target  characteristics  along  with 
individual  DE  gcais  nave  been  specified.  The  next  two 
sections  describe  the  objective  functions  and  constraining 
rules  that  can  be  used  in  BRIK. 


Objective  Functions.  There  are  three  objective  -functions 
available  in  BRIK  (Figure  A -6>  .  The  first  objective  -function 
attempts  to  m»et  all  o-f  the  user's  goals  with  the  available 
arsenal.  If  al i  goals  can  be  met,  the  model  will  hold  the 
remaining  weapons  in  reserve.  The  arsenal  characteristics 
used  in  the  allocation  will  depend  on  the  'extreme  goal". 
Extreme  goal  options  will  be  discussed  in  the  Constraining 
Rules  section. 


This  section  builds  your  objective  function.  The 
function  is  divided  into  3  Major  categories.  For 
complete  explanation  of  your  options  please  consult 
the  user  guide.  Please  select  one  of  the  following. 
Enter  1,  2 ,  or  3. 

1)  This  problex  requires  allocation  restricted 
to  the  available  arsenal. 

2)  This  problem  requires  all  goals  to  be  met 
regardless  of  the  available  forces. 

3)  This  problem  requires  allocation  restricted 
to  the  available  arsenal.  Also,  the  DE  goals 
are  converted  to  upper  bounds.  If  this  pro¬ 
blem  is  picked,  the  hedging  options  must  be 
used  to  drive  the  allocation. 


FIGURE  A -6.  Objective  Function  Menu 


The  second  objective  function  forces  all  DE  goals  to  be 
met  regardless  of  the  weapon  availability.  It  can  be  used  to 
determine  the  type  and  number  of  weapons  to  add  to  an  existing 
arsenal,  or  to  create  a  new  arsenal  to  meet  DE  goals.  The 
arsenal  created  or  expanded  will  only  consist  of  those  weapon 


classes  input  in  the  weapon  class  input  section. 

The  third  objective  function  bounds  the  problem  with  both 
weapon  availability  and  D£  goals.  Since  the  DE  goals  are  not 
used  to  drive  the  allocation,  hedges  must  be  used  (see  next 
section).  The  Target  DEs  will  no  longer  be  treated  as  goals. 
They  set  the  upper  limit  on  the  amount  of  destruction  to  any 
target  class.  The  only  way  an  allocation  can  be  run  using 
this  objective  function  is  to  add  hedging  constraints  to  the 
problem  or  by  using  the  fifth  extreme  goal  and  allocating  the 
entire  arsenal  *  If  no  hedges  are  added  or  one  of  the  other 
extreme  goals  is  used,  BRIK  will  not  perform  any  weapon 
al  1  oc at i on . 

Constraining  Rules.  The  following  options  allow  the  user  to 
select  a  variety  of  rules  to  constrain  the  allocation.  There 
are  three  types  of  rules  available  in  ERIK:  an  extreme  goal, 
inappropriate  weapon/target  assignments,  and  several  types  of 
hedging. 

Extreme  Goal:  This  is  the  lowest  priority  goal  in  each 
objective  function.  Its  purpose  is  to  drive  the  allocation  to 
a  unique  solution.  Sometimes,  however,  the  rule  will  have  no 
effect.  For  example,  using  the  first  objective  function,  if 
all  the  weapons  were  used  in  the  allocation,  this  rule  will  be 
ignored.  BRIK  has  five  different  extreme  goals  (Figure  A-7) . 
The  first  four  goals  can  be  used  with  any  of  the  three 
objective  functions,  while  the  fifth  goal  is  only  available 
with  the  third  objective  function. 
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Th#  following  is  the  lowest  priority  goal  for  this 
allocation .  Please  select  one*  Input  1,  2,  3*  4, 
or  5* 

1)  Minimize  warheads  used* 

2)  Minimize  aeoatonnage  used* 

3)  Minimize  counternilitary  potential* 

4)  Hininize  total  equivalent  negatonnage* 

5)  Use  as  much  of  the  renaining  arsenal  as 
possible* 


Figure  A -7.  List  of  BRIK  Extreme  Goals 


Inappropriate  Weapon/Target  Interaction i  After  the 
extreme  goal  ha*  been  selected,  BRIK  will  ask  for  any 
inappropriate  weapon/target  interactions.  This  function  can 
be  used  if  there  are  some  weapon  classes  that  should  not 
attack  certain  target  classes,  either  for  military  or 
political  reasons.  An  example  of  a  time  urgent  target  is  a 
bomber  base  on  alert.  It  may  be  inappropri ate  for  a  weapon 
system  which  will  not  arrive  until  all  the  bombers  have 
launched  to  be  allocated  against  it.  Range  limitation  is 
another  reason  to  use  this  allocation  rule.  For  example, 
suppose  that  a  target  class  had  some  targets  out  of  range  of  a 
specific  weapon  class.  The  target  class  could  be  divided  into 
subsets.  One  set  would  incorporate  all  of  the  targets 
that  could  not  be  reached  by  the  weapon  class,  and  the  rule  to 
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prohibit  the  weapon  class  from  attacking  the  subset  would  be 
added  to  the  problem. 

Hedges!  Hedges  are  additional  constraints  which  are  added 
to  further  control  the  allocation.  There  are  seven  types  of 
hedges  available  (Figure  A-8)  .  After  a  hedge  is  created,  the 
model  will  always  return  to  the  menu.  Also,  once  created,  tho 
hedge  cannot  be  deleted  from  the  problem,  so  the  user  should 
be  cautious. 

The  first  two  hedge  types  are  similar,  both  enforce 
minimum  acceptable  levels  of  damage  on  certain  target  classes. 
However,  they  differ  in  the  weapon  set  used  to  meet  the 
constraint.  Hedge  1  insures  the  minimum  DE  level  is 
accomplished  using  ary  of  the  available  arsenal,  while  hedge  2 
attempts  to  meet  the  minimum  OE  goal  with  some  user-designated 
subset  of  weapons.  These  minimum  DE  goals  should  not  be 
confused  with  the  OE  goals  input  in  the  DE  data  section. 

Since  hedges  appear  at  a  higher  priority  in  the  first  two 
objective  functions,  the  minimum  goals  will  all  be  met  (or 
attempt  to  be  met)  prior  to  any  allocation  to  satisfy  the  DE 
goals.  Also,  since  objective  function  3  only  uses  the  DE 
goals  to  provide  an  upper  DE  limit  for  all  target  classes, 
these  hedges  are  needed  to  help  drive  the  allocation. 
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This  section  permits  selection  of  user  defined  goals* 
There  are  seven  types  to  choose  from*  Select  any 
particular  type  by  inputting  the  appropriate  number* 
fir  further  information  about  hedging »  see  the  users 
guide* 

You  have  space  to  select  20  hedging  constraints* 

The  following  is  a  list  of  hedging  types  available* 
Type  0  to  exit  this  section* 

1)  Enforce  a  minimum  level  of  damage  on  a  particu¬ 
lar  target  class* 

2)  Enforce  a  minimum  level  of  damage  on  a  particu¬ 
lar  target  class  using  a  specific  set  of  weapons* 

3)  Enforce  an  upper  level  of  damage  on  a  particu¬ 
lar  target  class  resulting  from  a  specific  set 
of  weapons* 

4)  Restrict  the  number  of  a  class  of  weapons  which 
can  be  allocated  to  a  target  clas3» 

5)  Build  your  own  constraint* 

6)  Enforce  a  minimum  level  of  damage  on  a  particu¬ 
lar  set  of  target  classes* 

7)  Restrict  the  number  of  weapons  which  can  be  al¬ 
located  to  each  target  in  a  particular  target 
class* 


Figure  A-8.  BRIK  Hedging  Options 


The  third  hedge  is  designed  to  limit  the  damage  on  a 
particular  target  class  resulting  from  a  specific  set  of 
weapons.  For  example,  if  it  is  desired  that  an  I CBM  target 
class  receive  at  most  fifty  percent  of  its  DE  from  bomber  and 

SLBM  weapon  classes,  this  hedge  should  be  used.  _ _ 

The  fourth  hedge  restricts  the  number  of  weapons  from  a 
particular  weapon  class  that  can  be  allocated  to  a  designated 
target  class.  By  restricting  the  amount  of  weapons  from  a 
designated  weapon  class  that  can  be  allocated  to  a  specific 
target  class,  BRIK  is  forced  to  use  other  weapons  to  attain 
the  target  class  DE  goal.  This  allows  the  remaining  weapons 


128 


of  th*  designated  class  to  attack  other  targets. 

The  fifth  hedge  allows  the  user  to  build  his  own 
constraint.  This  hedge  should  be  used  only  if  there  is  a  good 
working  knowledge  of  how  BRIK  builds  constraints.  This  hedge 
allows  the  user  to  create  hedges  that  are  not  indigenous  to 
BRIK,  such  as  a  common  carrier  constraint.  If  two  weapons  are 
carried  aboard  a  common  carrier  in  a  specified  ratio,  the 
weapons  must  be  used  in  that  ratio.  For  example,  if  equal 
numbers  of  two  weapon  classes  had  to  be  carried  on  a  bomber, 
the  model  should  allocate  exactly  the  same  amount  of  each 
weapon  type  in  the  problem.  The  following  constraint  forces 
an  equal  amount  of  each  weapon  class  to  be  used  to  attack  two 
target  classes!  j 


xll+  X1 2“  x21”  x22*  0 


The  first  two  terms,  X11  and  X12  *  indicate  the  amount  of  the 


first  weapon  class  allocated  to  target  classes  1  and  2.  The 


second  two  terms,  X21  and  X22*  indicate  the  amount  of  the 
second  weapon  class  allocated  to  target  class  1  and  2.  If 


both  weapons  are  used  equally,  the  difference  between  the 


number  of  each  used  must  be  0.  If  the  bomber  could  carry  2 
type  2  weapons  for  each  type  1  weapon,  the  constraint  will 
change  toi 


XH+  X12-  1/2X21-  1/2X22-  8 

This  forces  one  type  1  weapon  to  be  used  for  every  two  type  2 
weapons  used. 
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The  sixth  hedge  sets  a  minimum  level  of  acceptable  damage 
against  a  set  o-f  target  classes.  This  constraint  should  be 
used  carefully.  For  example,  if  the  user  wanted  at  least  30Z 
destruction  against  a  set  of  three  target  classes  and  one  of 
these  classes  was  extremely  easy  to  kill,  BRIK  could  attack 
the  easy-to-kill  target  class  and  leave  the  other  two  classes 
untouched.  It  is  suggested  that  this  hedge  be  used  only  with 
the  third  objective  function  because  all  of  the  target  classes 
have  upper  bounds  on  the  amount  of  damage  inflicted  on  them. 

The  seventh  hedge  restricts  the  number  of  weapons  that  can 
be  allocated  against  each  target  in  a  target  class.  This 
constraint  could  be  used  to  force  an  upper  limit  of  total 
weapons  that  can  be  used  against  an  individual  target  class. 
Also,  if  one  of  these  constraints  were  added  for  each  target 
class  a  targeting  strategy  of  allowing  at  most  one  weapon  per 
target  for  every  target  in  the  entire  base  could  be  enforced. 
This  hedge  could  be  used  to  force  an  integer  solution  if  the 
DE  goals  are  set  extremely  high  and  there  are  enough  weapons 
in  the  weapon  base  to  meet  this  hedge. 

After  the  final  hedge  is  input,  "PLEASE  STANDBY"  will 
appear  on  the  monitor  screen.  This  indicates  that  BRIK  is 
creating  the  matrices  and  files  necessary  to  run  the  problem. 
Once  the  files  are  built,  BRIK  will  begin  the  allocation 
algorithm.  This  is  indicated  by  prompts  which  indicate  the 
priority  currently  being  worked.  The  allocation  function  can 
take  anywhere  from  a  few  seconds  to  several  tens  of  minutes, 
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depending  on  problem  size  and  the  time  sharing  system  used  on 
the  computer. 

Output .  The  purpose  of  this  section  is  to  explain  the  output. 
The  following  example  uses  the  target  and  weapon  data  listed 
in  Figures  A-2  and  A-4  with  a  type  one  objective  function. 

BR1K  displays  the  output  data  on  the  monitor  screen 
<Figure  A-?>  and  also  places  the  output  into  an  external  file 
called  "alloc"  (Figure  A-10) .  The  screen  output  gives  the 
user  three  types  of  information!  the  total  value  of  all  the 
targets  destroyed,  the  percentage  of  each  target  class 
destroyed  including  a  detailed  list  of  the  type  and  number  of 
each  weapon  class  allocated  against  each  target  class,  and  the 
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Figure  A-?.  Screen  Output 
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number  of  weapons  of  each  weapon  class  used  in  the  allocation. 
“Alloc*  contains  two  sets  of  information.  The  first  set 
contains  data  from  the  goal  programming  package  and  the  second 
set  is  a  duplication  of  the  screen  output.  Set  one  consists 
of  a  five-column  table.  The  first  column  lists  the  subscript 
number  of  the  priority  or  of  the  variable  in  that  row.  The 
second  column  contains  the  objective  value  for  each  priority. 
The  third  column  contains  the  values  of  the  decision 
variables,  and  the  fourth  and  fifth  columns  list-  the  values 
for  the  positive  and  negative  deviational  variables  for  each 
goal  respectively  (Figure  A-18)  .  The  length  of  the  table'  is  a 
function  of  the  greater  number  of  either  decision  variables  or 
goals  in  the  problem. 

The  first  table  in  "Alloc"  contains  a  lot  of  information. 
For  example,  when  a  problem  uses  a  type  1  or  a  type  2 
objective  function,  if  the  number  in  column  2  for  row  2  is 
non-zero,  then  at  least  one  of  the  hedging  rules  could  not  be 
met.  In  the  "X  OPT"  column,  the  specific  weapon/target 
assignments  can  be  found.  To  compute  the  location  of  the 
weapon/target  assignment,  use  the  following  formula 
<j  ♦  <i-l)  x  number  of  target  classes) 
where  j  is  the  target  class  and  i  is  the  weapon  class.  For 
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RUN  NUMBER  1 


THE  OPTIMIZATION  ENDED  ON  SUBPROBLEM  6 

THERE  WERE  9  CONSTRAINTS  IN  THE  FINAL  OPTIMAL  TABLEAU. 


**mmmmmmmmm**mmm**mmm*mm**m***mx*xx**xmx 

O  OUTPUT  SUMMARY 
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A  OPT 
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POS  DEV 
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.0000 

4 

.0000 
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.0000 
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Figure  A-18.  Output  File  "Alloc“ 


•x ample,  in  Figure  ?,  there  are  five  target  classes  and  three 
weapon  classes  in  the  problem.  To  find  out  how  many  weapons 
of  class  2  <slbm)  were  assigned  to  target  class  4  (pol),  the 
formula  yields 

<4  ♦  <2~1>  x  5)  -  9 

Looking  in  column  aX  OPT"  at  row  nine,  yields  383.34,  which 
agrees  with  the  screen  output.  The  deviation  from  each 
constraint's  right  hand  side  is  listed  in  the  "POS  DEL**  and 
“NEG  DEV"  columns.  The  order  of  appearance  for  the 
constraints  in  the  goal  programming  problem  is  as  follows:  the 
target  OE  goals,  weapon  availability  goals,  extreme  goal,  and 
the  hedging  options.  Since  there  are  five  target  classes  used 
in  the  samp! 4  output  in  Figure  9,  the  rows  which  contain  the 
deviations  for  the  weapons  availability  are  4  through  8.  The 
amount  of  unused  type  2  weapons  is  474.8,  which  agrees  with 
the  screen  output. 

Although  all  of  the  information  is  readily  available  and 
easier  to  interpret  from  the  screen  output,  the  goal 

i 

programming  package  data  is  useful  to  give  the  advanced  user  a 
further  understanding  of  the  allocation.  The  ability  to 
eliminate  the  first  set  of  data  from  "alloc”  is  available  in 
BRIK's  sensitivity  section,  which  will  be  discussed  next. 

SENSITIVITY . RSRW 

This  section  discusses  the  sensitivity  capability  of  BRIK. 
Sensitivity  consists  of  changing  any  of  five  types  of 
parameters  {Figure  A-ll>.  After  the  limitations  of  this 
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section  are  explained,  the  use  of  this  routine  will  be 
described. The  sensitivity  analysis  used  in  this  section  should 
not  be  confused  with  a  linear  progamming  type  of  sensitivity. 
Instead  of  using  the  existing  final  tableau,  the  model  reruns 
the  problem  after  the  parameter  changes  have  been  made. 


MENU 

SENSITIVITY  RERUN 

1)  Change  DE  goals 

2)  Change  weapon  availability 

3)  Change  target  weights 

4)  Change  target  paraeeters 

5)  Change  weapon  parameters 

Enter  your  choice.  Type  0  if  you  are  finished. 


Figure  A-ll.  Sensitivity  Rerun  Menu 

Also,  this  section  cannot  change  the  objective  function  or  the 
hedges  used  in  the  original  problem,  if  the  user  needs  to 
change  the  hedges  or  the  objective  function,  a  new  problem 
must  be  run. 

DE  Goals.  This  subsection  is  used  to  change  the  DE  goals 
of  any  of  the  target  classes.  It  can  be  used  to  determine  how 
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a  change  in  DE  affects  the  allocation. 

Weapon  Availability.  This  option  can  be  used  to  see  how  a 
change  in  the  number  of  weapons  affects  the  allocation.  If  a 
problem  is  initially  run  with  a  weapon  availability  of  zero, 
increasing  the  number  of  weapons  has  the  effect  of  adding  a 
new  weapon  class  to  the  problem. 

Taroet  Class  Weights.  This  subsection  has  meaning  only  if 
the  first  objective  function  was  used  in  the  original  problem. 
If  either  of  the  other  two  objective  functions  were  used,  a 
change  in  weight  will  have  no  effect  on  the  problem.  The 
target  class  weight  is  a  function  of  target  value  and  target 
class  OE.  The  user  can  simply  input  a  new  weight,  or  can 
compute  a  new  weight  based  on  a  change  in  any  of  the  three 
parameters.  The  formula  is  as  follows: 

Weight*Target  Value  x  £<1  -  Target  DE> 3 

Target  Par ^meters.  This  subsection  allows  the  user  to 
change  two  parameters,  target  vul nerabi 1 i ty  and  target 
diameter.  A  change  in  either  of  these  two  parameters  will 
affect  the  survivabi 1 i ty  of  the  target  class  attached  by  each 
weapon  class. 

Weapon  Parameters.  This  subsection  allows  the  user  to 
change  three  weapon  parameters:  weapon  reliability,  weapon 
CEP,  and  weapon  yield.  A  change  in  any  of  these  parameters 
changes  the  effectiveness  of  the  particular  weapon  class  and 
will  affect  the  SSPSs  for  those  target  classes  which  the 
particular  weapon  class  can  attack. 
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After  the  user  is  -finished  making  a  parameter  change ,  BRIK 


returns  to  the  main  sensitivity  menu.  Mhen  all  changes  are 
complete,  the  user  is  given  the  option  to  suppress  the 
expanded  output.  The  expanded  output  consists  of  all  data 
that  is  entered  in  'Alloc'  that  is  not  displayed  on  the 
screen.  Output  suppression  is  recommended  if  the  problem  is 
large  and  it  is  anticipated  that  many  runs  will  be  made. 
Finally,  each  time  the  sensitivity  section  is  entered  all  of 
the  files  are  reinitialized  to  the  original  problem.  That  is, 
for  any  subsequent  sensitivity  run,  all  of  the  changes  made 
during  the  previous  runs  are  automatically  removed  before  the 
user  is  given  the  option  to  make  further  changes. 

This  concludes  the  user  guide.  It  was  the  authors' 
intention  to  provide  a  short  summary  to  make  BRIK  useable 
without  requiring  the  user  to  read  the  entire  thesis.  For 
further  detail  about  the  inner  workings  of  BRIK,  the  user  is 
referenced  to  the  thesis. 
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APPENDIX  B 


BRIK  Subroutine  and  Function  Listing 

ADDT6T  —  Enables  the  user  to  add  a  target  class  to  the 
problem 

ADDWPN  —  Enables  the  user  to  add  a  weapon  cl«»s  to  the 
probl cm 

AIJIN  —  Builds  the  -first  (ntgts  +  nwpns)  rows  ~ f  the 
Aij  matrix,  using  the  number  of  targets, 
number  of  weapon-**  and  the  SPARSE  matrix 
previously  built 

BOUT  —  "Main*  Program  that  calls  all  subroutines  that 
are  used  in  'he  allocation  algorithm 

GALMGT  —  Calculates  the  value  for  a  single  target  in  a 

force  target  class  if  the  user  has  selected  this 
option 

CHANGE  --  Permits  minor  changes  to  the  original  problem. 

After  the  changes  have  been  made  the  now  problem 
is  rerun  in  BRIKOUT.  The  following  list  of 
changes  is  allowedi 
Change  DE  goals 
Change  weapon  availability 
Change  target  weights 

Change  specific  weapon  or  target  parameters 

CINDX  —  Computes  the  relative  cost  coefficients  for  each 
variable  in  the  current  tableau  and  the  objective 
function  value  at  the  current  priority 

DECHEK  —  Displays  the  damage  expectancy  values  on  the 
screen 

DEFILE  —  Reads  Damage  Expectancy  <DE>  values  from 
an  external  file 

DEINAC  —  Allows  the  user  to  manually  input  desired 
DE  values  for  each  target  ejass 

DESAVE  —  Allows  the  user  to  save  to  an  external  file  his 
DE  val ues 

FILEIN  —  After  all  computations  are  finished  and  the 
needed  matrices  are  built,  this  subroutine 
stores  the  data  into  the  appropriate  files 
< SPARSE ,  AIJ,  PAGPIN)  for  usej  in  BOUT 
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FIX 


—  (Function)  Places  -floating  point  values  that  are 
within  l.E-5  o-f  an  integer  to.  that  integer 

HEADER  —  Prints  the  name  o-f  the  program  on  the  screen 
when  the  program  is  turned  on 

HEDGE  — -  Allows  the  user  to  select  -from  seven  user- 
defined  goals: 

En-force  a  minimum  level  o-f  target  class 
damage 

En-force  a  minimum  level  of  damage  on 
particular  target  class  frcm  a 
specific  set  of  weapons 
Enforce  an  upper  level  of  damage  on  a 
particular  target  class  frcm  a 
specific  set  of  weapons 
Restrict  the  number  of  a  specific  class  of 
weapons  which  can  be  allocated  to  a 
certain  target  class 
Build  a  user-defined  constraint 
Enforce  a  minimum  level  of  damage  on  a 
certain  set  of  target  classes 
Restrict  the  number  of  weapons  which  can 
allocated  to  each  target  in  a 
particular  target  class 

LI5TGT  —  Prints  a  list  of  target  class  names  on  the  screen 

LISTWP  —  Prints  a  list  of  weapon  class  names  on  the  screen 

OBJECT  —  Builds  the  objective  function  for  the  problem 
to  either  restrict  allocation  to  the  available 
arsenal  or  meeting  all  DE  goals  even  to  the 
extent  of  increasing  the  size  of  the  weapon 
arsenal.  This  subroutine  also  forces  the  user  to 
select  an  allocation  that  minimizes  one  of  the 
fol 1  owing: 

Warheads 

Megatonnage 

CMP 

EM? 

PAGE  —  Pages  the  display  to  the  screen 
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PARCH  1  —  Called  when  the  user  wants  to  change  one  of  the 
■following  target  paramet..'s: 

Name 

Vulr erabi 1 i ty 
Category 

Number  of  targets  in  the  class 

Diameter 

Priori ty 

Value 

PARCH2  —  Called  when  the  user  wants  to  change  one  of  the 
following  force  target  parameters! 

Reliability 

Cep 

Yield 

Warheads/weapon  system 

PERM  --  Performs  the  pivot  operation  on  the  pivot  element 

PHSE1  —  Reads  in  any  real  constraints  and  performs  a 
simplex  procedure  to  find  an  initial  basic 
feasible  solution 

PK  —  (Function)  Calculates  the  single  shot  probability 
of  survival  for  each  weapon/target  pairing 
if  PSI  expresses  the  target  vulnerability 

PLACE  —  Puts  the  objective  function  weights  for  the 

deviation  variables  at  the  current  priority  in 
the  correct  position  in  the  tableau 

POUT  —  Prepares  and  prints  the  solution  information 

READ!  —  Reads  in  the  goal  constraints  and  objective 
function  terms  assigned  to  priority  one 

READ2  —  Reads  in  the  goal  constraints  and  objective 

function  terms  assigned  to  the  current  priority 

TEST  —  Determines  the  next  entering  variable's  column  and 
raw. 

TGTEDT  —  Allows  the  user  to  edit  his  target  parameters 
and  to  add  or  delete  a  target  class 

TGTEXC  —  Allows  the  user  to  replace  a  target  class  with 
another  class 

TGTFIL  —  Writes  all  target  class  data  to  an  external 
formated  file 


TGTINF  —  Reads  t  .'get  data  from  an  external  user-defined 
formated  file 

TGTINS  —  Allows  the  user  to  interactive! y  enter  target 
class  data 

TGTP1  —  Displays  the  data  that  PARCH1  manipulates  on 
the  screen 

TGTP2  —  Displays  the  data  that  PARCH2  manipulates  on 
the  screen 

VTK  —  < Function)  Computes  single  shot  probabi 1 i ty  of 
survival  if  a  VNTK  number  expresses  target 
vulnerability 

WEIGHT  —  Allows  the  user  to  correct  or  change  the  value 
of  any  target  class 

WPNCH1  —  Allows  the  user  to  change  the  following  weapon 
class  data* 

Number  of  weapons 

Number  of  warheads/weapon 

Rel iabi 1 i ty 

CEP 

Yield 

WPNCH2  —  Allows  the  user  to  change  the  following  weapon 
class  datat 

Daily  alert  rate 
Generated  alert  rate 

WPNEDT  —  Allows  the  user  to  edit  his  weapon  parameters 
and  to  either  add  or  delete  a  weapon  class 

WPNEXC  —  Allows  the  user  to  replace  a  target  class  with 
another  class 

WPNFIL  —  Writes  all  weapon  class  data  to  an  external  file 

WPNINF  — ■  Reads  weapon  data  from  an  external  file 

WPNINS  —  Allows  the  user  to  interactively  enter  weapon 
class  data 

WPNP1  —  Displays  to  the  monitor  the  data  that  WPNCH1 
manipulates 

WPNP2  —  Displays  to  the  monitor  the  data  that  WPNCH2 
manipulates 


WTINTR  - 

ZEROIZ  - 


-  A1 lows  the  user 
assignments  are 


to  specify  which  weapon/target 
inappropriate 


-  Initializes  all  variables  and  matrices 
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APPENDIX  C 


BRIK  Variable  Listing 


TRIABLE 


DEFINITION 


achiev< j) 


ai  j<  j  ,  i) 

AIJ2 

alloc 

cr. 


Damage  Expectancy  achievement  of 
the  target  class 

Coefficients  of  the  aij  matrix 
for  row  j  and  column  it  RHS  is 
located  in  column  1 

File  for  aij  when  problem  is 
rerun 

File  containing  allocation 
output 

Crater  radius 


dl 


Lethal  diameter 


emt 

hedgl t 

iconnp< j ,n> 

i count 
IND<  I) 

ipcnt 


Equivalent  megatonnage 

Number  of  available  hedging 
constraints 

Subscript  of  the  jth  constraint 
at  priori ty  n 

The  ith  run  of  the  problem 

Indicator  row  that  marks  the 
eligibility  of  a  variable  to 
enter  the  basis 

Counter  used  to  page  the  screen 


iprin<i) 
isub< i , n> 


i type< i ,n> 


JCOL< i ,  1) 


Priority  of  the  ith  constraint 

Subscript  of  the  nth  constraint 
at  priority  i 

The  type  of  the  nth  deviational 
variable  at  priority  i, 

3,  Positive  deviational  variable 

4,  Negative  deviational  variable 

Type  of  variable  in  column  i 
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JCOLCi ,2) 


JRONCi  ,  1) 


JROW<i ,2) 


mtgts 


Subscript  of  the  variable  in 
column  i 

The  type  of  basic  variable  in 
row  i,  whore  type  is: 

X  =  2 ;  p  =  3 ;  n  =  4 

Subscript  of  the  basic  variable 
in  row  i 

The  amount  of  additional  target 
classes  that  can  be  added 


mwpns  The  amount  of  additional  weapon 

classes  that  can  be  added 


nc<n)  Number  of  constraints  assigned 

to  priority  n 

NCOLI  Number  of  columns  in  the  current 

tabl eau 


ncon( j ,n) 
nhedg 

node 

NPRIC 

npr  i  t 

nrcon 

NROWI 

ntgtpr < i) 

ntgts 
ntgtwh< i) 


Subscript  of  the  jth  constraints 
at  pr ior i ty  n 

Number  of  hedging  constraints 
used  in  the  current  problem 

Indicates  a  change  has  been 
made  to  the  target  class 
data 

Priority  currently  being 
optimized 

Total  number  of  priorities  in 
probl em 

Number  of  real  constraints 

Number  of  rows  in  the  current 
working  tableau 

Priority  of  target  class  i 
<1  -  7  ) 

Number  of  target  classes 

Number  of  warheads  for  a  force 
type  target  i 


/ 

/ 
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ntof<n)  Number  of  terms  in  the  objective 

-function  for  priority  n 

numcon  Number  of  constraints 

numtgt(i)  Number  of  targets  in  target 

class  i 

numwpn(J)  Number  of  weapon  systems  in 

weapon  class  j 

numwrh(i)  Number  of  warheads/weapon  system 

in  weapon  class  i 

nvar  Number  of  decision  variables 

nwpns  Number  of  weapon  classes 


pkl 


Correction  factor  applied  to  pk 
for  area  targets  <Diameter>0> 


rhsCi) 


Right  hand  side  of  constraint  i 


rl 

sparse( i , j> 


TA<N> 

TB<  i) 

temp 
TE< i  ,  j) 


tgtcat<i) 

i 


Lethal  radius 

Single  shot  probability  of 
survival  for  target  class  i 
against  weapon  class  J 

Total  deviation  from  the  goals 
at  priority  n 

Right  hand  side  constant  of  the 
constraint  in  row  i 

Countermi 1 i tary  potential 

Coefficient  of  the  variable  in 
column  J  of  the  constraint  in 
row  i 

Category  of  target  class  i 
v^value  target 
m*military  target 
f*force  target 

CEP  of  weapon  in  force  target 
class  i 


Desired  level  of  destruction  to 
target  class  i  <  V.  ) 
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tgtdia< i) 


Diameter-  <NM>  of  a  target  in 
class  i 


tgtnam< i) 
tgtpsi < i) 


Name  o-f  target  class  i 

Vulnerability  o-f  target  class  i 
designated  by  a  WNTK  number  or 
PSI  hardness  value 


tgtrel < i) 


Reliability  o-f  -force  target 
class  i 


tgtrmn< i) 


Amount  o-f  targets  remaining  in 
class  i  after  the  allocation 


tgtval ( i > 


Value  o-f  one  target  in  target 
class  i 


tgtyld(i) 


TL< I ,N> 


TT(N,J> 


valdes 


Yield  of  a  weapon  in  -force 
target  class  i 

Weight  assigned  to  the  basic 
variable  in  row  i  at  priority  n 

Probability  o-f  kill 

Weight  o-f  the  variable  in 
column  j  at  priority  n 

Total  value  destroyed  by  the 
al location 


wght< i , n) 


wl  e-f  t<  i) 


Weight  <Va1ue>  o-f  the  nth 
deviational  variable  at 
priority  i 

Amount  of  unused  weapons  in 
class  i  after  the  allocation 


wpname< j) 
wpncep( j> 


wpnddat j> 


wpnga< j) 


Name  of  weapon  class  j 

CEP  of  a  weapon  warhead  in 
weapon  class  j 

Daily  alert  rate  for  weapon 
class  j 

Generated  alert  rate  for  weapon 
class  j 


wpnrel < j) 


Reliability  (Probability  of 


14  6 


wpnrel < j) 

Reliability  (Probability  of 
penetration)  for  a  warhead  in 
weapon  c 1  ass  j 

wpnyl d< j) 

Mar head  yield  for  a  weapon  in 
class  j 

wused( i) 

Amount  of  weapons  in  class  i 
used  in  the  allocation 

XXX 

Dummy  variable  which  reads  the 
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Indicates  if  the  program  will 
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Flowchart  -for  BRIK  Allocation  Algorithm 
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The  following  is  the  extreme  goal.  It  is  the  lowest  priority 
goal. 
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****  FUNCTION  FIX  BRINGS  FLOATING  POINT  VALUES  THAT  ARE  WITHIN 
****  OF  AN  INTEGER  TO  THAT  INTEGER. 
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